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AUS BESPRECHUNGEN DER 1. AUFL. 


7 ,, Die rasche Entwicklung in der Erzeugung neuer Plaste und anderer Kunststoffe sowie deren Anwendung 
I in der Praxis laBt bestimmte Biicher, die sich damit beschaftigen, nach verhaltnismaBig kurzer Zeit 
i schon deshalb veralten, weil sie nicht den neuesten Stand der Technik bringen. Diese Lticke wird durch 
“fl das neu erschienene Druckwerk in erfreulicher Weise ausgefiillt. Es werden die Rohstoffbasis, die Er- 
zeugung, die Verarbeitung und Anwendung sowie die Eigenschaften der verschiedensten Plaste und 
sonstigen Kunststoffe leicht verstandlich und griindlich behandelt. Das Buch vermittelt nicht nur den 
Studenten und Laien wertvolle Kenntnisse von neuen synthetischen Werkstoffen, sondern gibt auch dem 
Fachmann gute Anregungen, und zwar dadurch, da8 auf Einzelheitenin der Verarbeitung und auf zahl- 
reiche Anwendungsmoéglichkeiten, die nicht jedem gelaufig sind, hingewiesen wird. In zahlreichen Tafeln 
wird auf die Eigenschaften hingewiesen, wobeiauch Werte genannt sind, diein den meistbenutzten DIN- 
Blattern nicht enthalten sind. Insbesondere ist es angenehm, da auf Vergleiche der verschiedenen Bau- 
stoffe ausfiihrlich eingegangen wird und auch die Anderungen der Eigenschaften bei verschiedenen 
Temperaturen behandelt werden.“ Fertigungstechnik 
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Uber die zwischen Abgriffen an einer von Wechselstrom 
radial durchfiuteten planparallelen Kreisplatte auftretende Spannung 


~ Von 
HAnns-HEINz WoLFF, Matawan, New Jersey, USA 
Mit 1 Textabbildung 
(Eingegangen am 3. Mai 1959) 


Ubersicht. Es wird die zwischen den freien Enden zweier Abgriffe auftretende Spannung berechnet, 
___ welche auf einer von Wechselstrom radial durchfluteten planparallelen Kreisplatte senkrecht angeordnet sind. 
_ Die Berechnung erstreckt sich auf den Fall der Kreisplatte, deren Stromriickfiihrung beidseitig der Platte, 
_ parallel zu ihr, radialsymmetrisch und auf beiden Seiten mit gleicher Phase und Amplitude erfolgt, sowie auf 
_ den Fall einer Kreisplatte, deren Stromriickfiihrung ausschlieBlich auf einer Seite der Platte parallel zu ihr 
und radialsymmetrisch verlauft. 


= 

as. 

Be * A. Einfiihrung 

Bs .. Die Erdung von Schaltungsteilen einer Hochfrequenz fiihrenden Anordnung muB8 oft aus 
2a raumlichen Griinden an der Grundplatte in der Nahe einer Hochfrequenzerdung vorgenommen 
~ _ werden, die einen starken Hochfrequenzstrom fiihrt. Der Hauteffekt in der Grundplatte sowie 
x das magnetische Feld eines solchen starken Erdungsstromes kénnen bekanntlich stérende 
> Spannungen tuber Erdungen in die Anordnung einfiihren, falls die Erdungspunkte und die 
5. Anordnung der Erdungsleitungen nicht sorgfaltigst ausgewahlt sind. Obwohl ein gemeinsamer 
__ _Erdungspunkt den Einflu8 des Hauteffektes vermeidet, so bleiben doch die induktiven Stérun- 
gen bestehen und kapazitive Stérungen nehmen hier im allgemeinen zu. 

<< : Eine Berechnung der durch Hauteffekt und magnetisches Feld hervorgerufenen St6érungen 
2 ist im allgemeinen wegen der raumlichen Gestalt und Lage der Schaltungsteile nicht médglich. 


ba 


Die vorliegende Arbeit behandelt deshalb einen idealisierten Fall, der es wenigstens erméglicht, 
_ die auftretenden Stérspannungen abzuschatzen. 


ae 


Ya 

In einer fritheren Arbeit1, im folgenden kurz mit [I] bezeichnet, war die Spannung » | 

% 

_berechnet worden, die sich zwischen zwei Punkten der Oberflache einer von Wechselstrom 

radial durchfluteten planparallelen Kreisplatte unter der Wirkung des Hauteffektes einstellt. 

_ Die Berechnungen bezogen sich dort auf eine Anordnung, bei der die Stromriickfiihrung beid- 

seitig, parallel zur Platte und radialsymmetrisch, auBerdem auf beiden Seiten nach Amplitude 

und Phase gleich vorausgesetzt wurde (vgl. ? Gl. (1.69), (I.69a) und (I.70)), und ferner auf eine 

Anordnung, bei der die Stromriickfithrung ausschlieBlich auf einer Seite der Platte parallel 

_* zu ihr und radialsymmetrisch erfolgt (vgl. Gl. (I.115), (I.115a), (1.116), (I.117), (I.117a) und 
= (1.118). 

Nun interessiert praktisch im allgemeinen wie oben ausgefiihrt nicht der Potentialabfall 

—  gwischen zwei Punkten der Oberflache des Leiters, wie er in [I] berechnet wurde, sondern die 
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# Spannung, die sich zwischen zwei Abgriffen einstellt, die sich in den Abstanden 7, und 7; 
vom Mittelpunkt der Platte, d. h. der Stromzufithrung befinden und die eine endliche Lange 
| aufweisen. Diese Spannung setzt sich, wenn sie unbelastet ist, aus dem bereits in [I] berech- 
___neten Potentialabfall zwischen den Abgriffsfu8punkten und der in der durch die Plattenober- 
s eee 

3 1 Wo.rr, H.-H.: Arch. Elektrotechn. Bd. 36 (1942) S. 493. 


2 Soweit wir im folgenden auf Gleichungen aus anderen Arbeiten z. B. [I] zuriickgreifen, werden wir diese 
mit dem Vorsatz des Kurzzeichens der Arbeit z. B. I. versehen. 
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flache, die beiden Abgriffe selbst und die geradlinige Verbindung ihrer freien Enden gebildeten 
Abgriffsschleife durch das magnetische Feld induzierten Spannung zusammen. 

In einer weiteren Arbeit!, auf die wir unten zuriickkommen werden und die wir im folgenden 
mit [II] bezeichnen wollen, waren auch bereits die magnetischen GréBen, insbesondere die 
auBeren Induktivitaten der beiden erwahnten Anordnungen berechnet worden (vgl. Gl. (II.24), 
(II.25) und (II.46)). 


B. Die Platte mit beidseitiger Stromriickfiihrung 


Die magnetische Feldstarke |{5,| im AuBenraum einer von Wechselstrom radial durch- 
fluteten planparallelen Kreisplatte (Raum zwischen Platte und Riickleitung), deren Strom- 
riickfiihrung beidseitig, parallel zur Platte und radialsymmetrisch, auBerdem auf beiden 
Seiten mit gleicher Amplitude und Phase erfolgt, betragt (vgl. Gl. (II.23), (I-50) und (1.53)) 


l= 2 = isin (ot +8), () 


4nmr 


worin y den Radius des Aufpunktes (Abstand von der Mittelsenkrechten der Platte), 7, den 
Gesamtstrom der Platte, w seine Kreisfrequenz und & die Voreilung des Stromes tay gegen die 
Stromdichte |t|,_, in der Mittelebene der Platte bezeichnen (vgl. Gl. (1.45)). 

Nun war in einer weiteren Arbeit?, im folgenden mit [IIT] bezeichnet, die Spannung 4,,, 
(vgl. Gl. (III. 41)) berechnet worden, welche in einer trapezférmigen Schleife induziert wird, 
die sich in einem magnetischen Feld der Form 

n= Hf 2 Igy sin (wt + &) (2) 


24 274 


(vgl. Gl. (III.1)) befindet. Die parallelen Seiten des Trapezes standen dort auf dem Radiusy 
und dem Vektor der magnetischen Feldstarke § senkrecht. Ihre Radien waren mit 7, (innerer 
Radius) und 7, (auBerer Radius), ihre Langen mit J; und /, bezeichnet. 

Wir wollen entsprechend annehmen, da die Abgriffe 
senkrecht auf dem Radius und dem Vektor der magnetischen 
Feldstarke ,, d. h. senkrecht auf der Oberflache der Platte 
angeordnet sind und ihre Dicke vernachlassigbar gegeniiber 
den iibrigen Dimensionen der Schleife ist. Wenn wir nun 
fiir die Abgriffe die gleichen Bezeichnungen anwenden, welche 
fiir die parallelen Trapezseiten in [III] benutzt wurden, 
namlich die Lange des einen Abgriffs, der sich in der Ent- 
fernung 7, vom Plattenmittelpunkt befindet, mit /,, die Lange 
des zweiten Abgriffs, der sich in der Entfernung 7; vom 
, Plattenmittelpunkt befindet, mit J; bezeichnen und schlieB- 
auftiinuce und Ateritesaf de Dlettencbene, lich den Winkel zwischen den beiden Radiusstrahlen 6 nennen 

(Bild 1), dann zeigt ein Vergleich von Gl. (1) mit Gl. (2), daB 
die in dem aus der geraden Verbindungslinie der der Platte abgewandten Enden der Abgriffe, 
den Abgriffen selbst und der Verbindungsgeraden ihrer Fu8punkte gebildeten Trapez indu- 
zierte Spannung w#,,, sich aus Gl. (III.41) zu 


mn toalt == Pil Gar tae 1,—l; Y 
EE rat BAY [V2 Fre Po 2 : a In + 1,—I, 
7 
1+tgy 1 
Is 
l,— 1, y, sin y Tt : 
i eee Bats ey eke : 
= KE arc cos ms y = Upp COS (wt + &) (3) 
pa it otgy 
7; sin® y 


ergibt, wobei der Winkel y durch Gl. (III.4) bestimmt ist. 


1 Wotrr, H.-H.: Arch. Elektrotechn. Bd. 37 (1943) S. 302. 
2 Wo trr, H.-H.: Arch. Elektrotechn. Bd. 44 (1959) S. 373. 
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lq 
Fiir die unbelastete Gesamtspannung u | _, die sich langs der geraden Verbindungslinie 
Ij 
der der Platte abgewandten Enden der Abgriffe einstellt, gilt daher, wenn wir die positive 
Stromrichtung wie in den Gl. (I.69) und (3) als zum Mittelpunkt gerichtet annehmen, nach 
dem Induktionsgesetz 


Yq 
=| + 4 | Ss inde = eRe at (4) 
< l; % 
z=+d/2 


und damit folgt durch Einsetzen von Gl. (1.69) und Gl. (3) in GI. (4) 
r i 


Vq . . 
i, 2% | Gof kd—coskd * 1; ty OO! Pia) 
v 
Ig —1 
4 19-8 1/2 Hoe l, a — I Aah es 
4 TU a a 1 
tery ——1 
vy; sin" y 
I,—l1 ; : 
— ; oat are cos 22 ya t,cos(@t+é). (5) 
“ 
ait tgy 4 
v? sin? y 


C. Die Platte mit einseitiger Stromriickfiihrung 
Fir den Fall einer radial durchfluteten planparallelen Kreisplatte, bei der die Stromriick- 
fuhrung ausschlieBlich auf einer Seite der Platte parallel zu ihr und radialsymmetrisch erfolgt, 
ist die magnetische Feldstarke , im AuBenraum zwischen Platte und Stromriickfithrung 


[vgl. Gl. (11.45), (I-97) und (I.99)| 


= Z ES: sin (wt + é), (6) 
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wobei I,, den Effektivwert des Plattenstromes und & seinen Voreilungswinkel gegen die 
Stromdichte |i], , in der dem feldfreien Raum zugewandten Oberflache der Platte bezeichnet 
(vgl. Gl. (I.93)). " 

Ein Vergleich mit Gl. (2) zeigt die Ubereinstimmung des magnetischen Feldes mit der in 
{III] behandelten Anordnung, so daB wir fiir die sich langs der geraden Verbindungslinie zwi- 


schen den der Oberflache der Platte abgewandten Enden der Abgriffe einstellende unbelastete 
Iq 
Spannung « | , wenn die Abgriffe auf der der Stromriickfiihrung zugewandten Oberflache 


i 
der Platte angebracht sind, aus den Gl. (1.115) und (III.41) 


[oases 


in Ty, sin (wt + pp) 


i 1 Coj 2k D—cosz2kD ¥; 
4 
, Bes, Bei [pepe fe Uf, 
+10 —//2 a 1 L,+ - _ im see 
| 1 tee |/ z a ——1 
vz sin® y 
SS yh arc cos #92” __,, +)lz, COS (tS) 4 (7) 
2 
—i——1 +otgy 
¥; sin® y 


erhalten. 
Fiir Abgriffe, welche auf der dem magnetisch-feldfreien Raum zugewandten Seite der 


Platte (z = 0) angebracht sind, ist ,,, = 0, so daB hier Gl. (1.117) gilt. 
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Es wird die langs der geraden Verbindung der freien Enden zweier senkrecht auf einer 
radial durchfluteten planparallelen Kreisplatte angebrachten Abgriffe auftretende unbelastete 
Spannung aus dem Spannungsabfail langs der Oberflache der planparallelen Kreisplatte und 
der in der durch die Abgriffe, die Verbindungsgerade ihrer FuBpunkte und die Verbindungs- 
gerade der freien Enden gebildeten Schleife induzierten Spannung ermittelt. Die Gleichartig- 
keit des magnetischen Feldes im AuBenraum einer radial durchfluteten planparallelen Kreis- 
platte mit dem magnetischen Feld eines geraden Leiters gestattet fiir die Berechnung der 
induzierten Spannung Formeln heranzuziehen, wie sie fiir den Fall einer in dem Feld eines 
geraden stromdurchflossenen Leiters befindlichen trapezfoérmigen Schleife, deren parallele 
Seiten parallel zum Primarleiter liegen, bereits friiher berechnet wurden. 

Die Berechnung wird durchgefiihrt fiir eine planparallele Kreisplatte, deren Stromriick- 
fiihrung beidseitig der Platte, parallel zu ihr, radialsymmetrisch und auf beiden Seiten mit 
gleicher Phase und Amplitude erfolgt, ferner fiir den Fall einer planparallelen Kreisplatte, 
deren Stromriickfiithrung ausschlieBlich auf einer Seite der Platte parallel zu ihr und radial- 
symmetrisch verlauft. Die Formeln unterscheiden sich in beiden Fallen lediglich durch den 
Faktor 2. 

Dr.-Ing. Hanns H. Wotrr, The W. L. Maxson Corp., 460 West 34th Street, New York 1, N. Y. 
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Theoretische Untersuchungen iiber den komplexen Scheinwiderstand 
magnetostriktiver Ring-Dehnungsschwinger* 


Von 
A. FRAUDE, Backnang 


(Mitteilung der Telefunken-GmbH., Geschaftsbereich Anlagen, Weitverkehr) 
Mit 5 Textabbildungen 
(Eingegangen am 22. Juni 1959) 


Inhaltstibersicht : 1. Einleitung.—2. Die Magnetostriktionsgleichungen.— 3. Elastizitatstheoretische Grund- 
_ lagen.— 4. Das Randwertproblem der Ring-Dehnungsschwingungen.— 5. Der komplexe Scheinwiderstand t.— 
6. Die Bestimmung der Eigenwerte «, und Vereinfachungen der Scheinwiderstandsgleichung. — 7. Ersatz- 
schaltbild und Ortskurve. — 8. Experimentelle Ergebnisse. — 9. Zusammenfassung und Schlu8. — Literatur. 


1. Einleitung 


Ebenso wie die Schwingquarze gehéren auch magnetostriktive Ferritschwinger zu jener 
Gruppe von Bauelementen der Nachrichtentechnik, mit denen es méglich ist, die im allgemeinen 
sehr hohe Giite und den relativ kleinen Temperaturkoeffizienten mechanischer Resonanzkreise 
unmittelbar fiir elektrische Schaltungen nutzbar zu machen. Die Voraussetzung dafiir sind 
die sogenannten Wandler-Eigenschaften des Werkstoffs, d.h. sein Vermégen, elektrische 
Energie in Form von elektrischen oder magnetischen Feldern in mechanische Energie (elasti- 
sche Schwingungen) umzuwandeln und umgekehrt. Ein solches System laBt sich als ein Wand- 
ler-Vierpol darstellen, der zwei elektrische EingangsgréBen (Strom und Spannung) in zwei 
mechanische Ausgangsgr6Ben (Kraft und Geschwindigkeit) transformiert. Bei Verwendung 
in elektrischen Schaltungen interessiert meist nur der zwischen den elektrischen Klemmen 
meBbare komplexe Schweinwiderstand und seine Abhangigkeit von der Frequenz, und zwar 
bei ,,Leerlauf‘‘ oder ,,KurzschluB“ auf der mechanischen Seite, je nachdem, welche Analogie 
zwischen elektrischen und mechanischen GréBen verwendet worden ist. Wenn man demnach 
den elektrischen Eingangs-Scheinwiderstand eines solchen Vierpols analytisch darstellen will, 
so ist dazu auch die Kenntnis des mechanischen Verhaltens erforderlich. 

Im folgenden soll die Rechnung fiir einen magnetostriktiven Ring-Dehnungsschwinger 
durchgefiihrt werden. Die ersten Untersuchungen dieser Art wurden von PIERCE [1] im Jahre 
1929 und — speziell fiir Ringschwinger — von BUTTERWORTH und SMITH [2] 1931 ausgefiihrt, 
jedoch unter sehr vereinfachenden Annahmen. Auch das von Srxtus [3] angegebene Ersatz- 
schaltbild und die darin auftretenden Gréfen sind nur als Naherungen fiir diinnwandige Ringe 
in der Grundfrequenz aufzufassen. Wenn diese Naherungen praktisch meist auch ausreichen, 
so mag die hier angestrebte Lésung fiir den allgemeineren Fall des Ringschwingers mit beliebi- 
gem Durchmesserverhaltnis fiir manche Untersuchungen doch von Interesse sein, zumal aus 
ihr auch die Verteilung der Oberschwingungen des Ringes hervorgeht. 

Auf die physikalischen Eigenschaften der magnetostriktiven Ferrite und ihre praktische 
Verwendung soll im Rahmen dieser Arbeit nicht naher eingegangen werden. Dariiber ist in der 
neueren Literatur bereits ausfiihrlich berichtet worden, wie z.B. von VAN DER Burct [4], [5], 
[6], DIETHELM [7], SCHWEIZERHOF [8] und anderen. 


2. Die Magnetostriktionsgleichungen 


Bei der Beschreibung des Magnetostriktionseffektes liegt der Gedanke nahe, sich an den 
piezoelektrischen Kristallen zu orientieren, obwohl zwar das eigentliche Analogon zur Magneto- 
striktion nicht die Piezoelektrizitat sondern der Elektrostriktionseffekt ist. Unter einer ge- 
wissen Einschrankung ist es jedoch méglich, die Gleichungen, die den linearen piezoelektrischen 


* Auszug aus der von der Fakultat fiir Elektrotechnik der TH Darmstadt genehmigten Dissertation (D17). 


eee ; 5 Archiv fir 
400 A. FraupEe: Komplexer Scheinwiderstand magnetostriktiver Ring-Dehnungsschwinger Diciccteetialt 
Ee I 


Effekt beschreiben, in analoger Form auch auf den nichtlinearen Magnetostriktionseffekt anzu- 
wenden. Diese Einschrankung besteht darin, daB man nur kleine Anderungen der magneti- 
schen und mechanischen ZustandsgréBen betrachtet, die um einen Ausgangszustand herum 
erfolgen, welcher durch eine gegebene Vormagnetisierung bzw. mechanische Vorspannung 
erzeugt wird. Damit werden einerseits die Beziehungen zwischen den einzelnen Zustands- 
groBen linearisiert, zum anderen erhalt man dann die VerhAltnisse, wie sie beim magnetostrikti- 
ven Schwingungsvorgang tatsdchlich vorliegen. 

Dieser erfordert namlich stets eine bestimmte Vormagnetisierung, die entweder mit Hilfe 
eines konstanten Gleichstromes durch die Wicklung des Schwingers aufrechterhalten wird, 
oder durch die remanente Magnetisierung gegeben ist. Dieser Vormagnetisierung iiberlagert 
sich ein durch einen Wechselstrom erzeugtes, hinreichend kleines Wechselfeld zur Anregung 
der Schwingungen. 

Auf dem Gebiete der Piezoelektrizitat leitet man nun die zur Beschreibung des Effektes 
erforderlichen Gleichungen mit Hilfe einer Energiebetrachtung aus den Funktionen fiir die 
freie Energie, den thermodynamischen Potentialen, ab!. In analoger Weise kann man beim 
Magnetostriktionseffekt vorgehen. Die in dem Material pro Volumeneinheit enthaltene Energie 
setzt sich aus drei Anteilen zusammen, wenn man — bei adiabatischer Zustandsanderung — 
thermische Glieder in dem Ausdruck fiir die freie Energie vernachlassigt: rein elastische - 
Energie’, rein magnetische Energie und Magnetostriktionsenergie. Schematisch und ohne In- 
dizierung geschrieben, nimmt das thermodynamische Potential die Form an 


=i 58+ =p, H+ d-0-H. (1) 


Darin sind o [kp/cm?] die mechanische Spannung, H [A/cm] die magnetische Feldstarke, 
E,[kp/cm?] der Elastizitatsmodul bei konstanter Feldstarke H, uw, [Hy/cm] die reversible 
Permeabilitat bei konstanter mechanischer Spannung o und schlieBlich d [Gau8 cm?/kp] eine 
Magnetostriktionskonstante. 

Ausfiihrlich geschrieben wiirde die Gl. (1) in den drei Termen Produkte aus den sechs 
méglichen Spannungskomponenten und aus den drei Feldstarkekomponenten enthalten, mit 
einer entsprechenden Vielzahl von Koeffizienten, die sich jedoch z. T. reduziert, wenn man das 
Material als isotrop annimmt; ein isotropes festes Medium hat nur zwei Elastizitatskonstanten 
und eine Permeabilitatskonstante. Eine weitere Vereinfachung ergibt sich durch Beschrankung 
auf nur eine Koordinatenrichtung. 

Aus dem thermodynamischen Potential als Funktion der mechanischen Spannungen und 
der magnetischen Feldstarke ergeben sich durch partielle Differentiation nach den Spannungs- 
komponenten die Komponenten der relativen Dehnung e, und durch Differentiation nach den 
Feldstarkekomponenten die Komponenten der magnetischen Induktion B [GauB]. Es ist also, 
bei Beriicksichtigung wieder nur einer Koordinatenrichtung, 


SLA Tee . | 3 
pene En ot+d-H, (2a) 
oC — =— ° . 
ee ee (2b) 
und ferner 
ee e a fe mc 
60 0H ns nee Pa ele ages (2 c) 


Das ist eines der vier ineinander iiberfiithrbaren Gleichungssysteme zur Darstellung des 
Magnetostriktionseffektes (s. [4]). Die erste Gleichung beschreibt dabei den direkten, die zweite 


 Siehe dazu die ausfiihrliche Darstellung bei Capy [9] oder die kiirzere und tibersichtlichere Zusammen- 
stellung in den ,,Standards on Piezoelectric Crystals‘ [10a u. c.]. 
 Dabei ist die kinetische Energie nicht beriicksichtigt; die Gleichungen gelten zunachst also nur fiir stati- 


sche oder quasistatische Zustandsanderungen bei Frequenzen weit unterhalb der Resonanzfrequenz des be- 
treffenden Systems. 
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Wahrend nun im Falle der Piezoelektrizitat die entsprechenden Gleichungen in weiten 
Grenzen giiltig sind, trifft dies bei den Magnetostriktionsgleichungen nicht zu, und zwar nicht 
einmal fiir den Zusammenhang zwischen «¢ und o, wo infolge des JE-Effektes Abweichungen 
vom linearen Hooxeschen Gesetz zu beobachten sind. Die Koeffizienten des Gleichungs- 
systems (2a, b) sind keine Konstanten, sondern hangen wesentlich von der Vormagnetisierung 
ab (s. [4], [5], [6]). AuBerdem sind sie wegen der elastischen und magnetischen Dampfungs- 
verluste als komplexe GréBen anzusetzen; der komplexe E-Modul entspricht dabei einer ge- 
schwindigkeitsproportionalen elastischen Dampfung. 

Die einzelnen Kristallite des Ferrits zeigen an sich zwar anisotropes Verhalten, jedoch liegen 
die Ferrite als Sinterwerkstoffe niemals als Einkristalle vor, sondern haben stets ein poly- 
kristallines Gefiige. Infolge der statistischen Verteilung der einzelnen Kristalle und ihrer 
Achsen auf alle raumlichen Richtungen innerhalb eines Ferritkérpers tritt die Anisotropie 
auBerlich nicht in Erscheinung, der Werkstoff ist quasiisotrop. Dieses gilt jedoch nicht fiir 
ein vormagnetisiertes Material. In diesem Zustand ist ein ferromagnetischer Werkstoff keines- 
wegs mehr isotrop. Die beiden Koeffizienten, die am starksten von der Vormagnetisierung 
beeinflu8t werden, namlich die reversible Permeabilitat uw, (oder yu.) und die Magnetostriktions- 
konstante d, interessieren im Verlauf der Rechnung nur in einer Koordinatenrichtung; anderer- 
seits andern sich E-Modul und Querkontraktionszahl y nur relativ wenig mit dem Magneti- 
sierungszustand. Deshalb wird im folgenden trotz der notwendigerweise vorhandenen Vor- 
magnetisierung das Material als isotrop angenommen. 


3. Elastizitatstheoretische Grundlagen 


Mit den Gln. (2a, b) ist der Magnetostriktionseffekt an sich hinreichend beschrieben. Diese 
Gleichungen allein geniigen jedoch noch nicht, um die Differentialgleichung fiir die elastischen 
Schwingungen aufzustellen, zu denen ein gegebenes System durch Ausnutzung des Magneto- 
striktionseffektes angeregt werden kann; dazu sind einige Beziehungen erforderlich, die von 
der Elastizitatstheorie geliefert werden (s. SOMMERFELD [11)). 

Die elastischen Eigenschwingungen eines isotropen kraftefreien Festkérpers mit den 
beiden LAmEschen Materialkonstanten uw und J [kp/cm*] werden durch die vektorielle Differen- 
tialgleichung beschrieben 

vine (2 +A) grad div u —protrotu. (3) 
In dieser Gleichung ist u der Verriickungsvektor [cm], 9 die Dichte [gr/cm*] und ¢ die Zeit. 
In Zylinderkoordinaten 7, g, z gelten fiir den Zusammenhang zwischen den Spannungs- und 
Dehnungskomponenten und den drei Komponenten des Verriickungsvektors u = (u,, ug, 4,) 
die Relationen: 


6, = 2 me, -Adivu = 2p 4a(+ 2 (r u,) +e , (4a) 
tp = 22, bAdivu=2n(— S24“) 44(— 5 ru) + e+ Se), (4b) 
o,= 2 me, bAdivi = 2% 4 a(- 2 vu) +4 a+), (40) 
com anrn =e (e ete 8) a 
To, = 2M Vos = i += a). (4e) 


Ou, , OU, 
Tay = 2M Yay = (et HY. (4f) 
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Ferner miissen die einzelnen Spannungskomponenten in den drei Koordinatenrichtungen die 
inneren Gleichgewichtsbedingungen. erfiillen: 


. Y 7) th 

in r-Richtung wr 4 8 oF <P as pe petcy EEN OS (5a) 
if 

: . iG] 

in y-Richtung : e Sl “ee += aot ns =O, (5b) 
yr ap 

in z-Richtung a -- 7 a - ; oa 07 (5¢) 

2 


In Gl. (3) sind die Gln. (4) und (5) zusammengefaBt enthalten. Wenn der betreffende Kérper 
frei von duBeren Kraften ist, muB ferner seine Oberflache spannungsfrei sein; diese Rand- 
bedingung ergibt zusammen mit der DGI. (3) ein Randwertproblem, das sich jedoch nur in 
wenigen Sonderfallen streng lésen 1a8t. Dazu gehdren z.B. die Torsionsschwingungen von 
Staben oder Hohlzylindern mit kreisférmigem Querschnitt. Dagegen sind die Dehnungs- 
schwingungen solcher Kérper nicht streng berechenbar, weil in diesem Falle durch die For- 
derung der Spannungsfreiheit an den Oberflachen das Randwertproblem tiberbestimmt ist. 
Hierauf aber fiihrt das Problem des Ring-Dehnungsschwingers, wenn auch die achsialen 
Dehnungsschwingungen beriicksichtigt werden, die der Ring infolge der Querkontraktion aus- 
fiihrt. 

Nimmt man jedoch die Hohe des Ringes als hinreichend klein an, so kann man seine radialen 
Schwingungen mit Hilfe der Beziehungen des ebenen Spannungszustandes berechnen. Dabei 
wird die Querkontraktion in achsialer Richtung vernachlassigt und nur in der Ebene des Ringes 
beriicksichtigt. Der ebene Spannungszustand ist dadurch gekennzeichnet, daB alle Spannungs- 
komponenten in einer Ebene liegen (hier in der v-y-Ebene) und sich in der Richtung senkrecht 
dazu nicht andern (s. Féppr [12]). Es treten nur die Spannungskomponenten o,, o, und Te 
auf, die die beiden Gleichgewichtsbedingungen erfiillen miissen: 


1 OTyo 


in v-Richtung = Mae ill as 0; (6a) 
y i v Op 
in y-Richtung : “ fie a ce ee: (6b) 
vy Op v Y 


Wenn man statt der LAmEschen Konstanten die gebrauchlicheren Materialkonstanten Elasti- 
zitatsmodul F und Querkontraktionszahl » verwendet, so folgt aus den Gln. (4a—f) 


1a 

Op eee pct Wi eate (7a) 
E 

On = ae Sp Pes (7b) 
12, 

Tre = rete “Vro* (7c) 


4. Das Randwertproblem der Ring-Dehnungsschwingungen 


In Bild 4 ist ein Ferritring mit reckteckigem Querschnitt skizziert, der in einer zweiteiligen 
Schale nahezu dampfungsfrei gelagert ist, so daB er frei schwingen kann. Die Schale tragt 
auBen die Wicklung zur Erzeugung des die Schwingungen anregenden magnetischen Feldes. 
Dieses hat — ebenso wie die Vormagnetisierung — nur die y-Komponente, H,=w J Gta mY 
(J = Wechselstromamplitude, m = Kreisfrequenz der Anregung, w = Windungszahl). Damit 
liegen aber auch die Induktion B, sowie die durch das Feld hervorgerufenen Spannungs- und 
Dehnungskomponenten parallel zur Ringebene. AuBerdem sind die genannten GréBen in 
Richtung der Achse des Ringes konstant, so daB — zumindest in einem gewissen Frequenz- 
bereich — tatsdchlich der Fall des ebenen Spannungszustandes vorliegt. Dieser verlangt an 
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sich keine Beschrankung hinsichtlich der Hohe h des Ringes, wenn nur die genannten Voraus- 
setzungen erfiillt sind. Dennoch wird im folgenden die Hohe als klein gegeniiber dem AuBen- 
durchmesser angenommen (s. a. Abschn. 8). 


e Halterungsschale 
5 ee [mit Wicklung 


ea [oN 
SY Mes 


Bild 1. Querschnittsskizze eines Ferrit-Ringschwingers mit Schale und Wicklung. 


Wegen der Kreissymmetrie des Ringes und der Anregung (die Wicklung sei gleichmaBig 
_ tuber den Umfang verteilt) sind = = 0, 4 = 0 und auch t,, = 0, so daB nur die beiden 


Normalspannungen o, und o, und die zugehérigen Dehnungen e, und e, verbleiben. 

Nach Gl. (2a) setzt sich die gesamte Dehnung «,,,, aus einem rein elastischen und einem 
magnetostriktiven Anteil infolge des Feldes H, zusammen, wahrend die Dehnung e, von H, 
unbeeinfluBt bleibt: 


Uy 
fogs =p t d+ Hy=—, (9) 
ou, 
€, oe aye . (10) 


Nach dem d’ALEMBERTschen Prinzip 1aBt sich ein dynamisches Problem auf ein statisches 
zuriickfiihren, wenn in den Gleichgewichtsbedingungen die Tragheitskrafte beriicksichtigt 
werden, wie das auch in Gl. (3) geschieht. Die Bewegung jedes einzelnen Volumenelementes 
erfolgt hier nur in radialer Richtung. Somit folgt aus den Gln. (6a, b) — nach Wegkiirzen 
von 7- dr- dy — 


— 2 
ee (11) 


Setzt man nun die Gln. (9) und (10) in die Gln. (7a, b) ein, so erhalt man mit Gl. (11) unmittel- 
bar die Schwingungs-DGl. 
au 1 Bee oe 48) Q eu q. Holt) dw fet (12) 


or? ze Y OY ry EEC os Y 20 7 


(Statt wu, wird im folgenden stets nur noch wu gesetzt.) 

Da jetzt nur der eingeschwungene Zustand betrachtet werden soll, gilt fiir w die gleiche 
Zeitabhangigkeit wie fiir die Anregung, u(r, ¢) = u(r) e°*. Aus der partiellen DGI. (12) wird 
damit die gewohnliche inhomogene BEssELsche DGl. 


ou 1 Ou a \2 1 ee 1 
or” - PU a3 (=) ste | u s Pp (13) 
wenn gleichzeitig zur Abktirzung 
(14 — v?) = O20 (14) 
und . 
eee (15) 
27 


geschrieben wird, wobei « wegen des komplexen E-Moduls ebenfalls komplex ist. 7,, = = (7, +7) 
ist der mittlere Radius des Ringes. 

Zu der DGl. (13) treten zwei Randbedingungen, die sich aus Gleichgewichtsforderungen 
an der inneren und auBeren Begrenzungsflache des Ringes ergeben. Von den beiden Spannungs- 
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komponenten o, und o, muB fiir 7 = 7, und y = 7, (Innen- und AuBenradius) die radial ge- 
richtete Komponente a, zu Null werden, ‘da diese Flachen frei von auBeren Normalkraften 
sind'. Man gelangt so zu der inhomogenen Randwertaufgabe 


ye ae shy = 4 (j\u=—KS, (16a) 
y Y es y 
ie Pee Ree fire eee (16 b) 
iff Y Y= 7, 


(’ = Ableitung nach 7). 
Nach der Theorie der Randwertprobleme (s. INcE [13] und KAMKE [14]) hat die Aufgabe 


Lin, th) i Ob = PAT) scare oe mt a7) 
U;(u) = ¥; meer * 
die Lésung : 
u(r) = J G(r, &, a) + (6) do + 27; G(r) . (18) 


4 


Hierin ist G(r, ,«) die als Hilfsmittel verwendete GREENsche Funktion, die ihrerseits eine 
Lésung der homogenen Randwertaufgabe 


EG areas jaz... 
U(u) ='0 

darstellt, jedoch nur fiir den Fall, daB diese Aufgabe keine eigentliche Lésung besitzt, also fiir 

«~,, wenn mit «, die Eigenwerte der homogenen Aufgabe bezeichnet werden. Die Werte 

« =a, sind Polstellen der GREENschen Funktion. Die au8erdem in der Gl. (18) auftretenden 

Funktionen G,(v) sind Lésungen der Randwertaufgaben 


L(G) +o2G;=0, 
UG) =15' UG) =0 pt, w=1,2,...m 
Man kann nun die GREENsche Funktion in geschlossener Form aus der Grundlésung der DGl. 


ableiten oder aber die Darstellung in Form einer Eigenfunktionen-Reihe wahlen, was im vor- 
liegenden Falle vorteilhafter ist: 


(19) 


(20) 


co 


G(r, eo x) on >» 1 : Uy(7) + Vy(€) : (21) 


2 
a*—gh 1% 


fase S Un(§) + On(€) dé 


In dieser Gleichung sind die w,,(7) die Eigenfunktionen der Aufgabe Gl. (19) wahrend die », (7) 
diejenigen der zu der Randwertaufgabe (19) adjungierten Aufgabe 


TD) ets On : : 
7=1,2,...M 


Vio (22) 


darstellen. Diese Eigenfunktionen sind nun aufzufinden. 
In dem Ansatz mit der BEssELschen und NEumANNschen Funktion 1. Ordnung 
Pas Y 
u(r) = ¢ Ja(a #) + Cy N,(o 2 
Ym, Ym, 
kann von den beiden Konstanten aus den homogenen Randbedingungen Gl. (16b) nur eine 
bestimmt werden, die andere bleibt unbestimmt. Die Randbedingungen sind dann nur fiir 
a = a, erfiillbar, also eben fiir die Eigenwerte des Systems. Daraus folgt: 
C i(a2) + ca (a7) nl), 
Ym Ym 
Cy i(«22) + Cy n{a.?2) ON 
Yn, Ym 
* Das entspricht bei dem einleitend erwahnten Wandler-Vierpol einem ,,Leerlauf auf der (mechanischen) 


Ausgangsseite; zwischen mechanischen und elektrischen GréBen gilt hier die Analogie: Strom + Kraft, Span- 
nung «+ Geschwindigkeit. 
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wobei 7 (« *) und na “| die im folgenden haufig gebrauchten Abkiirzungen bedeuten 


Ym Ym 
yea ; J (2 4 
i(«) == Iolo) (1 — ») SaeraaTe (23a) 
Lo 
Ym 
ne: ee 
n(x 4 a No(s =| ( ia ”) Be es °) 
Ym 


= 24) nfo) A(x,) =0, (24) 


und fiir die Eigenfunktionen selbst ergibt sich bis auf einen konstanten Faktor der Ausdruck 


40) = FilonX)— 7 Milan). (25) 


Die zu Gl. (19) adjungierte Randwertaufgabe lautet ausfiihrlich geschrieben 


2 
v— i +(Shu=o, (26 a) 
v—(1—»)—=0 fir“. #, (26 b) 
if Vey, 


Sie hat als Lésung die Eigenfunktionen 


(7) 
4|&, — 
: )-——"™ ; N,| 


UA = a We Sea eae as 
m2 


mit den gleichen Eigenwerten «,, wie sie sich fiir die urspriingliche Aufgabe aus Gl. (24) er- 
geben. Setzt man nun die Gln. (25) und (27) in (21) ein, so ist die GREENsche Funktion bis 
auf das im Nenner noch aufzulésende Integral bestimmt. 
In der Lésungsgleichung (18) fehlen nun noch die beiden Funktionen G,(v) und G,(r) als 
Lésungen der Randwertaufgaben 
a \2 1 
Bic) 


o fir r= 7, 
ie ee 


ag) | = 7 ty( (27) 


1 


u a ' a \2 u al , 
G+26,+|(-] 7]6.=0, G44 


ees (28) 
ea es uk G+ = 


Omit ey 


Sie lassen sich ohne Schwierigkeiten angeben, und es ist 
rK() Gp G,) = es gees n(a “| a n(a 2 Tio”) 
v% Yq me V5 Ya . A(a) 'n Vn Ym 'mn Vin 


dame gralany agi ety byt 
— (122 (a2) oo (a) Nl )) (29) 


ZweckmaBigerweise entwickelt man auch diesen Ausdruck in eine Reihe nach den Eigen- 
funktionen u,(7) 


¥m 


y K(} G+= G,| z 54: u,(r) « (30) 
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Die Koeffizienten a, in dieser Reihe ergeben sich durch Multiplikation beider Seiten mit 
y-u,,(r) und nachfolgender Integration mit den Grenzen 7; bis 7,. Wegen der Orthogonalitat 
der Eigenfunktionen ist 


"a =o firnA~Am 
tg (0) + ty() 9 dt ee (31) 
r; SAO LUNI) are 
so daB sich die a, aus der Gleichung bestimmen 
J - Ge Gy gl”) dr 
V; Yq 
a,=yK.- 5 (32) 
fi -vdr 
Damit ist dann 
ile G, + = 6,| u,,(v) -¥ dr 
co "% ‘q 
it % : 
rK[LG4+26)—%-K. 2 a) 5) 
ie fae vy+dy 


Diese Gleichung wird nun zusammen mit der Reihe Gl. (21) fiir die GREENsche Funktion in die 
Lésungsgleichung (18) eingesetzt, so daB sich als Lésung der Randwertaufgabe Gl. (16) schrei- 
ben 1aBt: 


cates J un(8) : i Pp G(r) +7 640] un (t) 7+ dr 
WN So KD) eee gg ae re u() (4) 
{ ug(é) » dé J lan ved 


Mit der obigen Schreibweise ist bei der Integration der GREENschen Funktion bereits die 
gliedweise Integrierbarkeit angedeutet. Das setzt jedoch die Konvergenz dieser Reihe voraus. 
Sie 1aBt sich bei den beiden Reihen Gl. (21) und (33) leicht nachweisen, wenn man in den héheren 
Gliedern dieser Reihen statt der BEssELschen- und NEuMANNschen Funktionen deren asymptoti- 
sche Naherungen nimmt und die dann entstehenden Ausdriicke nach geringer vereinfachender 


Umformung mit den Gliedern der Reihe = als Majorante vergleicht. Uber die Konvergenz 
dieser Reihe s. KNopp [15]. 


5. Der komplexe Scheinwiderstand f 


Mit Gl. (34) ist — bis auf die Auswertung der Integrale — die Deformation des Ringes bei 
seinen Schwingungen als Funktion des Radius 7 beschrieben. Die zeitliche Abhangigkeit wird 
durch Multiplikation mit exp (jw?) beriicksichtigt. Mit den Gleichungen (7a, b), (9) und (10) 
erhalt man ferner Beziehungen fiir die mechanischen Dehnungen und Spannungen. Von den 
letzteren wird im weiteren Verlauf nur die Komponente o, gebraucht, da fiir die in der Wicklung 
des Schwingers induzierte Spannung Ul nur der magnetische FluB in y-Richtung maBgebend 
ist. Die nach Gl. (2b) magnetostriktiv erzeugte Induktionskomponente B, = d, - a, liefert zu 
diesem keinen Beitrag. Fir die somit allein interessierende Komponente B, gilt demnach mit 
den Gln. (2b), (9), (10) und (7b) 


&| by 


Sygiok 
eho (35) 
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Aus der Integration tiber den Querschnitt des Ringes erhalt man sodann fiir den magnetischen 
FluB8 in gy-Richtung 


f B, dr (36) 


Yr 


®, =f B,-dK =h. 
P 4 
mit df = h- dy (h = Hohe des Ringes). Damit ist dann 


o 4%: Enh [ Mar +» (u(r,) — 0) Jett (pp — SH) O°F* J 2. oie! 


@ 1— Y 1— 2 20 1; 


% 


In diese Gleichung wird nun der mit Gl. (34) gewonnene Ausdruck fiir w(r) eingesetzt: 


a a 
uy(é) + dé i LF Gn + = G,(r) bi, 0) ard 
a @.Ey-hew- Jel? yor? t; : a? — of @ or 

? — yp?) 2 1 2 Co— op '% am b hs ‘<— hetesS "a 

up(é) +E dé un(r) ¥ dy 
"% "% 

° a a _ @Ey PU fone i Ie g 

és J My(7) y ae: (14 ("a) za) ee c | ue 2m v; 


t 


Eine erste Vereinfachung dieser Gleichung erreicht man durch die Auswertung des Zahler- 
Integrals in Gl. (32) 


a2 


Ao [ [EG +E Gn) 409) <7 ar = gle.) — slr) 


Ym 4 


unter Verwendung der Formel 
(0? — BP) [Zp (0 x) Z,(B x) x dx = B xZ,(ax) Zp_1(B x) —axZy_,(ax)Z,(Bx), (38) 


wobei Z,(«x) und Z,(B x) zwei beliebige Zylinderfunktionen darstellen. Dabei erhalt man 
zunachst einen recht umfangreichen Ausdruck; diesen kann man jedoch mit 


= v; F z 
4 (+ 4) i (+n “| 
1 ml Ym 


(39) 


ae Ya 
Tee, Ni ey, —— 

Ym Ym 
wegen A(«,) =O nach Gl. (24) und mit Beriicksichtigung der WRonskIschen Determinante 


Bl Ta(2), N,(2)} = [2 Nal) 2 (40) 


Jz) Ny@)| 

relativ leicht vereinfachen. Ferner wird in den yon der GREENschen Funktion herriihrenden 
Integralen statt der Integrationsvariablen & nun ebenfalls ry gesetzt. Damit nimmt die 
Gleichung fiir ®, die einfachere Form an 


te 2 
) / Uy(7) - = dr +4 (uy (7a) as wi) 


@Eyew-h- Je)! = v2, 1; 
es Te 2S 


4 


2 (1i—v?)-22 — a2 — ae ", 
fairer 
" 
a Ex weJrheeit x, 
a be 5| = ] (41) 
Hg (t—¥?) 20 "5 
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Das Integral im Nenner des Bruches 14Bt sich mit Hilfe der bekannten F ormel 


| Z3(ox) x dx =~ [23 (0x) —22 Z,(a%) Zps(07) + se) (42) 


auflésen. Nach einfacher Zwischenrechnung erhalt man — wieder mit Gl. (39) und (40) — 


f [ral etal cate | 


(43) 


9 1 
u,(r) 1 dr = - 5 
eRe oe 
1; : Ym ‘mn 
Fir die Lésung des in Gl. (41) nun noch verbleibenden Integrals ist die Kenntnis zweier 
Funktionen erforderlich, die mit den Zylinderfunktionen eng verwandt sind; sie sind Lésungen 


von inhomogenen BeEssELschen Differentialgleichungen (s. WATSON [16)]). 
Die Lommetsche Funktion s, ,(z) geniigt der DGl. 


a? d 
VrSuo@) =(# S+ eh +e —vls,,0) =e, (44) 


wahrend die Funktion-E,(z) von H. F. WEBER die Lésung der DGl. 


Pp 


Vp E,(2) = ee (1 —cos px) —— (1 + cos pa) (45) 


7 
darstellt. Fiir die weitere formelmaBige Rechnung ist es belanglos, welche der beiden Funk- 
tionen verwendet wird; die WEBERsche Funktion bietet jedoch den Vorteil, daB sie bei JAHNKE- 
EmDE [17] vertafelt ist. 

Behandelt man die DGl. (45) nach der Methode der Variation der Konstanten, so erhalt man 
sogleich den Zusammenhang mit den Zylinderfunktionen: 


E,(z) = = Jol) -{[ + cos px +2 (1 —cos pm) N,(2) dz 


—£ Nye) | [1 + cos p2 +2 (1 cos pa Jy le) de (46) 
und speziell fir # = 1 . 
E.(z) = Srl) | Me) de — NAC) | EM dz. (46a) 
In Verbindung mit der Beziehung 
250) — £2) —+ By +2 


lassen sich nun leicht die beiden Integralformeln ableiten 


[2 dz = = 2 [fp — hh Bol —h (47 a) 


und 


N,( 
[At ea:™2—MEI—M. (47) 


Zz 


Diese Formeln geniigen zur Lésung des in Gl. (41) noch ausstehenden Integrals; fiir weitere 
Einzelheiten tiber die Lommersche und WEBERsche Funktion sei auf [16] verwiesen. Analog 
zu den Gln. (23a, b) wird auch fiir die WEBERschen Funktionen die Abkiirzung eingefiihrt 


e(a =| = Eq{a7)— (1 —») al ra) (48) 
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Damit folgt nach geringer Zwischenrechnung 
he ( fee e( Ta) 4 _2tm 
n 
| U,,(7) « ; -dr= ’m Tl hy V5 \  %m HEL ey 
ze . 
v 


v, y 
n [+ | n (: i) 
aA 
i Y v 


™m™ 


(49) 


Jetzt ist nur noch die Anwendung des Induktionsgesetzes 11 = jmaw@® und des Ohmschen 
Gesetzes ft = U/J - &°? erforderlich, um— in Zusammenfassung der Gln. (41), (37), (43) und 
(49) — die Gleichung fiir den komplexen Scheinwiderstand §{ des magnetostriktiven Ring- 
Dehnungsschwingers anschreiben zu kénnen: 
Sy 1 @ Eg h tg» @Ege wh Oy of 
R=jo(1— . Hy: ws Soin @— jo" ae e (50) 


LF fe 


_Darin steht A, mit der Bedeutung: 


(a va) s= (l=) ae (x = + (1— ») eeulth =|? 


sable ms (50a) 


eal (7) 


Die Verwendung der Funktionen 7(z) statt wie bisher (z) erfolgt mit Gl. (39) und ist besonders 
fiir die Anwendung von Naherungen bequemer. 

Die Gl. (50) besteht aus zwei Anteilen, von denen der erste der Induktivitat des Schwingers 
entspricht, wie sie unabhangig von den Schwingungen allein durch die Wicklung gegeben ist. 
Der zweite Term beschreibt das eigentliche magnetoelastische Schwingverhalten, wobei die 
einzelnen Glieder des Reihenausdrucks in der Reihenfolge der Summationsvariablen n der 
Grundschwingung bzw. den weiteren Oberschwingungen zugeordnet werden. Ohne Magneto- 
striktion, d. h. mit d = 0, verschwindet dieser Term; dann ist der Ring magnetostriktiv auch 
nicht zu Schwingungen anregbar. 


6. Die Bestimmung der Eigenwerte a, und Vereinfachungen 
der Scheinwiderstandsgleichung 


Die Eigenwerte «a, sind die Nullstellen der charakteristischen Determinante Gl. (24). Im 
allgemeinen werden sie numerisch oder graphisch ermittelt werden miissen, da die Herstellung 
von brauchbaren N&herungen einen erheblichen Aufwand verlangt; so lassen sich auch die 
Ergebnisse zweier Arbeiten von BucHuHoLz [18], [19] hier nicht unmittelbar anwenden. 

An dieser Stelle interessiert nun hauptsachlich, wie die Eigenwerte von den geometrischen 
Abmessungen des Ringes abhangen. Dies laBt sich an Hand einiger Naherungsformeln zeigen, 
deren Anwendungsbereich allerdings z. T. begrenzt ist. 

Zunichst sei noch bemerkt, daB die Eigenwerte «, samtlich reell sind; dieses folgt unmittel- 
bar aus der Orthogonalitat der Eigenfunktionen. Wegen des als komplex angesetzten E- 
Moduls erhaélt man damit dann komplexe Eigenfrequenzen a,,. 

Bei den héheren Eigenwerten nun werden auch die Argumente der einzelnen Funktionen 
in Gl. (24) groB, so daB man vorteilhaft die asymptotischen Naherungen der Zylinderfunktionen 
benutzen kann. Fiir die Ausdriicke z(z) und v(z) nach Gl. (23) ergibt sich dann 


w= Vz sin (2 Petr sl cos(z + =} (51) 


n(z) =|al-~ i | a (e+). (51b) 


und 
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Bezeichnet man mit b = 7, —7; die Breite des Ringes und mit 7,, = = (ry, + 7,) seinen mitt- 
leren Radius, so ist mit B = 0/2 7,, 


14 = 1m (1 + B) und 1; == Yq (L — B) . 


Damit erhalt man nach Einsetzen der Gln. (51) in Gl. (24) und etwas Zwischenrechnung 


See Ghee ah a jae an _ 
1— f2 \2a 


Die Lésung dieser Gleichung erfolgt entweder graphisch oder numerisch. Die Eigenwerte 


ergeben sich im wesentlichen aus den Schnittpunkten der Kurve tg x mit der Hyperbel C a : 


je groBer die Werte fiir. 2 « £ werden, desto naher riicken diese Schnittpunkte der beiden Kurven 
an die Werte (n —1) a heran. Daraus folgt schlieBlich : 


ceive Coed eh (53) 


2B 


Mit wachsender Ordnungszahl n folgen die Eigenwerte im Abstand 7/2 6 aufeinander. 

Fiir die interessanteren Eigenwerte a, und a, lassen sich in einfacher Weise keine Naherungs- 
formeln angeben. Lediglich fiir den Sonderfall des diinnwandigen Ringes kann man sich einen 
guten Uberblick verschaffen. Bei diesen Ringen ist 6 <1, so da8 sich die Funktionen - 
t [a(1 + B)] und [a(1 + §)]in eine TAyLorreihe entwickeln lassen, von der noch Glieder mit 
(x 6)? beriicksichtigt werden. Man erhalt dann eine Gleichung 4. Ordnung fiir «,, deren Lésung 
nach kurzer Zwischenrechnung ergibt: 


a ae [1-2 9) O +9)-B), (5.4) 


owl | 4 (1 ») 6). (55) 


Wegen der schlechten Konvergenz der verwendeten TAyLorreihe ist dieses Ergebnis nur fiir 
B <1 brauchbar. Bei gréBeren f-Werten miiBten in der Reihe entsprechend mehr Glieder 
verwendet werden. Dann ist es jedoch wieder wesentlich einfacher, die Eigenwerte numerisch 
nach der Tabelle zu berechnen. 

Wie man nun aus diesen Naherungsformeln ersieht, ist bei dinnwandigen Ringen der kleinste 


1 
2 


Eigenwert a, yi—v? (d. h. die Grundfrequenz ist umgekehrt proportional dem mittleren 
Durchmesser), wahrend fiir die folgenden Eigenwerte die Breite des Ringes maBgebend ist. 
Die Angabe ahnlicher Naherungen fiir den Faktor A, nach Gl. (50a), der fiir jeden Eigen- 
wert a, einen speziellen reellen Wert annimmt, ist bereits erheblich schwieriger. Deshalb soll 
sein Verhalten hier nur in groBen Ziigen erlautert werden. 
Nimmt man bei den héheren Eigenwerten auch fiir die WEBERschen Funktionen die asymp- 
totischen Naherungen, so ist 


e(z) w n(z) ——, (56) 


UZ 
wenn Glieder mit héherer Ordnung als 1/z vernachlassigt werden. Damit erhalt man fiir AS 


4 (a, (1 — B)] rw Pre eer B)) r 
AS 9 COMA ies aye)) Xy (1— B) le 
eee i? [ay (1— B)] —2? [a, (1 + B)) © (57) 


Bei Verwendung von Gl. (51a) werden in den trigoriometrischen Funktionen die Argumente 
a, ° (1 +) im allgemeinen keine ganzzahligen Vielfache von z/2 sein. Daher ergeben sich 
auch fiir die A, — zumal auch wegen der Differenzen in Zahler und Nenner — keine einsinnig 
mit wachsendem m abnehmenden Werte, sondern eine unregelmaBige, je nach dem Verhaltnis 
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von 7,/r, sich andernde Verteilung, bei der jedoch A, = 0 und A, =o nicht auftreten. Die 
einzelnen Oberwellen, die auch nur fiir y + o existieren kénnen, lassen sich demnach nicht 
gleich stark anregen. 

Im Falle diinnwandiger Ringe (f <1) liefert die schon im vorigen angewendete TAYLOR- 
- Entwicklung nur fiir die Grundschwingung ein brauchbares Ergebnis. Unter Beriicksichtigung, 


y, (2) 
6 
5 
e(z)+(1-/Z 
4 Y, (Z)= UZ) 
3 


ae | 


=| 


\ 
V=03 
-2 v=0,25 


-4 + 


—5/- 


Bild 2. Verlauf der Funktion y,(z) nach Gl. (59a). 


daB mit Gl. (54) 0, © /1 — ist, ergibt sich fiir A, nach einiger Zwischenrechnung die einfache 
Beziehung 
Ay ¥ (1— 04) 28. (58) 


SchlieBlich sei noch darauf hingewiesen, daB man sich eine gewisse Orientierung iiber die Ver- 
haltnisse bei breiteren Ringen an Hand der Diagramme in Bild 2 und 3 verschaffen kann. 
Die Gl. (50a) ist in einer Form geschrieben, bei der sich zwei Funktionen y,(z) und p,(z) erkennen 
lassen, 


e (2) + (Q—») — 
(2) = - ; (59a) 


und 


NS Sear eee (59 b) 
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die in den Diagrammen fiir zwei verschiedene Werte von » dargestellt sind. Daraus kénnen bei 
bekannten «, und f die zugehérigen Werte von y, und y, entnommen werden. Es zeigt sich 
dabei, daB bei gegebenem 7, und 7, die GréBe von A, u. U. auch von dem Wert von v abhangt, ~ 
und zwar besonders dann, wenn die Argumente «, - (1 +) in der Nahe der Pole der beiden 
Funktionen liegen. - 


Y,(2) 
60 

7 1-(1-4)4 
a HOI Tal 


40 


30 


20 


Bild 3. Verlauf der Funktion y,(z) nach Gl. (59b). 


7. Ersatzschaltbild und Ortskurve 


Wie sich im vorigen gezeigt hatte, lassen sich lediglich fiir die Grundschwingung bei 
diinnwandigen Ringen einfache Naherungsformeln sowohl fiir den Eigenwert a, als auch fiir 
A, angeben. Fiir diesen Fall nimmt dann die Scheinwiderstandsgleichung (bei Beriicksichtigung 
nur des ersten Reihengliedes in Gl. (50)) die schon bekannte Form an. 

Mit Gl. (58) und In (7,/7,) & b/r,, erhalt man zunachst 


2 it (Podeye, w2-b-h : w-b-h dEy a? 
R=jo(s ; ; ) ens — —] WU,- 5 4 1 
Sa 


= 5 no 
ZV y, ZUVy, Ho A — ay 


(60) 


Diese Gleichung kann unter Verwendung des magnetoelastischen Kopplungsfaktors k weiter 
vereinfacht werden. Die GréBe k? gibt das Verhaltnis an, in welchem die dem Schwinger 
zugefiihrte magnetische Energie in mechanische Energie umgewandelt werden kann, und zwar 
wieder bei Frequenzen weit unterhalb der Resonanzfrequenz, so daB die kinetische Energie 
unberiicksichtigt bleibt. 


ae lo|? 
(yeni! 
oe eras (61) 
= ABR 
eile 
Ferner gelten die Beziehungen 
1 1 
E ia A: @E 
Re H B Ho Me H (61 a) 
nae Mo. Lo 
Ey 
und 
7.) . . 
k = Re a) Cueene) eae ae 671 Che + eB) (61b) 
1— k? 1 ke Me: Ey 
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Dabei sind e = (= =(3") eine weitere Magnetostriktionskonstante, 6,. der Ver- 
0H] ¢=const de /H=const ¥ 
lustwinkel der Permeabilitaét uw, und 6; y der des Elastizitatsmoduls E,. AuBerdem ist 
eer, eb h 46 
ji 10U5, oe w 10, 
ene || srs é (62) 


die Induktivitat des eingespannten (daher y,), nicht schwingenden Ringkernes. Bei den diinn- 
wandigen Ringen kann man in der Grundschwingung den Einflu8 der Querkontraktion ver- 
nachlassigen und 1 — y? & 1 setzen. 

In dem Quotienten a} /(«® — a2) ist mit Gl. (14) die laufende Frequenz w der Anregung 
und die komplexe Frequenz w, der Grundschwingung enthalten. Mit Gl. (54) ist also 


1 |Ex| BOR, , A Ory 
Oy ale ( ce erat arene ean ea (03) 
und demnach 
at 10 
2 P 2 16 ee 
ie a ee oe Ey (64) 
w—at wr—w? 2S haiat (< =) 
+4 . as 
@, Of, \®, ©® 


In der nahen Umgebung der Resonanzfrequenz w;, die wegen des bei Ferriten sehr kleinen 
Verlustwinkels Ory (~ 10) extrem scharf ausgepragt ist, kann man ohne wesentlichen Fehler 
w/o, ~~ 1 setzen. SchlieBlich ist 1/E,, = Q, die mechanische Giite des Schwingers, so daB fiir 
den Scheinwiderstand die Beziehung resultiert 


—j26 
| One Me 


1— k? : o w,\" 
1 +108(5—2) 
(60) 


1 


ie : 
K=jo Lye Pete A (65) 


Ublicherweise wird nun mit L, |k?/(1 —.k2)| eine neue Induktivitat L,, definiert, bei der der 
Index m darauf hinweist, daB sie durch rein magnetostriktiv-mechanische Gré8en bestimmt 
wird; mit Gl. (61b) und (62) ist namlich (¢ = Querschnitt 0 - h) 


2 | 92 
Bees 


Ly = ; (66 a) 


2UVy, Eq 
setzt man nun w=1/L,,C, und R,, = w, L,, 0; = Q,/w, C,,,. so erhalt man die von 
SIxTus [3] angegebenen Beziehungen 


C. = 22a), i wens (66b, c) 
m 472 w? q lel? , m 0-2: Ory 4 72) 


Mit der Einfiihrung der eben genannten GroBen L,,, C,, und R,, wird bereits erkennbar, daB 
es sich bei dem zweiten Term der Gl. (65) um einen Parallelresonanzkreis mit der Resonanz- 
frequenz , = w; handelt, der den magnetostriktiv angeregten Schwingungen entspricht. 
Der komplexe Scheinwiderstand %{ eines diinnwandigen Ring-Dehnungsschwingers in der 
Grundfrequenz laBt sich demnach durch das in Bild 4a gezeichnete Ersatzschaltbild be- 
schreiben. Die zugehérige Ortskurve ist in Bild 4b gezeigt.1 Sie besteht im wesentlichen aus 
einem um 7 w, Ll, verschobenen Kreis, der dem Parallelresonanzkreis des Ersatzschaltbildes 
entspricht. Der durch den Resonanzpunkt w, gelegte Durchmesser des Kreises ist um den 
Winkel 2 6, gegen die Horizontale geneigt. Um diese Neigung auch im Ersatzschaltbild wieder- 
zugeben, miiBte in diesem exakterweise zwischen L, und dem Parallelkreis ein Ubertrager mit 
dem komplexen Ubersetzungsverhaltnis e:1 gezeichnet werden (s. [4] u. [6]). In der Praxis 
wird er jedoch meist fortgelassen, womit man die duale Form zu dem bekannten Ersatzschalt- 
bild des Quarzschwingers erhalt. 


1 Diese Darstellung gilt im Prinzip fiir alle magnetostriktiven Schwinger. 
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Wie aus der Ortskurve ersichtlich, tritt auBer der Parallelresonanzfrequenz @, noch eine 
zweite, die Reihenresonanzfrequenz w, auf. Aus dem Abstand der beiden Frequenzen, der sog. 
Bandbreite Aw, laBt sich in guter Naherung der Kopplungsfaktor bestimmen: 

US Ey my Ae e Pas Dy (67) 

Ihe Op Wp 
Bisher ist in Ersatzschaltbild und Ortskurve nur die Grundfrequenz beriicksichtigt. Die Ober- 
schwingungen ergeben, entsprechend der Reihe in Gl. (50), eine Reihenschaltung von (theore- 
tisch unendlich vielen) Parallelkreisen mit wechselnd groBem R,,, je nach der GréBe des je- 


Im) 


Bild 4. Ersatzschaltbild (a), Ortskurve (b) und Verlauf des Scheinwiderstandsbetrages (c) magnetostriktiver Schwinger. 


weiligen A,. In der Ortskurve werden mit steigender Frequenz eine Reihe von Resonanzkreisen 
durchlaufen, die, entsprechend der GréBe von R,,, verschieden groBe Durchmesser haben, 
und deren Mittelpunkte immer weiter in Richtung der imaginaren Achse wandern. Bei den 
héheren Oberschwingungen ist diese Verschiebung gréBer als der Durchmesser, so da kein 
Schnittpunkt mehr mit der reellen Achse zustande kommt. In Bild 4c schlieBlich ist der Ver- 
lauf von |i(@)| skizziert. 

Mit den soeben angegebenen Beziehungen ist der Anschlu8 an bereits bekannte Formeln 
hergestellt worden. Andererseits bieten sie die Méglichkeit zur Orientierung iiber den Inhalt 
der Gl. (50) besonders auch fiir die Falle, in denen die Herstellung einfacher Naherungen nicht 
méglich ist. 


8. Experimentelle Ergebnisse 


Bei der Erorterung der elastizitatstheoretischen Grundlagen im 3. Abschn. wurde fiir die 
dann folgende Rechnung der ebene Spannungszustand zugrunde gelegt und dabei voraus- 
gesetzt, daB sich die mechanischen Spannungen und Dehnungen in Richtung der Achse des 
Schwingers nicht andern. In Wirklichkeit treten aber wegen der Querkontraktion auch perio- 
dische Dehnungen in dieser Richtung auf. Ihre Beriicksichtigung in der DGl. fiir die elastischen 
Schwingungen hatte auf ein nur naherungsweise losbares Randwertproblem gefiihrt. Diese 
Komplikation trat bei der Annahme des ebenen Spannungszustandes zwar nicht auf, dafiir ist 
nun aber damit zu rechnen, daB die so gewonnenen Ergebnisse nur in einem gewissen Frequenz- 
bereich giiltig sind. Der entstehende Fehler ist an dem Unterschied zwischen den meBbaren 
und den sich aus der Rechnung ergebenden Eigenfrequenzen erkennbar. Um einen Uberblick 
zu gewinnen, wurden deshalb bei sechs verschiedenen Ringschwingern die Eigenfrequenzen 
gemessen. 

Mit einem PreBwerkzeug zur Herstellung von kreisférmigen Platten mit ca. 45% mm 
AuBendurchmesser wurden zunachst einige Scheiben aus Ferritschwinger-Material gepreBt und 
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anschlieBend auf einer Drehbank mit verschiedenem Innendurchmesser ausgebohrt. Die Ringe 
wurden dann gesintert und in passenden Schalen dampfungsfrei gelagert. Die erforderliche 
Vormagnetisierung erhielten sie mittels eines kurzzeitigen Gleichstromes durch eine Wicklung 
von 30 Windungen, wobei das Material bis anndhernd zur Sattigung aufmagnetisiert wurde. 
Der Arbeitspunkt ist dann der Remanenzpunkt der Magnetisierungskurve. Uber die gleiche 
Wicklung sind anschlieBend die Ringe mit einem Wechselstrom von hinreichend kleiner Ampli- 
tude zu Schwingungen angeregt worden. 

In Tabelle 1 sind in Spalte 1 die geometrischen Abmessungen der Ringe angegeben. Die 
Spalte 3 enthalt die zugehérigen Eigenwerte «,. Fiir die Grund- und die erste Oberschwin- 
gung (” = 2) sind die Werte aus Gl. (24) ausgerechnet worden. Die weiteren Eigenwerte 
wurden nach Gl. (52) bestimmt, also unter Verwendung der asymptotischen Naherungen. 
Um einen guten Ubergang zu haben, sind die Werte fiir die erste Oberschwingung wiederholt. 
Es zeigt sich, daB die Naherungen schon bei relativ kleinen Argumenten gut brauchbar sind. 
Alle Eigenwerte sind mit einem im allgemeinen auch fiir Ferrite zutreffenden Wert von vy = 0,3 
bestimmt worden. 

Die nachste Spalte enthalt dann die gemessenen Eigenfrequenzen, wie man sie in einfacher 
Weise aus der Messung von |t(w)| gewinnt. Bei exakter Ubereinstimmung zwischen Rechnung 
und Messung miiBte sich nun aus der Gl. (14) mit 


y= [El = (1 —_ py?) (2: mn 


Q Xn, 


immer die gleiche Schallgeschwindigkeit v = v, der longitudinalen Schallwellen ergeben. Das 
ist jedoch nicht der Fall; mit steigender Ordnungszahl der Eigenfrequenzen wird v kleiner, und 
zwar unterschiedlich bei den einzelnen Ringen. Um die v-Werte bequemer miteinander ver- 
gleichen zu konnen, sind sie jeweils auf den Wert fiir v, in der Grundschwingung bezogen 
worden (Spalte 5); aus allen sechs Ringen ergibt sich ein Mittelwert v,; = 5,8975 - 10° [cm/s] 
+0,2%. Zur besseren Ubersicht sind die jeweils zu einem Ring und seinen diskreten Eigen- 
frequenzen gehdrenden Werte v/v, durch einen Kurvenzug verbunden in Bild 5 dargestellt. 


V Ring Nr —3 
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Bild 5. Die Abhangigkeit der Schallgeschwindigkeit von der Frequenz bei Ring-Dehnungsschwingern mit verschiedenem Verhaltnis h/b. 


Die Anderung der Schallgeschwindigkeit mit der Frequenz ist bei Dehnungsschwingern 
eine bekannte Erscheinung. Sie wird dadurch hervorgerufen, daB die in den achsialen Deh- 
nungsschwingungen enthaltene Energie bei Annaherung an die Resonanzfrequenz dieser Schwin- 
gungen nicht mehr vernachlassigbar klein ist. Fiir die im Prinzip gleichen Verhaltnisse bei den 
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Langsschwingungen zylindrischer Stabe, bei denen der Durchmesser nicht mehr klein ist gegen- 
iiber der Linge, ist schon von RAYLEIGH die bekannte Korrekturformel fiir die Schallgeschwin- _ 
digkeit angegeben worden. 

Wie im Abschn. 6 gezeigt wurde, ist bei den radialen Oberschwingungen die Breite 6 die 
frequenzbestimmende Dimension des Ringes, wenn die Querabmessung, also die Hohe h, 
hinreichend klein ist. Je kleiner h/} ist, desto weniger andert sich auch die Schallgeschwindig- 
keit, weil dann bei einer gegebenen Frequenz unterhalb der Grundfrequenz der achsialen 
Dehnungsschwingungen der in diesen Querschwingungen enthaltene Anteil der kinetischen 
Energie entsprechend klein ist. Der Grenzfall h/b + 0 liefert den Fall des ebenen Spannungs- 
zustandes mit der konstanten Schallgeschwindigkeit v = v,, wie er in der Rechnung ange- 
nommen wurde. 

Die skizzierten Kurven nahern sich mit wachsender Frequenz asymptotisch dem Wert 
v7,/0, 0,6 (fiir y = 0,3), d. h. die Schallgeschwindigkeit geht von dem Wert fiir die Longi- 
tudinalwellen in den fiir die Transversal-(Scherungs-) Wellen tiber. 

In der achsialen Grundschwingung ist die Héhe des Ringes gleich der halben Wellenlange, 
h = j/2. Mit v = /- f ergeben sich aus den Eigenfrequenzen dieser Schwingungen, die bei allen 
Ringen tiber 500 kHz liegen, fiir die Schallgeschwindigkeit héhere Werte als der oben angegebe- 
ne Mittelwert fiir die radiale Grundschwingung; es ist dabei aber zu beriicksichtigen, daB 
diese Werte schon auf héheren Asten der Dispersionskurven liegen. 

Die Identifizierung dieser Schwingungen ist bereits recht schwierig, da bei den unter- 
suchten relativ breiten Ringen mit steigender Frequenz — ungefahr oberhalb 300 kHz — 
eine wachsende Vielzahl von Resonanzfrequenzen anderer Schwingungsformen auftritt, zu 
denen auch Torsions- und Scherungsschwingungen gehéren kénnen, sowie Kopplungen der 
verschiedenen Schwingungstypen. Aus diesem Grunde sind die Resonanzfrequenzen in Zahlen- 
tafel 1 nur bis 7m = 4 angegeben worden. 

In diesen unerwiinschten Nebenschwingungen bei héheren Frequenzen, deren Unter- 
driickung ein wesentliches Problem in der Praxis aller mechanischen Schwinger darstellt, 
liegt — auBer der genannten Anderung der Schallgeschwindigkeit — ein weiterer Grund fiir die 
begrenzte Anwendbarkeit der im theoretischen Teil dieser Arbeit gewonnenen Ergebnisse. 
Wie Bild 5 zeigt, gelten diese demnach streng fiir die Grundschwingung und auch noch fiir die 
erste Oberschwingung, sofern das Verhiltnis //b nicht zu groB ist. 

Fiir die Grundschwingung der sechs Ringe sind einige Me8werte in Zahlentafel 2 enthalten. 
Die Eigenfrequenzen /, zeigen eine gute Ubereinstimmung mit der Rechnung. Mit kleiner 
werdendem f nimmt auch der Wert fiir A, ab und ebenso auch R,,. Das bedeutet jedoch keine 
Verminderung der mechanischen Giite Q,, da L, sich ebenfalls verringert. Giite und Kopplungs- 
faktor k sind als reine WerkstoffgréBen unabhangig von den Dimensionen des Ringes. Sie 
werden nur durch die chemische Zusammensetzung, die Sinterbehandlung und den Magneti- 
sierungszustand bestimmt. Die in der Tabelle auftretenden Schwankungen dieser GréBen 
liegen innerhalb der erfahrungsgemaBen Grenzen. 


9. Zusammenfassung und SchluB 


Ziel der Arbeit ist die Bestimmung des komplexen Scheinwiderstandes von magnetostrikti- 
ven Ring-Dehnungsschwingern. Ausgehend von den magnetostriktiven Vierpolgleichungen 
und den Spannungs-Dehnungsgleichungen der Elastizitatstheorie wird unter Annahme des 
ebenen Spannungszustandes die DGl. fiir die erzwungenen elastischen Schwingungen des Ringes 
aufgestellt, die zusammen mit zwei Randbedingungen eine inhomogene Randwertaufgabe 
bildet. Diese wird mit Hilfe der Greenschen Funktion gelést, bei gleichzeitiger Entwicklung 
der Lésung in eine Reihe nach den Eigenfunktionen desSystems. Dadurch stellt sich auch der 
Scheinwiderstand i als ein Reihenausdruck dar, dessen einzelne Glieder der Grund- bzw. den 
Oberschwingungen entsprechen. Fiir den Sonderfall des diinnwandigen Ringes in seiner Grund- 
frequenz kann das Ergebnis auf bereits bekannte Formeln zuriickgefiihrt werden. 
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Zahlentafeli 


4 5 
Re ay {cm] | gemessene 
ae he Res. Frequ. hs 
hn 
[kHz] 
| 76,067 
i 4,390 2,1015 2AZ 322 0,9884 
0,640 2 2,1649 213,322 1,0096 
1 (ies 3 3.9646 =| = 356,288 0,9124 
0,593 ee 5,8455 | 506,756 0,8802 
0,391 5 7,7467 
1,350 1 0,7285 71,588 1,0000 
45395 2, 2,2733 220,630 0,9876 
: 0,645 2 2,2513 220,630 0,9973 
2 3 4,2208 369,950 0,8919 
0,530 4 6,2529 524,830 0,8541 
0,424 5 8,3007 
1,545 1 0,6940 68,448 1,0000 
4,385 2 2,3753 230,295 0,9829 
0,045 2 2,3058 230,295 0,9869 
3 3 4,5015 384,498 0,8660 
0,479 4 6,6883 546,455 0,8283 
0,454 See eo, OS0O 
1,775 1 0,6598 64,804 1,0000 
4,390 2 2,5369 245,188 0,9840 
0,640 2 2,5250 245,188 0,9884 
4 3 4,8669 402,076 0,8411 
0,424 4 7,2492 576,010 0,8090 
0,489 5 9,6418 | 
2,045 1 0,6214 61,189 1,0000 
4,385 2 2,7848 266,882 0,9732 
0,640 2 2,7807 266,882 0,9746 
5 3 5.4181 429,150 0,8043 
0,364 4 8,0872 601,445 0,7552 
0,547 5 10,7644 
2,485 1 0,5704 56,056 1,0000 
4,385 2 3737281), 3805530 0,9555 
0,635 2 BByOS || BUlOy 530 0,9555 
6 3 6,6459 | 474,903 90,7271 
0,276 4 9,9422 650,100 0,6654 
0,668 5 13,2440 


Zahlentafel 2 
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Da in der Rechnung die achsialen Dehnungsschwingungen, die infolge der Querkontraktion 
auftreten, nicht beriicksichtigt sind, ergeben sich Abweichungen bei der in Wirklichkeit fre- 
quenzabhangigen Schallgeschwindigkeit v, womit auch die errechneten Eigenwerte nicht mehr 
genau stimmen. Zur Abschatzung dieses Fehlers und des Anwendungsbereiches der Schein- 
widerstandsgleichung werden die MeBergebnisse von sechs Ringschwingern verschiedener Breite. 
diskutiert. 

Fiir die wertvollen Hinweise und Ratschlage, mit denen Herr Prof. Dr. H. BucHHoZLz, 
TH Darmstadt, das Zustandekommen dieser Arbeit geférdert hat, méchte der Verfasser auch 
an dieser Stelle seinen besten Dank sagen. Auch Herrn Dr. S. SCHWEIZERHOF im Hause Tele- 
funken GmbH. sei fiir klarende Diskussionen auf dem Gebiete der magnetischen Werkstoffe 
vielmals gedankt. 
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Allgemeine Grundlagen der Netzberechnung mit Inzidenzmatrizen 
Von 
HANS EDELMANN, Erlangen 
Mit 9 Textabbildungen 
(Eingegangen am 22. Juni 1959) 


Ubersicht. Zur Berechnung der Stréme und Spannungen in elektrischen Netzen bei gegebenen treibenden 
Spannungen verwendet man die Maschen- oder Knotenpunktsmethode. Wahrend die erste Methode die Stréme 
unmittelbar liefert, ergibt die zweite Methode die Spannungen unmittelbar. Zur iibersichtlichen Aufstellung 
der Gleichungen ist es zweckmaBig Indizenzmatrizen, welche die Topologie des Netzes eindeutig beschreiben, 
zu verwenden. AuBerdem kann die Theorie der Inzidenzmatrizen mit Vorteil dazu benutzt werden, das Auf- 
stellen der Gleichungen fiir die unbekannten Stréme bzw. Spannungen einem Rechenautomaten zu tiberlassen, 
wodurch es méglich wird, Netzberechnungen auch weniger getibten Hilfskraften zu tibertragen. Bei Ver- 
wendung eines digitalen Rechenautomaten sind dann die Netzdaten nur noch in Form einer Liste aller Impe- 
danzen und der mit ihnen verbundenen Knotenpunkte einzugeben. 

Um den reibungslosen Ablauf einer Netzberechnung zu garantieren, ist es ferner zweckmaBig, gewisse 
hinreichende Kriterien fiir die Auflésbarkeit der Netzwerksgleichungen zu beachten, die hier bewiesen werden. 


Liste der benutzten Symbole 


k = Anzahl der Knotenpunkte, Y/ =Z1=G+ 7D = Zweigadmittanz- 
Z = Anzahl der Zweige, matrix, 
2p = Anzahl der Baumzweige, 1 , ; 
zy = Anzahl der unabhangigen Zweige, c se f= SH ee es 
y = Anzahl der unabhangigen Maschen bzw. matrix. 
Zweige, i? = Vektor der Maschenstroéme, 


b = Gesamtzahl aller méglichen Baume, 

S = komplexe Leistung = P+7Q, 

P = Wirkleistung, 

Q = Blindleistung. 

0 . = Nullmatrix, -vektor, 

1 = E/inheitsmatrix, 

M_ = Quadratische oder rechteckige Matrix, 

M, = transponierte von M, 

M* = Matrix der konjugiert-komplexen Werte 


von M, 
My = hermitescher Teil von M, 
vy = Spaltenvektor, 
v;, = i-ter Zeilenvektor der Matrix V, 


Vv; = k-ter Spaltenvektor der Matrix V, 

i = Vektor der Zweigstréme, 

u,e = Vektor der Zweigspannungen, bzw. der 
Zweig-EMKe., 

Ze = diae(zq-. . 2.) = R +7 X = Zweig- 
Impedanz-Matrix, ; 


i = Vektor der iibrigen Zweigstréme, es ist 


j9 
also j = é ) 
i 


uv? = Vektor der Maschenspannungen, 

Vhs = Maschenimpedanzmatrix. 

K = Knotenpunktszweig-Inzidenzmatrix, 
K = (A, B) = reduzierte Knotenpunkts- 


Zweiginzidenzmatrix, Zeile des Bezugs- 
knotenpunkts gestrichen, 


, it = Vektor bzw. reduzierter Vektor der Kno- 
tenpunktsstréme, 
ut, ut = Vektor bzw. reduzierter Vektor der Kno- 
tenpunktsspannungen, 
Yt, Y+= Knotenpunkts-Admittanzmatrix bzw. 
reduzierte Knotenpunkts-Admittanz- 
matrix. 


1. Allgemeine Bemerkungen 


Im Zuge einer immer weiteren Vervollkommnung automatisierter Netzuntersuchungen der 
Energie-Ubertragungstechnik wird man bestrebt sein, diese ftir den damit Beauftragten so 
einfach wie méglich zu gestalten, so daB die Méglichkeit von Irrtiimern auf ein Minimum redu- 
ziert wird. Hierbei wird die Eingabe der Netzdaten weitgehend mit dem entsprechendenVor- 
gang bei einem Netzmodell iibereinstimmen. Auch dort beginnt man damit, alle Knoten- 
punkte und alle Zweige durchzunumerieren, alsdann kann man eine Liste anfertigen, welche 
die Netzzweige den beiden begrenzenden Knotenpunkten zuordnet. Mit der Unterscheidung 
zwischen einer 1. und 2. Knotenpunktsnummer des Zweiges, den sie begrenzen, ist auch eine 
Orientierung festgelegt. Sie kann bei der aufzustellenden Inzidenzmatrix verwendet werden. 
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Diese Liste enthilt eine vollstandige Information tiber die Topologie des betreffenden Netzes. 
Eine weitere Liste enthalt QuantitatsgréBen, und zwar die Impedanzen oder die Admittanzen 
der einzelnen Netzzweige nach den Zweignummern geordnet. SchlieBlich miissen noch die in 
gewissen Zweigen wirksamen EMKe angegeben werden. Dies kann in einer weiteren Liste nach 
Zweignummern geordnet geschehen. 

Fiir die topologische Charakterisierung des Netzes verwendet man zweckmaBig Inzidenz- 
matrizen. Derartige Matrizen sind algebraische Darstellungen der Topologie des Netzes und 
demnach invariant gegeniiber quantitativen Anderungen der Zweigimpedanzen (-admittanzen). 
Da die letzteren GréBen im Verlauf einer Untersuchung gelegentlich geandert werden miissen, 
ist es zweckmaBig, topologische und quantitative Informationen aus dem gegebenen Netz zu 
trennen. 

Ohne auf spezielle von der zur Verfiigung stehenden Rechenanlage abhangige Programme 
naher einzugehen, sollen hier nur allgemeine Grundlagen fiir die Aufstellung der Gleichungs- 
systeme zur Ermittlung der Spannungen bzw. Stréme behandelt werden, im wesentlichen also 
die beiden grundlegenden Methoden: Maschen- und Knotenpunktsmethode, Ermittlung eines 
vollstandigen Baumes aus der Zweig-Knotenpunkts-Matrix bzw. aus der Liste aller Zweige, 
der Zusammenhang der beiden Methoden und die Berechnung der Maschen-Zweig-Inzidenz- 
matrix aus einer Normalform der Knotenpunkts-Zweig-Inzidenzmatrix. SchlieBlich werden 
noch einige hinreichende Satze iiber die Existenz von eindeutigen Lésungen bewiesen. 

Im iibrigen soll hier nur der praktisch wichtigere Fall, daB das Netz ausschlieBlich von 
treibenden Spannungen erregt wird, behandelt werden. Eine Verallgemeinerung auf treibende 
Stréme macht dann keine Schwierigkeiten. AuBerdem wollen wir uns zunachst auf Netze ohne 
Transformatoren (RCL-Netze) beschranken. Die Beriicksichtigung von Transformatoren 
und die Ermittlung gewisser abgeleiteter Matrizen (bezogen auf ausgewahlte Klemmenpaare) 
wird in einer weiteren Arbeit gezeigt. 

Wir betrachten zundchst die beiden grundlegenden Methoden: Die Maschenmethode und 
die Knotenpunktsmethode. 


2. Die Maschenmethode 


In dem zu untersuchenden Netz (Streckenkomplex, Graph) mit k Knotenpunkten und 
z Zweigen (Strecken) ermittelt man einen vollstandigen Baum dadurch, daB man von einem 
beliebigen Knotenpunkt ausgehend immer nur solche Zweige hinzunimmt, die eine Verbindung 
zu einem neuen Knotenpunkt herstellen, bis alle Punkte miteinander einfach verbunden 
sind. (Schrittweise VergréBerung des Baumes durch Hinzunahme neuer ,,Aste‘‘!) Dies ist 
genau k —1mal méglich. Die Anzahl der Baumzweige (engl.: tree branches) ist also 


Zp =kR—1. 


Die tibrigen Zweige werden unabhangige Zweige (engl.: cotree links, air branches) genannt, 
da sie mit den Baumzweigen zusammen unabhangige Maschen (Cyklen) bilden. Der Baum 
selbst enthalt keine Maschen, denn durch die Hinzunahme nur solcher Zweige, welche zu 
einem neuen Knotenpunkt fiihren, wird die Bildung von Maschen verhindert. Ein zu einem 
unabhangigen Zweig gehériger Cyklus kann iiber den Baum nur auf eine Weise geschlossen 
werden, andere tells miiBte der Baum Maschen enthalten. Die Anzahl der unabhangigen 
Zweige zy ist also gleich der Anzahl der unabhangigen Maschen 7 


By Ts 
Fir die Gesamtzahl aller Zweige gilt also 
Z=%+ep=r+t+k—i1. 
Hat man einen vollstandigen Baum ermittelt (dies ist im allgemeinen auf mehrere Weise 
méglich), so numeriert man die Zweige mit den unabhangigen Zweigen beginnend durch und 


gibt jedem Zweig einen Orientierungspfeil (Definitionspfeil). Beginnend mit einem unab- 
hangigen Zweig zeichnet man die unabhangigen Cyklen ein; die Cyklenorientierung wird 


Crs eer H. EpELMANN: Allgemeine Grundlagen der Netzberechnung mit Inzidenzmatrizen 4921 


zweckmaBigerweise! in Ubereinstimmung mit der Orientierung des zugehérigen unabhangigen 
Zweiges gewahlt. 


Alsdann kann man folgende Matrix konstruieren: In ein rechteckiges Schema mit z Zeilen 
und 7 Spalten tragt man in der ¢-ten Zeile und o-ten Spalte Zahlen + 1, — 1 und 0 ein, je nach- 
dem, ob in dem o-ten Cyklus der ¢-te Zweig die gleiche, entgegengesetzte Orientierung hat 
oder nicht enthalten ist. Wir betrachten hierzu ein Beispiel. 

Beispiel 1: 


In Bild 1 sehen wir eine WHEATESTONE-Briicke, die diinn eingezeichneten Zweige (ig) 
sind die unabhangigen Zweige, die dick eingezeichneten (4, 5, 6) die Baumzweige. Es ist 


ara aye 


LSA, f= 3, 
r=2Z—k +1 =3 = 2. 


Das dem orientierten Graph eindeutig zugeordnete Schema lautet 


o= 4 2 3 
f= 
1 +1 (6) O 
B ) +1 ) 
z - He = 16 (1) 
Geer Patan S 
5 —1 —oO —1 
6 +41 1 +1 


Diese Matrix wird die Zweig-Maschen- oder auch Zweig-Cyklen-Inzidenzmatrix 
(Verkniipfungsmatrix) genannt. Mit Hilfe dieser Matrix laBt sich das vollstandige Gleichungs- 
system fiir die Maschenstréme aufstellen. Im folgenden sind 72,...7, (z = 6) die Zweigstréme, 
u,...u, die Zweigspannungen, 79...%7 (7 =3) die Maschenstréme, u?...u? die Maschenspannungen, 
ferner %,...2,, die Zweigimpedanzen und zugleich die Diagonalelemente der Zweigimpedanz- 
matrix Z, ¢,...¢, treibende Spannungen in den Zweigen. Es gelten dann die folgenden Be- 
ziehungen 


_=+7 
1p = + ty 
oe ae (2) 
1 +13+ 7 
Cs ij 13 — 45 
Wee 7, +49 
in Matrizen 
be Cute 
ferner 
ut = Uy, = hs aE Ug ; 
U= Uy +%y—4U;, 
Un = Uz + Us — Us + Ug 
in Matrizen 
Ue eC (3) 


1 Hierdurch wird im oberen Teil von C (durch Strichelung abgeteilt) eine Einheitsmatrix erzeugt. 
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My = 41h ah 

Uy = 2a Ue + € 

Us = 233 1 + é3 

UT Ges 24414 + eq 

Tea 25515 “es 

Ug = 266 %6 + & 


in Matrizen 
U = Zs ae" (4) 
Die Diagonalmatrix Z wird Zweigimpedanzmatrix genannt. 
Die KircHHOoFrFschen Knotenpunktsgleichungen sind in Gl. (1) enthalten, die Maschen- 
gleichungen in Gl. (3), wenn man 
u° = 0 (5) 
setzt, wahrend Gl. (4) die Anwendung des Ohmschen Gesetzes darstellt. Ausgehend von Gl. (3) 
gewinnt man durch Einsetzen von (2) und (4) eine Matrizengleichung 


Uses Ce(Z dee) 
5CZ,C Meee 3 

=/7° Gre 0% (6) 

wobei die Matrix 
7 Me (7) 
241 + %55 + 266 255 255 + 266 
oS 255 Z9 + 244 + 255 24a + 255 (8) 
255 + 266 aa + 255 233 + 244 + 255 + 266 


Maschen-Impedanzmatrix genannt wird und in Netzen, welche ausschlieBlich Impe- 
danzen (Zweipole) enthalten, immer symmetrisch (zur Hauptdiagonale) ist (Beweis!). Diese 


Ersatzschaltung 
eines Zweiges 


Unabhdngiger Zweig 


Baumzweig 


Unabhangige Masche 


Bild 1. Das topologische Bild einer WuEatstonEschen Briicke. 


Matrix kann man auch ohne Bildung des Produkts C, Z C direkt aus Bild 1 ablesen. Fiir die ein- 
zelnen Elemente gelten, wie man leicht nachweist, die folgenden Regeln: 

| 1. In der Hauptdiagonale ist das Element 2?; gleich der Summe aller Impedanzen, 
| welche der 7-ten Masche angehéren (Eigenimpedanzen). 

2. AuBerhalb der Hauptdiagonale ist das Element z°, = 2°, gleich der algebraischen 
| Summe aller Baumimpedanzen, welche gleichzeitig der i-ten und k-ten Masche angehoren. 
| Das Vorzeichen einer solchen Baumimpedanz in dieser Summe ist positiv oder negativ, 
| je nachdem, ob die beiden Maschen gleichsinnig oder gegensinnig iiber die betreffende 

Impedanz laufen (Koppelungsimpedanzen). 
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Die Orientierung der Baumzweige ist hier vollstandig verloren gegangen, nicht aber die 
Orientierung der Cyklen und zwar nur die gegenseitige iiber die Baumzweige. Man erkennt 
dies daran, da8 bei Umkehrung der Orientierung der Baumzweige die Matrix (8) sich nicht 
andert. Wenn man jedoch die Orientierung eines unabhangigen Zweigs und nach unserer 
Festsetzung gleichzeitig auch die Orientierung einer Masche andert, so kehren sich die Vor- 
zeichen aller auBerhalb der Hauptdiagonale liegenden Elemente der betreffenden Zeile und 
Spalte um. 


Die Transformation 


Zi G7 C (7) 


kann auch aus der Leistungsinvarianz beziiglich der Zweigstréme einerseits und den Maschen- 
strémen andererseits gefolgert werden. Die komplexe Leistung durch die Zweigstréme aus- 
gedriickt ist 


SPO ih =i Zit (9a) 
und durch die Maschenstréme ausgedriickt 
SaaS FAs (gb) 


Beide Leistungen miissen gleich sein, d.h. S = S®. Setzt man Gl. (2) in (ga) ein, so ergibt 
sich bei Vergleich mit (gb) 


ievh ene eZ 8 iors r? 1° (10) 
Nun folgt hieraus noch nicht die Beziehung (7), sondern erst daraus, daB die Maschen- 


stroéme 7), i2...4 (7 = z—k +1) linear unabhangig gewahlt werden diirfen. Wahlt man 
z. B. der Reihe nach die folgenden Maschenstromverteilungen 


ab @) 

(ae oO Oo il AN (0) 

iy _— A ’ Io =e . gt I, <7 . ’ (11) 
O O 1 


so 1aBt sich unter Zusammenfassung der Stromvektoren i°, i9, ... i$ zu einer Matrix |°, 


wobei 
° = 1 = Einheitsmatrix (12) 


ist, die Beziehung (7) herleiten.? 

Die Maschenstréme erhalt man eindeutig, falls die Determinante des Systems nicht ver- 
schwindet, bzw. gewisse Passivitatsbedingungen erfillt sind (vgl. Kap. 5.1), aus einem Glei- 
chungssystem mit y Unbekannten 


FLD = ae tS (13) 


Hieraus erhalt man die Zweigstréme aus Gl. (z) und die Zweigspannungen aus Gl. (4). 
Damit sind alle wesentlichen GréBen des Netzes ermittelt. 


3. Die Knotenpunktsmethode 


In dem zu untersuchenden Netz nach Bild 1 numerieren wir auch noch die Knoten- 
punkte (1...). Wir kénnen dann eine Knotenpunkt-Zweig-Inzidenzmatrix K aufstellen: 
In der x-ten Zeile und ¢-ten Spalte tragt man die Zahlen +1, —1 oder 0 ein je nachdem, ob 
der Orientierungspfeil des ¢-ten Zweigs auf den x-ten Knotenpunkt hinweist, wegweist oder 


1 Index ¢ bedeutet transponierte(r) Matrix (Spaltenvektor), Stern (*) konjugiert-komplexer Wert. S = 
komplexe Leistung, P = Wirkleistung, Q = Blindleistung. Die Betrage der 1-, bzw. u-Drehzeiger sind in 
Effektivwerten einzusetzen. 

2 Wegen der Deutung dieser verallgemeinerten Kongruenztransformation als Ausdruck der Leistungs- 


invarianz einer rechteckigen Ubertragermatrix vgl. [46]. 
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nicht mit ihm verbunden ist. In unserem Beispiel lautet die den orientierten Graph umkehr- 


bar eindeutig beschreibende Matrix 


¢= 1 2 3 4 5 6 
Sm 
(@) +1 +1 —1 O (o) 
= K. (14) 
2 +41 (0) fe) +41 +1 (6) 
(0) —1 (0) e) —1 —i1 
4 => il (0) == il () () +1 


In jeder Spalte kommt je einmal +1 und —1 vor. 


Die letzte Zeile ist gleich der negativen 


Summe aller iibrigen Zeilen, sie ist linear abhangig. Sie (oder irgendeine andere) kann bei der 
Aufstellung der notwendigen Gleichungen unberiicksichtigt bleiben. Bezeichnet man die durch 


Streichen der letzten Zeile entstandene Matrix mit K, so kann man die notwendigen Gleichun- 
gen fiir die Knotenpunktspotentiale folgendermaBen gewinnen: Fiir die Knotenpunktssummen- 


stréme gilt 


Tis + ty. + tg — tq . 
Lane n state ele . 
ae — ty —15;—1¢ 
ig SS us — ts + t6 
oder 
pee ie (15) 
Die Zweigspannungen sind darstellbar als Differenzen von Knotenpunktsspannungen 
= + ug — uf 
Uy, = + uf — uz 
us = + uf — uy, 
UW, = — ut + uf 
(Pie + uf — uf 
pee —uz + ut 
oder 
Ui Kus (16) 
Aus Gl. (4)ereibt sich, falls 2j7-25.40.. @.,.6= 0 
i = Z1(u —e) 
(17) 
= ¥ (u—e), 
wobei 
i 
—- fo) O 
411 
1 
0) — . fo) 
1 
(0) O ere a 


. die Zweigadmittanzmatrix ist. 


a ' 

+ Hierin unterscheiden sich die beiden Methoden: Wahrend bei der Maschenmethode verschwindende 
Impedanzen, jedoch keine verschwindende Admittanzen zugelassen sind, sind bei der Knotenpunktsmethode 
verschwindende Admittanzen, jedoch keine verschwindende Impedanzen zugelassen (wichtig fiir Grenzfalle 
wie Unterbrechungen und Kurzschliisse). In der Maschenmethode bedeutet eine verschwindende Admittanz / 
Wegfallen eines Zweiges, in der Knotenpunktsmethode verschwindende Impedanz Zusammenfallen zweier 
Knotenpunkte. In beiden Fallen miissen Eingriffe in die Topologie und damit in die Inzidenzmatrizen vor- 


genommen werden! 
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Die KircHHorrFschen Knotenpunktsgleichungen sind durch Gl. (15) gegeben, wenn man 
iO (19) 
setzt, die Maschengleichungen sind in Gl. (16) enthalten, wahrend Gl. (17) die Anwendung des 
ohmschen Gesetzes darstellt. Wir erhalten dann ausgehend von Gl. (15) durch Einsetzen von 
Gl. (16) und (17) die Beziehung 


t= KY K,ut+—KYe=0 (20) 
oder 


Key Keun =a Ye: (21) 


In diesem Gleichungssystem ist aber wegen der besonderen Eigenschaft von K die letzte 
Gleichung die negative Summe aller iibrigen und deswegen entbehrlich. Ferner kénnen wir, 
da die Knotenpunktspotentiale nur bis auf eine willkiirliche Konstante bestimmt sind, das 
Potential des letzten Knotenpunkts u, = 0 setzen (Bezugsknotenpunkt), damit entfallt in K 
die letzte Zeile und in K, die letzte Spalte. Das System, das jetzt keine tiberfliissigen Gleichun- 
gen mehr enthalt und eindeutig lésbar ist, falls gewisse Passivitatsbedingungen erfiillt sind, 
lautet ' 


KY ue hY¥e, (22) 
Neu Yeo (22a) 
In unserem Beispiel ergibt sich fiir K Y K, die folgende Matrix (mit — vi) 
244 
Yeo + Yaa + Yaa PERL Y — Yee : ay as 
KYK y+ aaa Yar + Van + V55 em-5 15 | => sam 
= = a (2 
: — Yee = 66 Yoo + Y55 + Yoo | Sa C6 (23) 
hy 8G aw At = Vee ¥11 + Vas + Yeo 


Sie ist ebenfalls symmetrisch (Beweis!) und auBerdem noch singular, da die Summe aller 
Spalten (Zeilen) null ergibt. Eine fiir Re y;; > o nichtsingulare Matrix erhalt man durch Strei- 
chen der letzten Zeile und Spalte. Hierdurch entsteht die durch die gestrichelten Geraden ab- 


geteilte Hauptuntermatrix (k—41)-ter Ordnung Yoo (Gy. K,. Auch diese Matrix hatte 


man ohne Bildung des Produkts K Y K, bzw. K Y K, direkt aus Bild 1 ablesen kénnen. Fir 
die einzelnen Elemente gelten, wie man leicht nachweist, die folgenden Regeln: 


1. In der Hauptdiagonale ist das Element y;; gleich der Summe aller Admittanzen, 
welche mit dem 7-ten Knotenpunkt verbunden sind. 

2. DasauBerhalb der Hauptdiagonale stehende Element y;, = yj; ist gleich der nega- 
tiven Summe aller Admittanzen, welche den 7-ten und k-ten Knotenpunkt miteinander 
verbinden. 


Die iibrigen Hauptuntermatrizen lassen sich als Y+ -Matrizen deuten, welche einem Netz 
entsprechen, in welchem die den gestrichenen Zeilen und Spalten entsprechenden Knoten- 
punkten mit dem Bezugsknotenpunkt identifiziert sind. 

Die Matrix (17) beschreibt das zugrundegelegte Netz eindeutig, d. h. man kann aus dieser 
Matrix das Netz jederzeit wieder eindeutig rekonstruieren. Die Information iiber die Orien- 
tierung der Zweige ist auch hier verlorengegangen. 

Aus Gleichungssystem (22) ergeben sich die Knotenpunktsspannungen bezogen auf den 
Bezugsknotenpunkt eindeutig, falls die Determinante des Systems nicht verschwindet. Hier- 
fiir ist das Erfiilltsein gewisser Passivitatsbedingungen hinreichend (vgl. Kap. 5.2). Aus Gl. 
(16) lassen sich dann alle Zweigspannungen ausrechnen und aus Gl. (17) alle Zweigstréme. 

Auch die Transformation 


Y=KYK, bews Ye=KYK, (23) 


426 H. Epetmann: Allgemeine Grundlagen der Netzberechnung mit Inzidenzmatrizen Cee 
a 


kann aus der Invarianz der Leistung analog dem entsprechenden Vorgang bei der Maschen- 
methode gefolgert werden. Hierzu denkt man sich treibende Spannungen U;* von auBen zu- 
gefiihrt, wahrend die inneren Spannungen verschwinden sollen. Dann ist 


SS UP IOY Kup un 0 (24) 

Nun darf man wieder k —1 linear unabhangige Vektoren uj--+-U%_,, auswadhlen zweck- 
maBigerweise wieder so, daB sich bei Zusammenfassung zu einer Matrix die Einheitsmatrix 
U+ = 1 (25) 


ergibt. Alsdann ergibt sich aus Gleichung (24) die Transformation Y+=KY K,. Da das 
Potential des Bezugsknotenpunktes ebenfalls frei wahlbar ist, kann man auch ein System von 
k linear unabhangigen Vektoren ut---u% so aufstellen, daB sich eine k-reihige Einheits- 
matrix U+ = 1 ergibt. Wendet man dieses System von Vektoren auf eine der Gleichung (24) 
entsprechende Beziehung fiir die absoluten Potentiale an, $0 folgt auch Yt = K Y K,,.? 


4. Besondere Eigenschaften der K- und C-Matrizen 


7 


4.1 Ermittlung eines Baumes und der unabhangigen Zweige aus der K-Matrix 


Wir betrachten eine bereits lexikografisch geordnete K-Matrix, d.h. eine K-Matrix, die 
durch Umordnung der Spalten (= Vertauschung der Zweignummern) derart, daB die Indizes 
der Knotenpunktsnummern yp, » mit w < y lexikografisch angeordnet sind. Man kénnte auch 
von einer nicht geordneten Matrix ausgehen, doch erscheint das entsprechende Diagramm 
iibersichtlicher, wenn die Matrix lexikografisch geordnet ist. Anstelle der GréBen +1 oder 
—1 denke man sich zunachst entsprechend einem nichtorientierten Graphen +1 gesetzt oder 
noch einfacher: Kreuze (x). Die Orientierung soll spater eingefiithrt werden und zwar grund- 
sdtzlich mit w< y von w nach v, so daB der Pfeil immer von einer kleineren Knotenpunkts- 
nummer zu einer gréBeren zeigt. In der K-Matrix hat das zur Folge, daB spater anstelle der 
beiden in jeder Spalte auftretenden Kreuze oben eine ,,—1“ und unten eine ,,+1“‘ zu setzen 
ist. Hierzu soll die konstruktive Definition des Begriffs Baum in einem Netz herangezogen 
werden: 

Man ermittelt einen vollstandigen Baum in einem gegebenen Netz dadurch, daB man von 
einem beliebigen Knotenpunkt ausgehend immer nur solche Zweige des Netzes hinzunimmt, 
die eine Verbindung zu einem neuen Knotenpunkt herstellen, bis alle Punkte miteinander 
verbunden sind, d. h. jede Hinzunahme eines weiteren Zweiges mu8 auch die Anzahl der Knoten- 
punkte um eins erhohen. 

M. a. W. um einen Baum zu ermitteln kann man ein Belegungssystem fiir die Knoten- 
punkte anwenden. Bei jeder Neubelegung eines Punktes ist darauf zu achten, daB ein bereits 
belegter Punkt mit einem bisher noch unbelegten Punkt durch einen im Netz vorkommenden 
Zweig, der dann Baumzweig wird, verbunden wird. Die iibrigen Zweige sind Verbindungs- 
zweige, auch unabhangige Zweige genannt. Um alle Zweige vollstandig zu erfassen, geniigt es 
darauf zu achten, da® der Reihe nach fiir jeden Knotenpunkt saémtliche von ihm ausgehenden 
Zweige daraufhin untersucht werden, ob sie zu einem nocht nicht belegten Knotenpunkt 
fihren und also zu den Baumzweigen zu rechnen sind oder zu den unabhangigen Zweigen. 
Sind samtliche Zweige eines Knotenpunkts erfaBt, so wird dieser abgehakt. 


Beispiel 2: 


Als Beispiel fiir eine abstrakte Ermittlung eines Baumes, d.h. ohne Zuhilfenahme der 
zugrundegelegten Figur (Bild 2) des Streckenkomplexes, aus der K-Matrix soll ein etwas kom- 
plizierteres Netzwerk betrachtet werden, das auBerdem auch Parallelzweige (8 u. g) enthilt. 


1 Auch diese verallgemeinerte Kongruenztransformation kann durch eine Ubertragermatrix [46] anschau- 
lich dargestellt werden. 
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Zum Aufbau des Baumes beginnen wir z. B. mit dem Knotenpunkt 1 (1. Zeile) und be- 
legen sdmtliche mit dem Knotenpunkt 1 tiber einen Zweig verbundenen Knotenpunkte ein- 
fach, d.h. sind Mehrfachverbindungen vorhanden, so belegt die 1. Verbindung den betreffen- 
den Knotenpunkt, die iibrigen Verbindungen erscheinen als unabhangige Zweige. In dem 


1} —1 —1 —1 fo) ° ° fe) fo) 
2 +1 fo) o —1 —1 fo) ° ° ° fe) ° 
3 o +1 o +1 o —1i —1 fo) fo) fo} fo) 
b 4 ° fo) oO ° OF 1 o —1 —1 —1 —1 0 
5 fo) o +1 fe) fo) o +1 +1 +1 fo) Onno: 
6 fo} fo) {0} o +1 fo} fo} fe) o +1 o—1 
GA ° fe) fo) ° 0) ° ro) ° ° o +1 +41 


Bild 2. a Zur Ermittlung eines vollstandigen Baumes. Orientierungspfeile weisen von einer 
niedrigeren zur hOheren Knotenpunktsnummer. b Zugehérige K-Matrix. Entsprechend der hier 
gewahlten Orientierung steht eine ,,—1“‘ tiber einer ,, 1“. 


Schema (Bild 3) erfolgt die Belegung durch kleine Kreise (links neben der Nummer). In un- © 
serem Beispiel werden auf diese Weise durch die Zweige 1, 2 und 3 die Knotenpunkte 2, 3 und 
5 belegt. Knotenpunkt 1 als Ausgangsknotenpunkt gilt auch als belegt. Damit gilt Knoten- 
punkt 1 als erledigt, gekennzeichnet durch Abhaken rechts. 


Cae ty) 4 


cog CER Se AA OG Em We ESS 


Bild 3. Belegungs- und Abhakschema unter Verwendung der nichtorientierten 
K-Matrix fiir den Streckenkomplex (Graph) nach Bild 2a. 


Als nachster Knotenpunkt kann nur ein bereits belegter Knotenpunkt gewahlt werden. 
Der nachstliegende ist Knotenpunkt 2. Die mit ihm verbundenen Zweige 4 und 5 fiihren auf 
die Knotenpunkte 3 und 6, die zu belegen sind. Punkt 3 ist schon belegt und Zweig 4 ergibt 
daher einen unabhangigen Zweig. Alsdann kann Punkt 2 abgehakt werden. 

Der niachste bereits belegte Knotenpunkt ist Punkt 3. Von ihm aus kann iiber Zweig 6 
der noch nicht belegte Punkt 4 erreicht werden (Baumzweig), wahrend tiber Zweig 7 der 
bereits belegte Punkt 5 erreicht wird (unabhangiger Zweig). Nachdem also Punkt 4 belegt ist, 
kann Knotenpunkt 3 abgehakt werden. 

Wir kénnen nun zu Knotenpunkt 4 iibergehen, der soeben belegt wurde. Die Parallel- 
zweige 8 und 9 fiithren auf Punkt 5, wahrend Zweig 10 auf Punkt 6 fiihrt. Da die Punkte 5 und 
6 belegt sind, ergeben sich unabhangige Zweige. Zweig 11 hingegen fiihrt auf einen noch nicht 
belegten Punkt 7, der jetzt als letzter belegt wird. Zweig 11 ist somit auch der letzte Baum- 
zweig. Damit kann Knotenpunkt 4 abgehakt werden. Knotenpunkt 5 ist schon erledigt und 
kann auch abgehakt werden. 

Knotenpunkt 6 liefert nur einen abhangigen Zweig 12. Alsdann kann Punkt 6 und schlieB- 
lich auch Punkt 7 abgehakt werden. Alle Knotenpunkte sind belegt und abgehakt, somit sind 
auch alle Zweige erledigt, das Verfahren ist beendet. 
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Zahlentafel 1 


Ermittlung eines Baumes aus der Liste der Zweige. Ubergang zu einer Nor- 
malnumerierung zur Erzeugung einer Dreiecks-Baummatrix 


Ausfiihrliche Schreibweise Kurzschreibweise 
| aufi Vorlaufi Unabhangiger ‘ Neue 
Sa oa Reval enpunicts: aga - Neue Zweig-| knotenpunkts- 
nummer nummer oder Baum Bienen nummern 
ee ee SS eet 
1 yt 2 B, Gf ete 
2 1u.3 B, c= Pat 93) 
3 1u.5 Bs 9 — 1u.4 
4 Zales U, 1 = 7505-3) 
5 2u.6 B, 10 — 25 
: 6 verbindet die 3uU. 4 B; i == 3u.6 
Zweig 
Gy Knotenpunkte 3u.5 U, 2 — 3u.4 
8 4uU. 5 Us 3 — 4u. 6* 
9 4uU. 5 U, 4 — 4 u. 6* 
10 4u.6 U; 5 — 5 u. 6* 
11 Ata 7 Bg 12 — One7, 
12 J rues? U; 6 — Gyles 9 
* Reihenfolge der Knotenpunkte vertauscht (Umorientierung). 
Zahlentafel1ia 
Hilfstabellezum BelegenundAbhakender Knotenpunkte 
Belegung ( ) Abhaken Endgiiltige Belegung ( ) Abhaken Endgiiltige 
(Vorl. Numerierung) (Vv) Numerierung \Vorl. Numerierung) (Vv) Numerierung 
1 Vv 1 6 V 5 
2 Vv 2 4 Vv 6 
3 Vv 3 7 Vv 7 
5 Vv 4 


Aus dieser an einem besonderen Beispiel gezeigten Methode zur Ermittlung eines Baumes 
und der unabhangigen Zweige aus der K-Matrix kénnen jetzt leicht gewisse allgemeingiiltige Be- 
ziehungen hergeleitet werden. 

Der gleiche ProzeB laBt sich auch an Hand der Liste aller Zweige und einer Hilfstabelle 
durchfiihren (vgl. ByERLY, Lone u. Kinc [56]). Wie man in Zahlentafel 1 und 1a erkennt, 
macht es keine Schwierigkeiten, den Baum auch aus der Liste der Zweige zu ermitteln. 


4.2 Satze titber Unterdeterminaten der K-Matrix 


Ordnet man die Zeilen der K-Matrix in der Reihenfolge der Knotenpunktsbelegung und 
die Spalten der Baumzweige ebenfalls in der Reihenfolge der zugeordneten Zweigbelegungen 
an, so entsteht bei Streichen der ersten dem Ausgangspunkt zugeordneten Zeile eine obere 
Dreiecksmatrix mit Einsen in der Hauptdiagonale. Wahlt man die Orientierung der Zweige 
auBerdem auch noch in Richtung der Belegung, so daB in jeder Spalte iiber der Haupt- 
diagonale eine —1 und in der Hauptdiagonale eine +1 steht, so ergibt die Determinante +1. 
Vergleiche hierzu den umnumerierten und umorientierten Streckenkomplex in Bild 4. Wir 
k6nnen also den Satz aussprechen: 
| Wahlt man in einem Netz die Orientierung der Baumzweige in Belegungsrichtung 

und die Knotenpunkts- und Zweignumerierung in der Reihenfolge der Belegung (Normal- 

| numerierung), so ergibt die (k — 1)-reihige Baumunterdeterminante, die durch Streichen 

| der Zeile des Ausgangspunkts und der den unabhangigen Zweigen zugeordneten Spalten 
entsteht, den Wert + 1. 

Da im allgemeinen diese Ordnungsbedingungen nicht erfiillt sind, so kann sich der Wert 
der Baumdeterminante in anderen Fallen nur durch das Vorzeichen unterscheiden. Der Wert 
jeder Baumdeterminante auch von Teilbaumen ist also +1. Es ist hierbei auch gleichgiiltig, 
welche Zeile man einem Bezugsknotenpunkt zuordnet und streicht. Will man das Vorzeichen 
bestimmen, so denke man sich den ausgewahlten Baum schrittweise in der angegebenen Weise 
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aufgebaut. Hierbei ergibt sich bei Berticksichtigung des gewahlten Bezugsknotenpunkts eine 
gewisse Ordnung der Knotenpunkte und der Zweige und auBerdem eine Orientierung der Zweige. 


1 {e) (o) fe) (o) fe) OR a ed (o) (e) {o) 
Bp |-i fo) fo) +1 fe) @ sl ce) fe) 
2) ek fo} fe) o +1 fo) o —1 fo} 
De, Oy J SS Sn O @ bi ° Ono 
5 fc) fo) ° on 1 >—1, 0 fo) al Q@  ) 
6 ° o +1 +1 41 Gee fo) ° o +1 —1 
7 ° ° roy ° Cie Te eh i) ) o-+1 


Bild 4.a Umorientierter Streckenkomplex nach Bild 2a zur Herstellung einer Normalnumerie- 

rung der Knotenpunkte und Zweige einer Normalorientierung. b Zugehdrige K-Matrix. Ge- 

strichelt: Aufteilung in den Untermatrizen A und B. B erscheint infolge der Normalnumerierung 

als untere Dreiecksmatrix und wegen der Normalorientierung sind samtliche Hauptdiagonal- 
elemente gleich +1. 


“Ist die Anzahl der erforderlichen Zeilen- d. h. Knotenpunktsumordnungen a, der Spalten- 
d. h. Zweigumordnungen f und die Anzahl der Umorientierungen y so ist das Vorzeichen und 
damit auch der Determinantenwert bestimmt durch 

det B = (—1)*t8+?, (26) 
Man kann nun auch Unterdeterminanten mit Zweigen bilden, welche Zyklen enthalten. 
Die Spaltenvektoren von Zyklen sind aber linear abhangig und somit verschwinden solche 
Determinanten. Der Rang der Matrix B ist also genau k —1. Es gilt also der Satz: 

Die Untermatrizen der vollstandigen Baume sind Rangmatrizen. 

Durch Fortnahme von Zweigen aus einem Baum verbleiben einer oder mehrere Unterbaume. 
Einem solchen Unterbaum mit k’ Knotenpunkten und k’ —1 Zweigen kann man bei gleich- 
zeitiger Auswahl eines Bezugsknotenpunkts eine k’—1-reihige Unterdeterminante von K zu- 
ordnen. Diese hat den Wert +1. Alle iibrigen Unterdeterminanten, d.h. solche, die nicht 
einem Unterbaum entsprechen, verschwinden. Hieraus ergibt sich der folgende Satz: 

| Die Unterdeterminanten der Matrix K kénnen nur die Werteo, +1 oder —1 annehmen. 

| Die nichtverschwindenden Unterdeterminanten entsprechen Baumen bzw. Unterbaumen. 


Matrizen, deren Elemente und Unterdeterminanten nur die Werte 0, +1 oder —1 anneh- 
men kénnen, werden E-Matrizen genannt und sind hauptsadchlich von VEBLEN [3] und CEDER- 
BAUM [51] untersucht worden. 


4.3 Zusammenhang zwischen K- und C-Matrizen 


Hat man in K eine spezielle Baummatrix ausgewahlt, so kann jetzt gezeigt werden, daB 
die hierzu gehérige C-Matrix auch durch Rechnung aus der K-Matrix gefunden werden kann. 
Wir wahlen einen willkiirlichen Bezugsknotenpunkt und streichen aus der K-Matrix die zu- 
gehorige (linear abhangige) Zeile. Ferner ordnen wir die darin enthaltene quadratische Baum- 
matrix B rechts an, die verbleibende den abhangigen Zweigen zugeordnete Untermatrix A 
links, so daB die reduzierte Matrix mit diesen Untermatrizen geschrieben werden kann 


Ke AB) (27) 
Nach Gl. (9) und (14) gilt nicht nur 
i= Ki =-0" (28) 
sondern auch eat 
K i= A B) i= 0: (29) 


Die in den ersten ry = z—k + 1, namlich den unabhangigen Zweigen flieBenden Stréme sind 
zugleich die Maschenstréme mit dem Stromvektor i°, die iibrigen, in den Baumzweigen flie- 
Benden Stréme seien in dem Vektor i’ zusammengefaBt. Damit l4Bt sich Gl. (29) schreiben 


(A8)-(S)= Ae + Bi =O. G0) 


30m 
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Diese kann wegen det B = +1 40 aufgeldst werden nach den Baumstrémen i’. Es ist 

ii =—B'A??. (31) 
Nun kann auch eine Beziehung zwischen dem vollstandigen Vektor 7 und demjenigen der 
Maschenstréme aufgestellt werden. Es ist 


a () ig & a7 d z e a) ae ) 


Da diese Matrix eindeutig ist, ergibt ein Vergleich mit Gl. (2) 


1 
C=(_ peal: (33) 


Hierbei ist es nicht erforderlich B-! tatsidchlich vorher auszurechnen, sondern es geniigt, das 
Gleichungssystem (31) in der Form 


Bi.=—Ai° (34) 


aufzulésen, was wiederum keine Schwierigkeiten macht, wenn die Baummatrix die Form 
einer oberen Dreiecksmatrix hat!. Numeriert man also die Knotenpunkte und Zweige in der 
Reihenfolge, in welcher der Baum aufgebaut wird (Normalnumerierung) und wahlt man die 
Orientierung der Baumzweige in der Weise, daB der Pfeil von einem Punkt niederer zu einem 
héherer Nummer zeigt (diese sog. Normalorientierung ist zwar nicht notwendig, vereinfacht 
aber die Vorzeichenfestsetzung wesentlich), so kann das Gleichungssystem von unten her 
aufgerollt werden. Das Gleichungssystem unseres Beispiels lautet dann namlich folgender- 
maBen 
ty — to =+%4 
ts has Sk te 
ty ar teeta V4 
416 = + 15 + % 
yy — ty, = =a tg ae 

a Bi 
Dieses Gleichungssystem ist nach den Unbekannten 77...1,. aufzulésen. Die Lésung fiir 1,, 
steht bereits da. In der vorletzten Gleichung kann man i,, auf die rechte Seite bringen und 
einsetzen usw. In jeder Gleichung treten auBer der Unbekannten, nach der aufzulésen ist, 


nur bereits ermittelte Unbekannte auf (eine Folge der Dreiecksgestalt des Systems). Als 
Lésung ergibt sich 


ty = +h 4y.sbt¢ 
1g = —4 +1, —1, — 144 — 15 — 1 
1 a ty tg ty (36) 
to a + ts + % 
ty =a ye Upp ae bg 
Lal Pee menace 
oder 
1 ) O (0) +1 +1 
—1 oe —1 —1 —1 —i1 
Bs Ano Wee ae ‘ eae (37) 
) ) O (0) + 1 +1 
1) O alk == il —— il a 
(@) O O 10) O =i ‘ 


* Offenbar ist es in einer maschinenmaBigen Behandlung emnpfehlenswert, die K-Matrix sogleich in ihrer 
Normalform, d.h. unter Zugrundelegung einer Normalnumerierung der Zweige und Knotenpunkte, einzu- 
fiihren und die Vorarbeit: Festlegen eines Baumes und Umnumerierung anhand der Liste aller Zweige durch- 
zufiihren, damit die wegen der relativ groBen Zahl enthaltener Nullen nur geringe Informationsdichte auf- 
weisende K-Matrix nicht zweimal die Speicherplatze belegt. 
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und die C-Matrix lautet somit 


+1 (6) O O O 
O +1 O ( (6) (6) 
fe) (0) + 1 (0) ) 
fe) 0) (6) +1 (e) O 
(e) O (e) +1 ) 
eS ee 
22 Val +41 ) 0 Oho? gear +41 
Z —1 +41 —1 —1 —1 —1 
: fo) —1 +41 +41 O fe) 
fo) (e) O (0) +1 +1 
0 O —1 —i1 —1 —1 
O O 0) 0) ) —i1 


Natiirlich 14Bt sich diese C-Matrix, wie schon gezeigt, auch direkt aus dem Netzwerk ablesen. 
Wenn wir aber diese Aufgabe einem digitalen Rechenautomaten iiberlassen wollen, ist es not- 
wendig, bei Vorgabe der K-Matrix diese Aufgabe rein logisch und algebraisch zu bewaltigen. 
Die Lésung der logischen Aufgabe ist durch die Konstruktion des Baumes beschrieben. Die 
Realisierung hangt dabei noch von gewissen Gegebenheiten der benutzten Maschine ab. Nach 
der ggf. noch erforderlichen Umnumerierung bedeutet die geschilderte Aufrollung des dreiecks- 
formig gestaffelten Gleichungssystems keine Schwierigkeiten mehr. Damit ist ein weiteres 
gegentiber der Arbeit von ByerLy, Lone und Krnc [56] unabhangiges Verfahren zur Bildung 
der C-Matrix skizziert. 


4.4 Herleitung der Gleichung KC =0 


Aus Gl. (33) 1aBt sich noch eine weitere wichtige Beziehung herleiten. Es ist namlich 


Kc=(AB)(_ 9, |-A-A=0. (39) 
Fiigt man in K die gestrichene Zeile — nehmen wir an, es sei die letzte (k-te) Zeile — 
(toe ek! (40) 
so 1aBt sich die Gesamtmatrix darstellen in der Form 
K = ( : . (41) 
(—1—1---—1)K 
Ose 520 
0) £ 0 
se Ne | KK (42) 
fe) Onn ee 
—1 —1 —1 


und die entsprechende Beziehung fiir die Gesamtmatrix lautet 


OF. O 
O ieee O 
KC= : 5 OG eau (43) 
O OR. 1 
Sb ea el 


1 Es ist namlich jede k-te Zeile, welche die K-Matrix zur vollen Matrix K erganzt, gleich der negativen 
Summe aller tibrigen Zeilen. 
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Diese Beziehung 14Bt sich natiirlich auch direkt herleiten. Mit den Gln. (2), (15) und (19) 


it=Ki=0 (15), (19) 
i= Ci? (2) 
ergibt sich at ee ae (44) 


Nun folgt aus dieser Beziehung noch nicht ohne weiteres K C = 0, sondern erst in Verbindung 
mit der Tatsache, daB die Maschenstréme linear unabhangig gewahlt werden diirfen. Setzt 
man demzufolge dann einmal i?, = 1 und die iibrigen Maschenstréme gleich null, dann ee 
und wieder die iibrigen Maschenstréme gleich null u.s.f. Der hierdurch entstehende Satz 
von Gleichungssystemen l4B8t sich mit Hilfe dieser i¥-Spalten-Matrizen 


/1 ) O 

A O A al A 10) 
n=. |, ih=[V]e- =]. (45). 

O 10) il 


ub Ome O 
TOE A Lue) 
eA ba ape ae : ee aoe (46) 
(0) Or. at: 
schreiben KCl@=-KC=0.1 (47) 


Man kann diese Beziehung auch aus entsprechenden Spannungsgleichungen erhalten. Hier 
wird die Tatsache benutzt, daB die Summe aller Spannungen in einer Maschine gleich null ist 
und daB die Knotenpunktspotentiale linear unabhangig gewahlt werden kénnen. Analog zu 
den bereits gemachten Uberlegungen ergibt sich aus den Gln. (3), (5) und (16) 

Uareon GU ged) (48) 


und wegen der Linear-Unabhiangigkeit der wt-Komponenten 


1 LONGO 

; aerate 8 
CK Ure Gis : ae (49) 

Or tns.On ches od 
oder auch =O Nes (50a) 
KC=0. (50b) 


Diese Gleichung kann man als Probe fiir die richtige Bildung von C aus K verwenden. 


4.5 Ermittlung der vollstandigen Liste aller Baume aus det Yt 


Bildet man det Yt, so ergibt sich eine algebraische Summe von Produkten der Zweig- 
admittanzen mit k —1 Faktoren entsprechend der Definition einer Determinante. Man kann 
nun zeigen, daB nach Wegstreichen aller sich aufhebenden Glieder nur noch solche mit posi- 
tivem Vorzeichen tibrigbleiben und nur diejenigen Kombinationen, welche einem Baum ent- 
sprechen®. Hierbei kommt jeder Baum einmal vor. Man kann also durch Ausrechnen dieser 
Determinante eine vollstandige Liste aller Baume finden. Der Beweis dieses Satzes beruht 
auf der Tatsache, daB von den (k —1)-reihigen Unterdeterminanten von K nur die Baum- 
unterdeterminanten nicht verschwinden und gleich +1 sind. Nach einem Satz von CAUCHY 
und BINET (vgl. z. B. A. C. AITKEN [9] S. 86, H. DOrRIE [11] S. 45 oder W. GROBER [43] S. 101) 
ist die Determinante eines Produkts einer Matrix F von m Zeilen und m Spalten mit einer 
Matrix G von n Zeilen und m Spalten, hierbei sei m < n, gleich der Summe von Determinanten- 


1 Vgl. hierzu den analogen Beweis bei W. CavER [12], 2. Aufl., S. 70. 

* Auf diesen implizit in der kombinatorischen Regel von G. KircHuoFF [1] zur Ermittlung der Zweig- 
strome enthaltenen und offenbar véllig in Vergessenheit geratenen Satz haben verschiedene Autoren ausdriick- 
lich hingewiesen [14], [22], [26], [29], [30], [31], [35], [36], [50] und [56]. 
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produkten von Unterdeterminanten maximal méglicher Reihenzahl n x n (sog. Maioren) aus 
den beiden Matrizen; hierbei sind alle méglichen Spaltenkombinationen aus F und Zeilen- 
kombinationen aus G durchzuvariieren und in homologer Anordnung miteinander zu multi- 
plizieren. D.h. wenn H = F- G und ferner F{ und G% die n-reihigen Unterdeterminanten 
und jeweils zueinander homologen Spalten- ‘bzw. Zeilenkombinationen variiert tiber allen 


ie eee ("" moglichen Kombinationen, so laBt sich die Determinante des Produkts schreiben 
det H = > F™ = G®)., (51) 
a=1 
Im Falle m = n reduziert sich die Summe auf ein Glied, und es ist 
det H = det F- detG. (51a) 
Hingegen mu8 H im Falle m > n singular werden, und es ist : 
det Ho. (51b) 


Da nun aber 


so kann man die Diagonalmatrix Y mit einem linken oder rechten Faktor zusammenfassen. 
Nehmen wir an, wir haben K Y gebildet, so erscheint sie als eine spaltenweise mit den zugeord- 


neten Zweigadmittanzen durchmultiplizierte 
a2 NS v K-Matrix. 
Die nichtverschwindenden (k — 1)-reihi- 
gen Unterdeterminanten bestehen dann aus 
den Produkten der Zweigadmittanzen eines 


Baumes und haben positives oder negatives 
Vorzeichen. Das gleiche Vorzeichen und den 
Betrag 1 hat die homologe Baumunterdeter- 
minante von K,. Dies ist aber gerade die- 


= jenige, mit der die Unterdeterminante aus KY 
multipliziert werden muB. Bei den nach 


| | { uy CAUCH Y-BINET (25] durchzuvariierenden 


Spalten bzw. Zeilenkombinationen kommen 
somit alle méglichen Kombinationen zu 
Baumen als positive Produkte der beteiligten 
Zweigadmittanzen in dieser Summe vor. Es 
| 4 i A, ist also 
LS UY cee eV yrs ear a a ey ee oe 
366 456 Oe einen 


Bild 5. ee es aller moglichen vollstandigen Baume aller Baumzweige (53) 
der WHEATSTONE-Briicke. 


Fiir das Beispiel der WHEATSTONEschen Briicke nach Bild 1 ergibt sich 


Yoo + 33 + Vaa ay ad — Vee | 
det yt ae rien Via a Vea 55 Oda 
— vee set: UE ar od Aa Yel 


= + Vi1 Von Vas + Vir Yoo Yaa + Vir Yoo V55 + Vir Yoo Yoo 
+ ¥11¥s3 Vss + Via Yas Yon + Vii Yaa Y55 + V11 Yaa Yoo 
+ Vo2¥33 Yaa + Yoo Yaa ¥s5 + Yoo Yaa Yoo + Yao V55 Veo 
+ Y33 Yaa Vs5 + Yas Yaa Vou + V33 Y55 You + Yaa V55 Veo: 


Bild 5 zeigt die sich hieraus ergebenden Baumkonfigurationen. 


1 Bei Determinanten héherer Ordnung empfiehlt es sich kombinatorisch vorzugehen, vgl.z. B.S. J. MAson [45], 
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Auch die Gesamtzahl aller méglichen Baume 1aBt sich auch aus einer Determinante er- 
rechnen. Man muB nur saimtliche Zweigadmittanzen y,, gleich 1 setzen. Dann ergibt sich fur 
jeden méglichen Baum eine ,,1‘ als Determinante. Da alle Terme positives Vorzeichen 
haben, ergibt die Gleichung fiir die Gesamtzahl aller méglichen Baume 


= det K1 K, = det K K,. (55) 


5. Hinreichende Bedingungen fiir Auflésbarkeit der Maschen- 
und Knotenpunktsgleichungen 

Die Frage nach der eindeutigen Auflésbarkeit der Gleichungssysteme (13) und (22) ist nicht 
nur von akademischer Bedeutung, sondern sie hat in dem Falle, daB das Gleichungssystem 
von einer programmgesteuerten Maschine gerechnet wird, die sehr reale Bedeutung, da8 man 
nur in dem Falle, daB das Gleichungssystem eindeutig auflésbar ist, von vornherein sicher 
ist, daB die Maschine in allen Fallen ohne Blockierung durch Charakteristikiiberlauf (z. B. bei 
Division durch o aber auch durch sehr kleine Zahlen) den Algorithmus zu Ende rechnet. 


5.1 Eindeutige Auflésbarkeit der Maschengleichungen 


Das System (13) der Maschengleichungen, das aus ebenso viel, namlich r Gleichungen be- 
steht wie unbekannte Maschenstréme auftreten, ist eindeutig auflésbar, wenn die Determinante 
des Systems nicht verschwindet. In diesem Fall hat auch das homogene System 


Z0 j0 — 0 (56) 


nur die Triviallésungen 7?...77 = 0. Das bedeutet aber physikalisch, daB ohne treibende 
Spannungen stationdre Maschenstréme nicht bestehen kénnen. Wenn z. B. in jedem unab- 
hangigen Zweig ein positiver ohmscher Widerstand enthalten ist, so wiirde dies zu einer 
dauernden Energieabgabe fiihren; die Lésungen des entsprechenden homogenen Gleichungs- 
systems kénnen also nur gleich null sein. Es gilt somit der folgende Satz: 


Damit das System der Maschengleichungen 
Z° ie =—C,e (13) 
eindeutige Lésung hat, ist hinreichend, daB in dem zugrundegelegten System un- 


abhangiger Zweige samtliche Impedanzen positive Realteile besitzen (und die 
| Baumzweig-Impedanzen keine negativen Realteile)?. 


Der diesem physikalischen Sachverhalt entsprechende algebraische Beweis verlauft 
folgendermaBen: Man zeigt, daB der hermitesche Teil der Matrix namlich 


Z% = —(Z° + Z)), (57) 


der tibrigens wegen der Symmetrie von Z reell ist, unter den gemachten Voraussetzungen 
positiv-definit ist und also nur negativ reelle Eigenwerte hat. Der kleinste Eigenwert von Z?, 
SEl Myyin, Aer groBte w,,,,. Dann gilt nach einem Satz von A. Hirscu? fiir die Eigenwerte A; von 
Z° die Ungleichung 


Ewin = Re A; = Emax * (58) 


Wenn nun aber die Eigenwerte w,; von Z?, samtlich negativ reell sind, so haben auch die Eigen- 


werte von Z° sdmtlich negative Realteile. Keiner der Eigenwerte kann dann verschwinden, 
und es ist 


det 2° = DA, xa. (59) 


+ Einen Beweis fiir ein Netz mit ausschlieBlich ohmschen Widerstanden (reeller Sonderfall) findet man bei 
HERMANN WEYL [2]. 


* Sur les racines d’une équation fondamentale. Acta mathematica Bd. 25 (1902) S. 367—370. 
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Wir haben nun noch zu zeigen, daB Z, positiv definit ist. Zu diesem Zweck gehen wir aus 
von der gesamten im Netz verbrauchten Wirkleistung 


P = Re i?* Z? i? = — (i9* Z° i9 + 19 Z°* i*) , 


2 
We 5 : ; 
Se (he Zeya oi")s (60) 
= if? Z4 iP, 
ral bag Ors a Oe be 
Fiir symmetrische, insbesondere auch fiir Diagonalmatrizen gilt aber 
"22 
Lan Re ZR ne . (61) 
x 
0) Br 
Es ist also 
j irra bade Sad ee Gay 
= ap Ri 
z . 
Ey Pag \tals (62) 
w=1 
Y Zz r 2 
= Era liP + ES tan|E ue tl (63) 
“u=1 w=rt1 \y=1 | 
=O 


Der erste Summenterm gibt die gesamte in den unabhangigen Zweigen, der zweite Summen- 
term die gesamte in den Baumzweigen verbrauchte und in Warme umgesetzte Energie an. 
Sind in den unabhangigen Zweigen die ohmschen Widerstande samtlich positiv, so ver- 


a 
schwindet P nur dann, wenn ¥ |i)/? = 0. Damit ist P eine positiv definite hermitesche 
u=1 
Form in den Variablen 7{, if... 7? w. z. b. w. 
Man erkennt auch, daB es wesentlich war, die Positivitat der Realteile in den unab- 


hangigen Zweigen zu fordern. Die Positivitat der Realteile der Baumzweigimpedanzen 


r 
ergibt bei > |z,|? > o nur Nichtnegativitat von P, d.h. Semidefinitheit der entspre- 
=1 


chenden hermiteschen Form in den Veranderlichen 1, (u=1...7), da sich in dem zweiten 
Summenterm die Veranderlichen (Baumstréme) 
tee Cayity fh er Ata 2 (64) 
v=1 


zu null erganzen kénnen. Das bedeutet, daB iiber den be- 
treffenden Baumzweig kein Strom flie8t und trotz Positivitat 
von 7,,, kein positiver Beitrag zu P geliefert wird. 

Man kann sogar Netze angeben, in welchen trotz Positivitat 
der Realteile samtlicher Baumzweigimpedanzen das homogene 
System 


Zr i= oO (56) 

i,t eine nichttriviale Lésung hat und also det Z° = Oist. Betrachten 

Bild 6. Ein Netz ohne Spannungsquellen, Wir z. B. einen sternférmigen Baum mit verlustbehafteten 

in ee een eae ee ls Zweigen (Bild 6). Die unabhangigen Zweige, welche die End- 

pe eee punkte sdmtlicher Baumzweige verbinden, seien widerstands- 

los. Dann bilden die unabhangigen Zweige einen geschlossenen Weg. Langs dieses Weges kann 
ein dauernder Kreisstrom flieBen. In den Baumzweigen flieBt kein Strom. Obwohl 

i 
25 ful = [eal 0 (65) 


u=1 
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ist nun P = 0, weil einerseits diez, flirw—=1...7 verschwinden und andererseits die Baum- 
stréme verschwinden. Die quadratische Form ist semidefinit. Da es offenbar nur darauf an- 
kommt, daB die homogene Gleichung keine nichtverschwindenden Lésungen hat, -nichtver- 
schwindende Lésungen aber auch dadurch vermieden werden, daB ein anderes als das fiir die 
Bildung von Z° zugrundegelegte System unabhangiger Zweige Impedanzen mit ausschlieB- 
lich positiven Realteilen besitzt, kann der bewiesene Satz folgendermaBen erweitert werden: 


| Damit das System der Maschengleichungen 
Zilie Gre (13) 


| eindeutige Lésungen hat, ist hinreichend, daB inirgendeinemSystem unabhangiger 
| Zweige sdmtliche Impedanzen positive Realteile besitzen (und die Baumzweige 
| keine negativen Realteile). 


In gewissen Fallen ist nun auch noch interessant, den Fall der um die r-dimensionale 
Mannigfaltigkeit mehrdeutigen Lésungen zu untersuchen, insbesondere den Fall, in welchen 
die den freien Variablen entsprechenden Kreisstréme als verschwindend angesetzt werden. 
In der Praxis sind zwar immer geniigend Verluste vorhanden, welche ein dauerndes FlieBen 
von Kreisstromen verhindern. In der Rechnung jedoch beriicksichtigt man diese Verluste 
aus Griinden der Einfachheit nicht, und es ergibt sich auch theoretisch entweder eine Mehr- 
deutigkeit, und im Falle von Rundungsfehlern in der Berechnung der rechten Seite des Glei- 
chungssystems kommt es sogar zu unvertraglichen Gleichungen. Im letzteren Fall wiirde also 
keine Lésung existieren, d.h. man wiirde formal ftir einige Unbekannte den Wert oo er- 
halten, und praktisch hatte man im Verlauf des Rechenprozesses entweder durch null, oder in- 
folge von Rundungsfehlern durch sehr kleine Zahlen zu dividieren. Im Fall eines Charakteri- 
stik-Uberlaufs wiirde die Maschine stoppen, hingegen im Falle, daB dieser Uberlauf nicht zu- 
standekommt, hielte man eine durch die von den Zwangsbedingungen erzeugten groBen 
Stréme véllig unbrauchbare Lésung. Es soll hier aber nicht untersucht werden, welche Vor- 
sichtsmaBregeln man in solchen Fallen treffen mu8, um zu verniinftigen Lésungen zu kom- 
men. Gliicklicherweise hat man es normalerweise nur mit dem Regelfall, in welchem der 
Héchstrang von Z® von vornherein gesichert ist, zu tun. 


5.2 Eindeutige Auflésbarkeit der Knotenpunktsgleichungen 


Damit das System der Knotenpunktsgleichungen 
Y¥rou=skK ye (22) 


mit k — 1 Unbekannten und ebensoviel Gleichungen eindeutig auflésbar ist, muB also auch 
hier die Determinante des Systems 


det Y+ 4 0 (66) 
sein. In Kap. 4.5 wurde gezeigt, daB diese Determinante gleich der Summe der Produkte 
von allen k —1 Baumadmittanzen aller méglichen Baume des Netzes ist. Wir schlieBen nun 


wieder aus den Eigenwerten von Yj auf diejenigen von Y+. Falls alle Baumadmittanzen 
positive Realteile besitzen, ist zunachst 


det Yju= det'(Re Yt) = det Gt = Ds 
My ss Hp 
= Kombination 
aller Baumzweige 


oe Reyes, 
da dann wenigstens ein Glied der Summe nicht verschwindet und auBerdem positiv ist, wobei 
wir natiirlich noch stillschweigend vorsetzen miissen, daB in den unabhangigen Zweigen keine 
Admittanzen mit negativen Realteilen vorkommen. Es sind aber auch alle Hauptabschnitts- 
unterdeterminanten positiv, fiir welche die gleiche Regel fiir die Summendarstellung gilt. 
Der gestrichenen Knotenpunktszeile und -spalte entspricht namlich nicht nur ein gestrichener 


Bx, DH IS pai Oa Onc fo eee Wp een) (67) 
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Knotenpunkt, sondern auch ein gestrichener Baumzweig, und zwar derjenige (es gibt nur 
einen), der im Weg zum Bezugsknotenpunkt enthalten ist. Es ergibt sich auch hier wenigstens 
ein positives Glied in der Summe, das einem Unterbaum entspricht. Damit sind die Eigen- 
werte von Yj sdmtlich negativ und nach dem bereits erwahnten Satz von A. Hirscu ebenso 
die Realteile der Eigenwerte von Y*. Keiner der Eigenwerte verschwindet und somit ist 
det Yt4o. 
Damit gilt der folgende Satz: 


Dafiir, daB die Knotenpunktsgleichungen 
YoU Ker (22) 


eindeutig auflésbar sind, ist hinreichend, daB das Netz wenigstens einen Baum enthalt, 
dessen Zweigadmittanzen ausschlieBlich positive Realteile besitzen (und die iibrigen 
| keine negativen Realteile). 


6. Anwendung der Existenzsatze auf Energieiibertragungsnetze 

Offenbar kommt es in den Netzberechnungen nur darauf an, daB in der Maschenmethode 
Kreisstréme vermieden werden und in der Knotenpunktsmethode keine Spannungen bestehen 
bleiben kénnen, die nicht durch endliche Impedanzen kurzgeschlossen werden. Da die Be- 
dingungen nur hinreichend sind, und man auch leicht Beispiele mit verschwindenden ohm- 
schen Verlusten angeben kann, in welchen die Gleichungsdeterminante (fiir eine festgehaltene 
Frequenz) nicht verschwindet, ware es sinnvoll, nach weiter reichenden zumindest hinreichen- 
den Bedingungen zu suchen. Doch zeigt sich in der Praxis, daB man mit den angegebenen Be- 
dingungen in den meisten Fallen auskommt. Man muB nur die in dem gegebenen Netz auch 
wirklich auftretenden Verluste beriicksichtigen. Insbesondere sind die Verluste in den induk- 
tiven Zweigen normalerweise so groB, daB wenigstens die GréBenordnung bekannt ist, was bei 
Kapazitaten oft nicht der Fall ist. Man muB8 nun also im Fall der Maschenmethode nur das 
entsprechende System unabhangiger Zweige bzw. im Fall der Knotenpunktsmethode einen 
vollstandigen Baum finden, in welchem die Serien- bzw. Parallelverlustwiderstande auftreten. 


Kraftwerk1 Transf.1 Leitung1 Leitung 2 Leitung 3 Transf.2 Kraftwerke 


Transformator 3 


Kraftwerk 3 


Bild 7. Energiegraph eines einfachen Drehstromverbundnetzes. 


Betrachten wir z. B. ein Drehstromverbundnetz, dargestellt durch seinen Energie-Graphen? 
in Bild 7. Fiihrt man die Nennspannungen als Bezugsspannungen ein und nimmt man weiter 
an, daB die Transformatoren mit ihren Nenn-Ubersetzungsverhaltnissen arbeiten, so laBt sich 
dar Netz iibertragerfrei darstellen? (entsprechend der tiblichen Darstellung auf einem Netz- 


1 Energie-Graph’ = Graph der Energiefliisse, mit Vierpolen als Zweige des Graphen, nach einem Vor- 
schlag von Herrn Dipl-Math. G. NEEs. 

2 Enthalt der Energie-Graph Maschen, so kann es vorkommen, da8 im Fall, da8 das Produkt der Uber- 
setzungsverhaltnisse in einer Masche nicht 1 ist, ein Restiibertrager bleibt. In diesem Fall kannn der Ubertrager 
zusammen mit der nachzubildenden Streuimpedanz durch einen aquivalenten RLC-Vierpol nachgebildet 
werden. 
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modell). Dann enthalt das Netz eine durchgeschaltete Nullschiene (welche der Erdriick- 
leitung entspricht). 

Mit diesem Knotenpunkt als Mittelpunkt kann man unter Verwendung samtlicher Quer- 
schnitte einen sternférmigen Baum auswahlen (Bild 8). Das System unabhangiger Zweige 
enthalt von den treibenden Spannungen abgesehen, verlustbehaftete Induktivitaten. Die 
Realteile der Impedanzen in den unabhangigen Zweigen sind sdmtlich positiv (dartiber hinaus 
auch die Realteile einiger Baumzweig-Impedanzen, welche Wirkbelastungen entsprechen). 
Nach den Sdtzen in Kap. 5.1 existiert somit eine eindeutige Lésung der Maschengleichungen. 


Sele ae goile 
Xe XC, X, Xe, 


124 


Bild 8. Ein System unabhangiger Zweige fiir das Netz nach Bild 7, das nur verlustbehaftete 
Iuduktivitaten enthalt, so daB die Maschengleichungen eindeutig auflosbar sind. 


JK, ® My fu, JX Ff, 


4a, Ki, 


Bild 9. Ein vollstandiger Baum aus dem Netz nach Bild 7 mit ausschlieBlich verlustbehafteten Admittanzen, 
so daB die Knotenpunktsgleichungen eindeutig auflésbar sind. 


Auch fiir die Knotenpunktsgleichungen l48t sich ein vollstandiger Baum so bestimmen, 


daB in jedem Zweig positive Realteile der Admittanzen vorhanden sind, z. B. dadurch, daB 
man als Baum ein System von Langszweigen verwendet, das also dem entsprechenden Energie- 
graph entspricht, und bei welchem im Falle einer Vermaschung noch Zweige herausgenom- 
men werden miissen und ein verlustbehafteter Zweig, z. B. ein solcher, der eine Energiequelle 
enthalt, hinzuzufiigen ist, um die Verbindung mit der Nullschiene herzustellen (Bild 9)?. 


 Fiir die Knotenpunktsmethode benétigt man die Baumzweig-Konzeption eigentlich nicht, sondern nur 
um die Eindeutigkeit der Lésungen nachzuweisen. 
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Ein solches System von Baumzweigen enthalt in jedem Zweig Verluste, d.h. die Baum- 
zweigadmittanz besitzt einen positiven Realteil und nach dem Satz in Kap. 5.2 existiert eine 
eindeutige Lésung fiir die Knotenpunktsgleichungen. 


Zusammenfassung 


Ausgehend von der Maschen- und der Knotenpunktsmethode und ihren Darstellungen in 
Matrizenform wurden Wege fiir die Aufstellung der die Topologie des Netzes charakteri- 
sierenden Inzidenzmatrizen gegeben. Wahrend die Aufstellung der Knotenpunkt-Zweig- 
Inzidenzmatrix keine besonderen Schwierigkeiten macht und sich sofort aus dem aufgezeich- 
neten Netz oder der Liste aller Zweige ergibt, muB der Aufstellung der Maschen-Zweig-Inzi- 
denzmatrix die Festlegung eines vollstandigen Baumes vorausgehen. Dies ist aus dem auf- 
gezeichneten Netz relativ leicht méglich. Fiir die Rechnung auf einem digitalen Rechenauto- 
maten benétigt man aber einen Weg, dies aus der Liste aller Zweige oder aus der Knotenpunkt- 
Zweig-Inzidenzmatrix heraus durchzufiihren. Aus der Definition des vollstandigen Baumes 
1aBt sich ein entsprechendes logisches Programm leicht herleiten. Der nachste Schritt ist die 
Ermittlung der Maschen-Zweig-Indizenzmatrix. Sie kann nach Einfiihrung einer Normal- 
numerierung der Knotenpunkte und Zweige aus der Knotenpunkt-Zweig-Inzidenzmatrix 
in einem besonders einfachen EliminationsprozeB (Aufrollen eines Gleichungssystems mit 
den Koeffizienten 0, +1, —1 in einer Dreiecksmatrix) errechnet werden. Bei diesen Be- 
trachtungen ergeben sich nebenbei einige wichtige Satze der allgemeinen Netzwerktheorie, 
die schlieBlich auch in Beweisen fiir die Existenz eindeutiger Lésungen bei Erfiilltsein gewisser 
Passivitatsbedingungen Verwendung finden. 
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A. Ejinleitung 
1. Allgemeines 

Infolge seiner vielfaltigen Anwendungen in der Technik beanspruchte das magnetische 
Feld von jeher einen wichtigen Platz in der Elektrodynamik. Die genaue Kenntnis der differen- 
tiellen Struktur magnetischer Felder war daher bald ein dringendes Problem. Mathematische 
Hilfsmittel, wie verschiedene krummlinige orthogonale Koordinaten, die GREENschen Satze- 
allgemeine Funktionentheorie und konforme Abbildung sind nur fiir Felder konstanter Permea, 
bilitat brauchbar und versagen in wichtigen Fallen. Besonders fiir die Berechnung stark 
gesdttigter magnetischer Kreise wurden daher schon friith weitere Méglichkeiten gesucht, 
und diese teilweise in verschiedenen Analogiemethoden und graphischen Verfahren gefunden. 


2. Die Differentialgleichungen fiir das magnetische Potential 


a) Das Skalarpotential 
Zu untersuchen ist die differentielle Struktur des magnetischen Feldes in stromfreien 
' Gebieten. Da die elektrische Stromdichte G verschwindet, so mu8 nach der 1. MaxweELLschen 
Gleichung die magnetische Feldstarke H wirbelfrei sein,<d: hi. 


rot H =o. (A1) 
Fiir diese Gleichung lautet die allgemeine Lésung: 
je grad V 


worin V,, die skalare magnetische Potential-Ortsfunktion bedeutet. Aus dem Zusammenhang 


(A2) 


m? 


zwischen B und H in seiner speziellen Form: 


Bit = eu H (A3) 
erhalten wir unter Beriicksichtigung des Kontinuitatsgesetzes fiir die magnetische Induktion, 
divB=o, (A4) 


die Differentialgleichung fiir V,,: 
gradu grad V,, + wdAV,=0. (A5) 


Speziell fiir Felder mit einer skalaren konstanten Permeabilitat w gilt, daB grad w= 0 mit 
pu % 0, so daB Gl. (A5) in die Laplacesche Gleichung iibergeht: 


AV,, =0. (A6) 
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Damit ist die Analogie mit stationaren elektrischen Stroémungsfeldern gegeben. Fiir diese gilt 
bekanntlich: 


0, 
5. ees (A7) 
0 ; 


b) Das Vektorpotential 


Im wirbelfreien Feld ist der Feldvektor H als Gradient einer skalaren Ortsfunktion Von 
darstellbar. Mit dem magnetischen Skalarpotential kénnen damit nur stromfreie Gebiete 


beschrieben werden. 
Zur Darstellung stationaérer Magnetfelder beliebiger Art folgt als Lésung der universell 


giiltigen Kontinuitatsdifferentialgleichung (A4) fiir die magnetische Induktion B der Vektor 
Vm, aus dem B nach 
(A8) 


abgeleitet wird. 
Zusammen mit der Gl. (A3) und 


> 


Tot Hie : (AQ) 


ergibt sich die Differentialgleichung fiir das magnetische Vektorpotential unter Annahme einer 
konstanten Permeabilitat : 
rot rot V, = Ge 
oder: 
AV, — grad div-V,, 1G =o. (A 10) 


Ist wim betrachteten Gebiet eine skalare Veranderliche, so ergibt sich: 
AV, +p rot as x grad 7 — grad div VG aon (A11) 


Diese Gleichung soll zur Berechnung von Vn dienen. G kann mittels Gl. (Ag) aus H gefunden 
werden. Ist G gegeben, so mu8 fiir die Bestimmung von H die Gleichung (Ag) integriert 
werden. Ein allgemeiner Ansatz ist Gleichung (A8) in Verbindung mit Gl. (A11), welche 
V,, mit G verkniipft. Damit das Integral Gl. (A8) eindeutig wird, soll das Vektorpotential 
selber quellenfrei sein: 


div. V,, =0., (A12) 
Damit folgt nun fiir konstante Permeabilitat die Gleichung: 
AV, =—pG, (A13) 
und fiir variable Permeabilitat : 
AVR rot Vy, x gradZ|+ wG = 0, (A14) 


wobei die Gl. (A13) genau der Potssonschen Gleichung im elektrostatischen Fall entspricht. 


B. Analogienetzwerkmethode 


Mit dem Analogienetzwerk kénnen inhomogene mtagnetische Felder nachgebildet werden, 
in welchen die Permeabilitat eine stetige oder unstetige Funktion des Ortes ist. Dies ist méglich 
sowohl fiir das magnetische Vektorpotential wie auch fiir das Skalarpotential, doch soll sich 
die vorliegende Untersuchung auf den letzteren Fall beschranken. 
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1. Zweidimensionaler Fall 
a) Konstante (feldunabhdngige) Permeabilitat 


Die Feldstruktur gehorcht der Gl. (A6). Zuerst sei der Fall einer als Funktion des Ortes 
stetigen Permeabilitat betrachtet. 

Nach Bild 1 unterteilen wir das untersuchte Gebiet in Aquidistante Maschenpunkte mit der 
Maschenweite d. Der LapLacEsche Operator A wird in Differenzenform angeschrieben; der 
Vergleich der ,,magnetischen Ebene“ mit der analogen ,,elektrischen Ebene“ ergibt: 


1=1 


4 4 
== Ve > 3 Kan = 0; a ys .— = 
2 ( 4 .) Ry ios (V; Vo) a (B1) 


Die Konstante ky verkniipft das magnetische Potential V,, mit dem elektrischen Potential V, 
d. h. den physikalischen Raum mit dem Modell. ky hat die Dimension des elektrischen Wider- 
standes. 


Bild 1. Vergleich der ,,magnetischen Ebene“ (links) mit der analogen ,,elektrischen Ebene“ 
(rechts) bei stetiger und konstanter Permeabilitat py. 


BS . 
mw Me, dh | eq 
WN 


\\ as I 


Bild 2. Analogie entsprechend Bild 1 bei unstetiger Permeabilitat yu. 


Die Komponenten Z,, H, der magnetischen Feldstarke berechnen sich zu: 

Vin Va, = 1 
d ho a 

Vian eae 
d fig a 


(Vo aiTG V2) , 


H, = — grad, V,, = 


x 


(B2) 


H, = — grad, V,, = (Vo — V4) . 
Bei Unstetigkeit von mw sind an der Sprungflache zwischen dem Medium I und dem Medium II 
die Normalkomponenten B,, der magnetischen Induktion B und die Tangentialkomponenten H, 

- der magnetischen Feldstarke H stetig: 


BeBe era, =o (B3) 


Ein aus identischen Widerstanden zusammengesetztes Netzwerk mu8 am Rande mit dem 
doppelten Netzwiderstand abgeschlossen sein, damit keine StoBstelle auftritt [1]. 
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Betrachten wir nun den speziellen Fall eines Uberganges von einem Medium I mit der 
relativen Permeabilitat w, auf ein Medium II mit w, = 1, so ergibt sich gemaB Bild 2 der 
Widerstand R, zu: 


: (B4) 
4 Re Ree 


Durch Betrachtung der Differenzenflachenelemente ergibt sich, daB die Bedingung H, = H. 


trr 


a priori erfiillt ist, denn der Grenzstrom I (I; bzw. I,) ist derselbe fiir beide angrenzenden 

aR, aR “ORI Flachenelemente, und nach hyd H, 
= RJ kann an der Grenzflache nur 
@% eine und dieselbe y-Komponente 


von H herrschen.. 


Die Bedingung B,, = B,,, ergibt: 


Je = bo Hy dels: = Mo Has ’ 


Ro I, oF Ro Ip oc ae 
2ky Pola i Saag Srl Orang , 


aR} M, Ryo = Ry, 


oder, falls der Widerstand Rj als Be- 
zugswert genommen wird: 


Ry Sth ae (B5) 


Bild 3 zeigt die Modellanordnung fiir 

ein Winkelprofil, einen in der Praxis 

‘ afiny SPSS SSS sat tee ah y ah des Transformatorenbaues wichtigen 

Pua ce noranins Oe ee oor Wieolp aite ase Tees Fall [2]. Der auBere Teil (schraffiert) 

ist nicht nachzubilden, da dort der 

Streuflu8 vernachlassigbar ist. Der Widerstand R, an den Grenzflachen betragt R, = 2 Rj 

(1 + u,)7, wahrend kj O = — V den Faconmnenhane der anzulegenden elektrischen opianeee 
V mit der Durchflutung @ ergibt. 


b) Verdnderliche Permeabilitat 
Fiir diesen Fall lautet die durch my dividierte Gl. (A5): 
grad u, grad V,, + uw, AV,, = 0. 
Zuerst behandeln wir wieder den Fall der 6rtlichen Stetigkeit von w. Nach Bild 4 folgt: 


4 
grad u, grad-‘V,, + 4& on 7 = (Vg eV pe) (B6) 
oder mit den Komponenten von grad mw, und grad V,,: 
Mr, br 
grad, mu, = —, 
SP Ld Id 
grad, u, = ee 
erad Votes Maes i ; 
2d 
Ving ay Ving Y 
grad, V,, = Soe 


grad yu, grad V,, = Bose { pce a) re uh = Ving 
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ergibt sich nach kleiner Zwischenrechnung: 


(Fy — Vin) (hy — =) (Vong — Vin) (Hy + HE) Vg, — Pru) (1, — MH 


+ (Vn, = Vin) (1p + MMe) = 0. (B7) 


Darin bedeutet yw, die Permeabilitat 
im Potentialpunkt V,,,. Es ist nun: 


Mr, — Me 
Ly, — '3 ; 1 = U,,, ; 


Lr, — by 
bret. Ly = = bry » 


Pra Pry __ 
1 2 ae ga 


ue we Lr, — Erg ei, 
f 4 Koa ss Bild 4. Analogie entsprechend Bild 1 bei drtlich variabler Permeabilitat pu. 


wobei y,,, die auf die Maschenweite bezogene relative Permeabilitat zwischen den Potential- 
punkten V,, und V,,, bedeutet. Gl. (B7) schreibt sich dann: 


4 
S bry, (Vo Vind =n (B8) 


Fiir das elektrische Analogon folgt: 
4 


F=V,—V,) =0. (Bo) 


Multiplikation von Gl. (B8) mit dem Produkt k, k>+, in welchem k, eine Konstante von der 
Dimension des elektrischen Widerstandes ist, und Vergleich mit Gl. (Bg) fiihrt auf die Be- 
ziehung : 

he 


tere i 


R.= 


4 


(B10) 


Mit der Annahme der Konstanten k, ist damit die jeweilige GréBe der Netzwerkwiderstande 
festgelegt. Die Komponenten der magnetischen Feldstaérke gehorchen den Gln. (B2). 

Fir eine Unstetigkeit von m stellen sich bei einem Ubergang von Medium I\ (mu) u,) auf 
Medium II(u,) die Bedingungen der Gln. (B3): 


regal ee wea, DILOLE Eli. 


ty eae 
oe = Bees Mo M7, Ay, = lo Ho. , 
Ly, Ry = Ry ’ 
d.h. es wird gemaB Annahme R’ = hy. 


Als Anwendung betrachte man Bild 3; darin muB lediglich fiir die Nachbildung des py u,- 
Raums der konstante Widerstand R, durch den von der relativen Permeabilitat My, ab- 


hangigen Widerstand Rk, nach Gl. (B10) ersetzt werden. 


2. Eindimensionaler Fall 


Der eindimensionale Fall kann mittels einer einfachen Widerstandskette gelést werden. 
Die Differentialgleichung (A5) reduziert sich auf: 


2 
en eee ae (<" == Qi, (B11) 
dx? pe \dx/} \ dx 


31% 
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Der allgemeinste Fall, der mit einer ebenen Kette gelést werden kann, ist die Hintereinander- 
schaltung von n Medien der relativen (feldabhangigen) Permeabilitaten w,. Da sich diese 
jedoch in y- und z-Richtung ins Unendliche erstrecken miissen, (d. h. die Grenzflachen miissen 
geometrische Ebenen sein), ergibt sich tiber x keine stetige Variation von yw, sondern nur Strek- 
ken konstanter Permeabilitat und Spriinge an den Grenzflachen. Damit vereinfacht sich 
Gl. (Bi1) auf: 

evens O. (B12) 


An einer Unstetigkeitsflache gilt gemaB Bild 5: 
R, = Bo fe. 
"es pee 


Fiir den Spezialfall der Hintereinanderschaltung von zwei Medien der Langen uN = md und 
1, = nd und der relativen Permeabilitaten uw, und y,, ergibt sich Bild 6, das vereinfacht nach 


Bild 7 dargestellt wird. Dabei gilt: 


Vi—V_ _ Vi—Va 
mR, + nRy> Ry Ry” 
R, + Rg 


an (Ve) 2a 
mR, +n, M Uy + ™ Me 


Die reduzierte Spannung (V; — V,) gestattet das Ersetzen der Kette mit m +  Widerstanden 
durch eine Kette mit 2 Widerstanden. 


(Vo Valea Va) 


SS Hy 


Bild 6. Nachbildung der Hintereinanderschaltung von Medium py, Bild 7. Reduktion der Widerstandskette 
der Lange /, mit Medium yw, der Lange /,. aus Bild 6 auf 2 Widerstande. 


Der allgemeine Zusammenhang zwischen dem magnetischen Flu8 ® und dem elektrischen 
Strom J im Modell ergibt sich folgendermaB8en: Fiir den FluB gilt: 


O=f[B-dA=B.-A. (B13) 
Das Flachenelement A hat die Breite d und die Tiefe 1. Es wird damit: 
Dio 
@ = fot, HA = Mo H, aa 
0 
a Ray | 
Ho ho ? | (B 14) 
k 
O = UT, . 
Ho 5 J 


Es sei hier noch bemerkt, daB man alle hier betrachteten Falle statt mit variablen Wider- 
standen auch mit gebietweise konstanten Widerstanden nachbilden kann. Im allgemeinsten 
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Fall der 6rtlich variablen Permeabilitat gemaB Bild 4 ergibt sich der Einspeisestrom fiir den 
Maschenpunkt V, = ky V,,, als Funktion der 4 umliegenden Potentialwerte V,...V, durch 
entsprechende Zusammenfassung der Gl. (B7). 


3. Praktische Anwendung der Netzwerkmethode 
a) Allgemeines 


Die Maschenweite d des Netzwerkes ist im Prinzip ein frei wahlbarer Parameter. Nach oben 
ist sie jedoch durch den maximal zulassigen Fehler, nach unten durch die Gestehungskosten 
des Modells und den Zeitaufwand zur Messung begrenzt. (Ein 
rechteckiges Netzwerk mit m-n Punkten hat exakt 2mn-(m + n) 
Verbindungswiderstande.) 

In jedem Fall ist fiir die Dimensionierung von d eine opti- 
male Lésung anzustreben im Hinblick auf Aufwand und er- 
forderliche Genauigkeit, d.h. man wahlt fiir einen maximalen 
Fehler die gréBtmégliche Maschenweite. 

Die Probleme der Widerstandstoleranzen, der Festlegung der 
Maschenweite, der Randwerte, der Widerstandssdimensionierung, ; : 
der Randabschliisse bei krummlinigen Randern usw. sollen hier Fa on ae Tee 
iibergangen werden, da sie den Rahmen der vorliegenden Arbeit aire 
tibersteigen [3]. 
_ Das Vorgehen zur Lésung eines Feldproblems ist nun das folgende: Die Durch finer 
O=V,,,—V,, wird als Spannung V = —k,@ an das Netzwerk angelegt, wobei alle 
variablen Pease einen gewissen, gleich groBen Wert haben, der von Fall zu Fall zwischen 


‘ i gewahlt wird (kg = Rj ist der Widerstand zur Nachbildung des Mediums mit 
VMAX ymin 
Hb, = 1). Dann wird das entstehende Potentialfeld ausgemessen und fiir die einzelnen Maschen- 


mittelpunkte P die dort herrschende Feldstarke H, wie folgt berechnet (Bild 8): 


wa YC, 


Y /P 
5 2 2 (B 1 5) 
H $e Vina ee Ving Ving a5 Ving— Vin Fe Ving 
f= ESS ie eee 
da 
| tC) im (Gries Ving 
|\ 8% /p d | 
il 1 | 
(Fe) | zy Ving ap Ving alee a Vina + Ving | 
ay }p| + -|- ees: | 


Aus der Magnetisierungskurve uw, = f/(H) erhalt man yp, Be P. Mit diesem Wert y,,, wird der 
Netzwiderstand R, nach der abgeleiteten Formel Ry = a korrigiert. Dies wird solange 


wiederholt, bis die Konvergenzgeschwindigkeit Null wird. Zwischen dem gesamten, ins Netz- 
werk flieBenden Strom J und dem Flu8 @ im physikalischen Raum besteht dann Proportionali- 
tat, so daB man direkt einen Punkt der Kurve ® = @(@) d. h. der magnetischen Charakteristik 
fiir den betrachteten Kreis, erhalt. 


b) Praktisches Beispiel 
Betrachtet sei ein magnetischer Kreis, der sich aus je einer Strecke Ferromagnetikum (I) 


eset B 
und Luft (II) gleicher Lange zusammensetzt. Bild 9 gibt die fiktiven Kurven y,(H) und @ (H) 
0 


des Mediums I, Bild 10 die resultierende magnetische Charakteristik. Aus der durch die Permea- 
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bilitat yz) des leeren Raumes dividierten Induktion B kann durch eine einfache MaBstabs- 
anderung der magnetische FluB @ erhalten werden: 


B 
©=BA=(*) yd. 


Die magnetische Charakteristik hat das erwartete Aussehen: Fiir kleine Erregungen @ ist 
sie linear; der Einflu8 der Sattigung macht sich erst bei gr6Beren @ bemerkbar. ~ 


Oo 
0 100 200 300 A 400 


Bild 10 Magnetische Charakteristik der Serieschaltung 
Ferromagnetikum—Luft. 


Bild 9 (links). Kurven y, (H) und B(#) eines fiktiven Ferromagnetikums. 
Up (H) ist angenommen und B(H) daraus berechnet. Der ,,spezielle“ 
Verlauf der Permeabilitatskurve y(H) ist festgelegt im Hinblick auf 
Konvergenzuntersuchungen mit verschiedenen Permeabilitatskurven, 


nae usw. konvergenzverschlechternd 


wirken. 


welche durch Spitzen, groBe 


c) Konvergenz 


Theoretisch ist tiber die Frage der Konvergenz fiir die Analogie-Netzwerke und die ihnen 
zugrundeliegende Relaxationsrechnung wenig bekannt. In unserem Fall wird eine analytische 
Untersuchung noch dadurch erschwert, daB von einem Konvergenzschritt zum anderen iiber 
eine komplizierte Kurve wu, = u,(H) vorgegangen wird. 


GemaB FISHER [4] kann die prinzipielle Konvergenz bei den betrachteten magnetischen 
Feldproblemen als gesichert angenommen werden. So bleibt nur noch die Giite der Konvergenz 
zu untersuchen. Dies geschieht am besten an einem einfachen Fall wie z.B. dem weiter oben 
betrachteten. 


Wahlen wir etwa den Konvergenz-Punkt 0 = 100 A (Bild 10) und bestimmen die not- 
wendige Anzahl Iterationsschritte bei Annahme verschiedener nullter Naherungen: 


Nullte Naherung: H = 0 (1p sa) Schrittes?7 
kA 
Nullte Naherung: H=13—~- pm, = 2 Schritte: 8 


kA 
m 


Nullte Naherung: H =1 20 Schritte: 5 


Die Konvergenz ist bei der dritten Annahme am besten. Dies ist fiir den betrachteten 


5 f : ; kA 
Punkt nicht weiter verwunderlich, da die Annahme H = 1 5 dem Konvergenzwert H = 591 ic 
am nachsten liegt. ne 


In der Tat ergibt sich allgemein, daB die Giite der Konvergenz allein von der Toleranz in 


der nullten Naherung abhangig ist, wahrend der Ort des Konvergenzpunktes auf der Magneti- 
sierungskurve kaum einen EinfluB hat. . 


Bei eingehenden Untersuchungen mit verschiedenen Permeabilitatskurven u,({H) konnte 


festgestellt werden, daB bei negativer reversibler Permeabilitat im betreffenden Intervall von H 
keine Konvergenz vorhanden war. Solche Falle, bei denen die Magnetisierungskurve ein Gebiet 
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Tae : ; ; 
mit mir <o, d.h. einen fallenden Ast aufweist, miissen deshalb bei den Anwendungen aus- 


geschlossen werden. Dies ist weiter keine Einschrankung, da solche Magnetisierungskurven in 
der Praxis gar nicht vorkommen. 


C. Graphisches Verfahren 


Graphische Verfahren zur Bestimmung von magnetischen Feldern sind seit langem bekannt. 
Fiir die Praxis sind sie deshalb attraktiv, weil sie bei geringem Aufwand an Material und Zeit 
meist brauchbare Lésungen ergeben. Die Entwicklung der graphischen Verfahren fiir Magnet- 
felder nahm ihren Anfang durch die Arbeiten von LEHMANN [5] zu Beginn unseres Jahrhun- 
derts. In neuester Zeit beschrieb BENEDIKT [6] ein nomographisches Verfahren, das iiber die 
Bestimmung der magnetischen Potentialverteilung bei gegebener geometrischer Kon- 
figuration und gegebener Induktionsverteilung hinausgeht. Es erlaubt ndmlich, zwei 
beliebige dieser drei GréBen anzunehmen und daraus die dritte zu bestimmen. 

Der vorliegende Beitrag fiihrt in Erginzung und Erweiterung der Lehmannschen Verfahren 
auf die graphische Bestimmung des magnetischen Feldes im Innern eines Ferromagneti- 
kums bei beliebiger Geometrie. 


1. Grundlagen 


Die Methode beruht auf der Zerlegung des zu untersuchenden Gebietes in einzelne FluB- 
rohren. Zwischen der mittleren Lange Al eines Abschnittes der Réhre und ihrer mittleren 
Breite 4d bzw. ihrem mittlern Querschnitt 4A besteht dann ein aus der Magnetisierungs- 
kurve 6(H) folgender funktionaler Zusammenhang, anhand dessen das Feldbild bestimmt 
werden kann. 


Vin= konst 


Bild 11. Kraftréhrenabschnitt mit dem TeilfluB 49 Bild 12. Zur Anwendung der graphischen Methode auf ~ 


und der Durchflutung 4 0 = VinsV m, einen radial durchflossenen Kreisring. 


Betrachtet man einen Kraftréhrenabschnitt nach Bild 11, so folgt aus den Gln. (A2), 
(A3), (B13) und der Beziehung 


@=-/Hd=—Yv, (C1) 
folgendes: Der Teilflu8 A@ ist innerhalb der Kraftréhre konstant: 
NO— BNA —— Bb iAbid cs, (C2) 


worin d die Abmessung des Ferromagnetikums normal zur Zeichenebene bedeutet; ebenso ist 
zwischen zwei Potentialflachen die Teildurchflutung 40 konstant: 


NO AV ss Holl 21K. (C3) 
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Die GréBen K, (Dimension: Ampére) und K, (Dimension: Voltsekunden) beschreiben den FluB- 
rdhrenabschnitt und zusammen’ mit der Magnetisierungskurve B(H) den Feldverlauf. Aus den 
Gln. (C3) und (C2) folgt nun unter Beriicksichtigung der Gl. (A3): 
kK. K. K. 
Ab - = See A Ad, | ¢ 
H B K, 
; (C4) 


Ky, 
AL K, GieeANUm 
Damit ist der Zusammenhang zwischen Al und Ab gegeben. Darin ist speziell: 
K, 
= u(B) = 1 (32). (C5) 
Ka, Ket Parameter. 


2. Anwendung auf Kreisring 


Betrachtet wird ein Sektor nach Bild 12 mit einem Zentriwinkel a = 10° + 7/18, bestehend 
aus Si-Fe-Blech mit einer Dicke d = 0,5 mm und einer Verlustleistung von 1,7 W/kg, dessen 


B B 
26 
V. 
¥y 2,0 
2A Vs 
a4 ™ 
20 4,5 
18 
46 
40 3 
14 
42 
40 05 
08 
H 
0894 1 40 700 KA 1000 Z 
™ 0 kA 0,4 
Bild 13. Gleichstrommagnetisierungskurve fiir Si-Fe-Blech mit 1,7 W/kg : Ht e ie Th 
Verlustleistung. 0 4 8 12 16 
EE ee ee | ee 
Bild 14 (rechts), Wie Bild 13, jedoch mit linearem MaBstab fiir H. 0 ia s ie) yy 


Gleichstrom-Magnetisierungskurve in Bild 13 und 14 gegeben ist. Aus dem angenommenen 
B, = 2,1 Vsm™, entsprechend H, = 10,6- 104 Am™ (Bild 13, 14), ergibt sich mit Bild 12: 


K, = A®, = B,- AA, = B,},d = B, 37, d = 2,747- 10+ Vs; 


b= eile 2,618 cm, 
17 
b= mate 2,965 cm 
l Yy— ¥ 
Annahme: Al, = — = “~—* =0,4mm, 
5° 5° 


it BAL 42 ae 


Mit K, und K, kann der Zusammenhang Al = Al(b) berechnet werden wie folgt: 
1. Annahme von 0, 


B,b 
2. Berechnung von B = ; ,, . 
3. Ablesen von H(B) aus Bild 13, 14 
PRT ge ES Bay ocak er 
4: BES Bae AR get ae et pee oa 
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Die berechneten Werte sind in der Zahlentafel 1 zusammengestellt und in Bild 15 auf- 
getragen. 


Zahlentafel 1 


b cm 2,60 | 2,65 2,70 | 2,75 2,80 | 2,85 | 2,90 | 2,95 | 3,00 
| 
| Vs 

B a 2,112] 2,073] 2,035 | 1,995] 1,960] 1,925) 1,893/ 1,862) 1,830 
| | 
logs | 

H | 108 11,6 85 00a 052 4,6 3,4 2,65 | 2,23 ThOOm els OO 
| | | | | | 
[es | 

1 mm 0,365 | 0,493 | 0,6845| 0,944/ 1,246)| 1,6 1,90 | 2,23 | 2,65 


Die Auswertung geschieht gemaB der Beziehung AJ. H = Al, H, = K, durch Aufteilung 
von / in einzelne Kraftréhrenabschnitte mit in Richtung 7 konstanter magnetischer Potential- 
differenz in Bild 16. Dabei wird jeweils die mittlere Breite b jedes folgenden Elementes ab- 


b 
3,0 
cm 
2,5 
Al 
30 
se 20 
20 
2,0 45 
fi) 
1,0 
1,0 
0,5 
05 
b 0 ie 
05 37 28 29 om 30 15 16 cm 17 
Bild 15. Nichtlinearer Zusammenhang 4/(b) fiir das kreisring- Bild 16. Kartesische Darstellung b = b(7) des Kreisring- 
formige Si-Fe-Blech nach Bild 12 und 13 mit K,; = 42,4A sektors von Bild iz zur Aufsummierung der K,-Abschnitte. 


und K, = 2,747:+ 10-4 Vs. 


geschatzt und die dazugehérige Lange A/ des Abschnittes Bild 15 entnommen und abgetragen. 
Die gesamte Durchflutung ergibt sich damit zu: 
CO 21-4800. A. 
’Wird die Durchflutung aus dem linearen Mittelwert der Induktion berechnet, ergibt sich 
(Seba 12);mit = by = (2:1 Venr?; 
Bae easier es 2 ( “) = 1,976 Vsm=, 
2 


2 2 


Pie—s3,0 sO; ATi (sils-bildit4) 
der Wert 0 = H,, (re — 7) = 780 A. Die auf den Wert 890 A bezogene Abweichung betragt 
ae 10) 
Tal ales der Stmpsonschen Regel erhalt man: 
0 = ees! (72 — 1) 
Sey OROIG AS = gee tee (17 — 15) 10-2 = 932 A, 
d. h. eine bezogene Abweichung von 5%. 


é Archiv fiir 
452 P.A. Tscuopp u. A. H.. Fret: Analogienetzwerkmethode und graphisches Verfahren ... ic ree ema 
en aie rtaieinanbtiat Some eee Serbs Jian a de STC NS eR 


SchlieBlich folgt bei Anwendung der doppelten Srmpsonschen Regel: 


By = 21 VYsm* HN "FOrO. 1O-sAti es 
Bae 208 H; = 6,9 +104 
Bz = 1,97 H; = 3,05 - 10% 
Be= a Ok Hy = '2,45*> 10" 
Boi "1,055 Hyx=- 1,806 10. ' 
0 = 0,4 0,2 74ers (24) + Bata” (54) = 879A, 


mithin eine bezogene Abweichung von — 1,3%. 


3. Anwendung auf Winkelprofil 


Untersucht wird das Winkelprofil Bild17. Fiir die Parameter K, und K, ergibt sich 
daraus: 
Kio Be AA =e 107d. = 20,25 alOr avs 
Ka Hye Ale AAG) ae ae 
Der Zusammenhang Al (Ab) wird wie bei der Zahlentafel 1 berechnet ; er ist in Zahlentafel 2 zu- 
sammengefaBt und in Bild 18 aufgetragen. 


Al 


4b 
0 Sea 20 30 «mm 40 
Bild 18. Nichtlinearer Zusammenhang 4/ = Al(Ab) 


fiir das Winkelprofil Bild 17 bei K,=2,1A und 
Ke = 16,2515 Lome Ss 


Bild 17. Rechtwinkliges Winkelprofilaus Si-Fe-Blech nach Bild 13. Blechdicke d = 0,5 mm. Beziiglich A h— Ab 
sei das Profil symmetrisch, im Schnitt Ag — Ag das Feld homogen. 


Zahlentafel 2 


Ab mm 6 8 10 12,5 15 20 30 40 50 
Vs 

B mm? 2,085 1,563 1,250 | 1,000] 0,834] 0,625] 0,417 |0,3125| 0,250 
A 

a m 9,5+10* | 4,63-10%| 500 | 210 | 144 | 94 |.59,5 | 43 34 

Al mm 0,022 0,454 4,20 10 | 14,58 | 22,35 | 35,30 | 48,80 | 61,75 


Das Vorgehen zur Bestimmung der Feldlinien und. der Potentiallinien ist nun iterativ. 
Die einzelnen Iterationsschritte sind die folgenden: 


Nullter Schritt: Aufteilen des Profils in =; 


= 8 FluBréhren durch Annahme von 
plausiblen Feldlinien (Bild 19). ‘ 
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Erster Schritt: Unterteilen jeder FluBrohre anhand des Zusammenhangs Al = Al(Ab) 
aus Bild 18 in Abschnitte konstanter magnetischer Potentialdifferenz, beginnend im homogenen 
Bereich. Es resultieren Treppenfunktionen fiir V,, (Bild 19) ; 
offenbar weicht die nullte Naherung vom wahren Feldver- 
lauf ab. 


Zweiter Schritt: Einzeichnen von Potentiallinien 
(Bild 21) gem&B Bild 20 als mittlere stetige Kurven der 
V,,-lreppenfunktionen aus Bild 19. Die Potentiallinien stehen 
normal auf den Feldlinien und dem Rand. 


Dritter Schritt: Einzeichnen von treppenférmigen 
B-Linien (Bild 21) anhand der Funktion 4b = 4b(Al) aus 
Bild 18. 


Vierter Schritt: Ausgleichen der B-Treppenfunktion 
analog zum zweiten Schritt. 


Weitere Schritte: Wiederholung der Schritte, 1, 2, 
3, 4, wobei die Spriinge in den Treppenfunktionen fiir V,, 
und B immer kleiner werden. Das iterative Verfahren wird 
abgebrochen, sobald die Treppenspriinge in die GréBenordnung 
der Bleistift-Strichdicke gelangen. 


Bild 19. Angenommene Feldlinien (nullter 
ZweckmaBig wird fiir das Aufzeichnen der fortlaufenden —"*rationsschritt)_ und VTreppenfunktion 


rster Schritt). 
Schritte Null und Eins, Zwei und Drei, usw. Pauspapier be- oe a 
nutzt, so daB das Einzeichnen der mittleren stetigen Kurven aus den Treppenfunktionen ein- 
fach wird. Bild 22 zeigt den Feldverlauf nach Abbruch des Verfahrens. Die Durchflutung 
betragt OG age IG = 2753-4 iir 
das halbe Profil bei einem totalen 
Flu8 von 


oo Midas ea ees 
P= Fy Ke = 5 102 Ns: 


Der Feldverlauf kann auch mittels 
eines zweiten Iterationsverfahrens 
bestimmt werden. Dazu wird der zu 
untersuchende Bereich durch Feld- 
linien und Potentiallinien in plau- 
sibler Weise in Kraftroéhrenab- 
schnitte aufgeteilt. Fiir jeden Ab- 
schnitt wird das Verhaltnis A1/Ab 


Bild 20. Bestimmung der resultierenden stetigen Bild 21. Mittlere stetige Potentiallinien gemaB Bild 22. Feldverlauf im 
Potentiallinie aus der Treppenfunktion fiir Vy. Bild 20 (zweiter Iterationsschritt) aus der Vyy,- Winkelprofil fiir ® = 5-+-10-5 Vs 
Treppenfunktion von Bild 19 und treppenfor- und © = 27,3 A. 


mige B-Linien (dritter Schritt). 


ausgerechnet und anhand von Bild 18 festgestellt, ob der betreffende Abschnitt breiter oder 
kiirzer bzw. schmdler oder langer gemacht werden muB. Dieses Verfahren wird in einzelnen 
Schritten solange durchgefiihrt, bis der Zusammenhang Al = Al(Ab) in jedem Punkt erfiillt ist. 
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4. SchluBbemerkungen 


Das beschriebene graphische Verfahren ist wohl sehr rasch, aber nicht tibermaBig genau. 
Sofern die Parameter K, und K,, d.h. die angenommenen Abmessungen Al, und 4d, dem 
Problem angepaBt sind, kann der Fehler leicht innerhalb der Grenze von 5% gehalten werden. 
Die durch die Annahme von Ab, und Al, festgelegte GréBe der Kraftréhrenabschnitte ist nach 
oben durch die erforderliche Genauigkeit, nach unten durch den Zeitaufwand fiir die Aus- 
wertung begrenzt. Starke Nichtlinearitat und/oder komplizierte Rander bedingen dabei kleine 
Abschnitte, so daB einige Mehrarbeit erforderlich ist. Weiterhin kénnen starke Nichtlinearita- 
ten das Iterationsverfahren verzégern, da bei den einzelnen Schritten iiberkorrigiert wird. 
Anderseits korrigieren sich Fehler im Lauf des Verfahrens von selbst, wie etwa falsches Ablesen 
von Al oder Ad aus der Kurve Al(Ab). SchlieBlich sei noch erwahnt, daB diese fiir die Gl. (A5) 


AV. + ; grad werad V0 


des magnetischen Skalarpotentials V,, gefundene graphische Lésung nicht allgemein auf die 
Behandlung ‘der Porssonschen Gleichung 


= 6 
erweitert werden kann. AP + {(P) 2 (C6) 
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technik der Eidgenéssischen Technischen Hochschule, fiir die Anregung zu vorliegender Arbeit. » 
Diese wurde erméglicht durch Mittelaus dem Jubilaumsfonds und dem Centenarfonds der ETH so- 
wie durch einen Kredit des Aluminiumfonds Neuhausen. Allen diesen Stellen sei herzlich gedankt. 


Zusammenfassung 

Die vorliegende Arbeit beschreibt zwei Verfahren zur Bestimmung von magnetischen 
Feldern im Innern von ferromagnetischen Kérpern bei variabler Permeabilitat ~ = u(#). 
Nach einleitender Ableitung der Differentialgleichungen fiir das magnetische Potential folgt 
die Beschreibung einer Analogiemethode mittels eines elektrischen Widerstandsnetzes. Fiir 
dieses werden die GréBen der Netzwerkelemente berechnet und die praktische Anwendung 
anhand einiger Beispiele erlautert. Der zweite Teil der Arbeit ist eine Erweiterung des bekann- 
ten graphischen Verfahrens von LEHMANN auf Felder mit 6rtlich variabler Permeabilitat bei 
beliebigen Randern. Nach Aufteilen des untersuchten Feldgebietes in einzelne Kraftréhren- 
abschnitte kann zwischen der mittlern Breite und der mittlern Lange eines Abschnittes anhand 
der Magnetisierungskurve B(H) ein funktionaler Zusammenhang gefunden werden, mit dessen 
Hilfe durch iteratives Vorgehen das magnetische Feld bestimmt wird. Das Verfahren wird 
auf einen radial durchflossenen Kreisring und auf ein 90°-Winkelprofil angewandt. 


Literatur 


[1a] LizpMann, G.: Resistance-Network Analogues. Proc. Internat. Comput. Meet. Brussels (1955) 

S. 346—369. — [1b] Lrepmann, G.: Resistance-Network Analogues. Rev. E. Soc. Belg. Electr. Bd. 1 (1956) - 
Nr. 12, S. 2874f. — [1c] Gremognik, G., A. FREI u. M. J. O. Strutt: Neue Anwendungen von Impedanznetz- 
werken als Analogierechengerate. Scientia Electrica Bd. 3 (1957) S. 37—53. — [2] Brecuna, H.: EinfluB von 
Form und Bauart der Transformator-Eisenkérper auf die Ummagnetisierungsverluste. Scientia Electrica Bd. 5 
(1959) S. 1—18. — [3a] LiEBMANN, G.: Resistance-network analogues with unequal meshes or subdivided 
meshes. Brit. J. Appl. Phys. Bd. 5 (1954) S. 362—366. — [3b] Oremo&nik, G.: Ermittlung von Raumladungs- 
feldern und Kapazitaten mit Hilfe eines Widerstandsnetzwerkes mit Stromquellen. Diss. Eidg. Techn. Hoch- 
schule, (1955), Prom. Nr. 2667. — [3c] Cremognix, G., u. M. J. O. Strutt: Praktische Anwendungen eines 
Widerstandsnetzes zur Bestimmung eines ebenen Potential- oder Raumladungsfeldes. Arch. Elektrotechn. 
Bd. 43 (1957) S. 177186. — [4] Fisuer, M. E.: Stability and convergence limitations on the use of analogue 
computors with resistance network analogues. Brit. J. Appl. Phys. Bd. 9 (1958) S. 288—291. — [5a] Lrn- 
MANN, Tu.: Graphische Methode zur Bestimmung des Kraftlinienverlaufes in Luft. ETZ Bd. 30 (1909) S. 995 
bis 998, 1019—1022. — [5b] Lenmann, Tu.: Détermination graphique des champs magnétiques laplaciens et 
tourbillonnaires lignes de flux planes. Rev. gén. Electricité Bd. XIV (1923) S. 347—357, 395—403, 530. — 
[5c] Lrnmann, Tu.: La détermination des champs magnétiques partiels et resultants dans les dynamos sa- 
turées. Rev. gén. Electricité Bd. XXII (1927) S. 1105—1119. — [6] BENEDIKT, O.: Eine neue Berechnungs- 
methode fiir komplizierte, stark gesattigte magnetische Kreise. Acta Technica Academiae Scientiarum Hun- 
garicae Bd. 19 (1957) S. 169—191. 

Dipl.ing. P. TscHopp und Dr. sc. techn. A. FREI, 

Institut fiir Héhere Elektrotechnik, Eidgen. Techn. Hochschule, Gloriastr. 35, Ziirich 7. 


XLIV. Band 
Heft 7 — 1959 


Or 


E. Nasser: Die zeitliche Entwicklung positiver Entladungskanile . . . 45 


Die zeitliche Entwicklung positiver Entladungskanidle 
langs der Oberflache von Isolierfilmen 


Von 
EssAmM NASSER 


(Mitteilung aus dem Hochspannungsinstitut der Technischen Universitat Berlin) 


Mit 7 Textabbildungen 
(Eingegangen am 9g. Juli 1959) 
I. Einfihrung 
1.Der Entladungsaufbau in der Luft 


Zur Untersuchung des raumlichen und zeitlichen Entladungsaufbaues im Feldraum 
einer in atmospharischer Luft befindlichen positiven Spitze wurde in jiingster Zeit eine neue 
Untersuchungsmethode verwendet: Ein photographischer Film wurde zwischen den Elektro- 
den einer Spitze-Platte-Funkenstrecke parallel zur Plattenelektrode angeordnet. Bild 1 
zeigt schematisch die Versuchsanordnung. Beim Anlegen einer geniigend hohen positiven 
StoBspannung wachsen von der Spitze aus Entladungskanile vor, die beim Auftreffen auf den 
Film in dessen Ebene umgelenkt werden. Sie hinterlassen auf ihm ein latentes Bild, das nach 
der Entwicklung des Filmes Auftreffpunkte und davon ausgehende Gleitentladungsspuren auf- 
weist. Bild 2 zeigt eine schematische Darstellung einer solchen Entladungsaufnahme und die Me- 
thode ihrer Auswertung. Durch Verschieben des Filmesin axialer Richtung wurden Entladungsbil- 
der im ganzen Feldraum bei SpannungsstéBen verschiedener Héhe und Zeitdauer aufgenommen. 

Diese Untersuchungen erbrachten neue grundsatzliche Auffassungen itiber den Ent- 
ladungsaufbau im inhomogenen Feld. Uberraschend war zunachst die Feststellung der er- 
heblichen raumlichen Ausbreitung, sowie der erstaunlich groBen Anzahl von Entladungs- 
kanalen, und zwar sowohl in axialer Richtung bis zur Plattenelektrode (Kathode), als auch in 
radialer Richtung bis zum Plattenrand?. Dieses ausgedehnte Entladungsgebiet wurde bereits 
bei StoBspannungen festgestellt, die viel niedriger als die statische Durchschlagspannung der 
untersuchten Anordnung lagen. Diese Beobachung, daB Entladungskandle beide Elektroden 
uberbriicken, ohne einen Durchschlag hervorzurufen, widerspricht den friiheren Auffassungen. 
Die Untersuchung der verschiedenen Abhangigkeiten und die Deutung der einzelnen Beobach- 
tungen wurden in einem friiheren Aufsatz [1] eingehend behandelt. 

Die Untersuchung der zeitlichen Entwicklung dieses groBen komplizierten Entladungs- 
gebildes, d. h. die Untersuchung der Frage, wie schnell die Entladung sich entwickelt und aus- 
breitet, erfolgte unter Verwendung steiler Rechteck-SpannungsstéBe veranderlicher Dauer 
[2]. Es zeigte sich ein duBerst friihzeitiger Entladungseinsatz (etwa 3,5 -107%s nach Anlegen 
der StoBspannung) und ein sehr rasches Vorwachsen der Entladungskanale von der Spitze 
zur negativen Plattenelektrode. Vorwachsgeschwindigkeiten von 1—6-10% cm/s wurden 
dabei gemessen. Die gréBte in einem Entladungsbild vorhandene Spurenlange der Gleitent- 
ladungen wurde in der ersten Arbeit so ausgewertet, wie in Bild 2 gezeigt ist, und dann zur 
Klarung der physikalischen Vorgange im Entladungskanal in der atmospharischen Luft heran- 
gezogen [1]. Die zeitliche Entwicklung dieser GréBe wurde in der zweiten Arbeit zwar aus- 
gewertet, auf ihre Deutung und auf die Deutung der sich an der Oberflache des Isolierfilmes 
abspielenden Vorgange wurde jedoch nicht eingegangen, weil die Entwicklung dieser Ioni- 
sierungsvorgange auf der Filmoberflache fiir die Untersuchung der Vorgange im Gasraum 
nicht weiter interessierte. 

In dem vorliegenden Aufsatz sollen nun an Hand des in beiden genannten Arbeiten® ge- 
wonnenen reichen Filmmaterials sowohl die verschiedenen Bildgr6Ben — in Hinblick auf die 
Vorgange auf der Filmoberflache — als auch die Einzelheiten der erhaltenen Entladungs- 
bilder ausgewertet werden, woraus dann die zeitliche Entwicklung der Entladungskanale 

rs Be Loes berichtet in ,,Report of the 3. Intern. Congress on the Ionisation of Gases‘ 1957, Venedig, 
S. 646, daB G. G. Hupson das Auftreffen von Entladungskandlen auf der Kathodenplatte erst bei Spaunungen 


knapp unterhalb der statischen Durchschlagspannung photographieren konnte. 
2 Fiir eine Zusammenfassung der beiden Arbeiten s. ETZ-A Bd 80 (1959) H. 19, S. 674. 
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langs der Oberflache des Filmes sich ermitteln 1a8t. Diese Ergebnisse werden dann in Zusam- 
menhang mit densich im Gasraum abspielenden Ionisierungsvorgangen wie7.B. (Ionisierung durch 
ElektronenstoB, Entwicklung und Aufbau des Entladungskanals und der von diesem ausge- 
sendeten ionisierenden Strahlung) gebracht und unter Heranziehung unserer heutigen Vor- 
stellung iiber den raum-zeitlichen Entladungsaufbau in der atmospharischen Luft bei positiver 
Spitze gedeutet. 

2. Die Gleitentladung 


Es sei zunichst eine Ubersicht iiber die dltere Literatur gegeben. Die Entladungen langs 
der Oberflaiche von Isolierstoffen, als Gleitentladungen bezeichnet, wurden nach ihrer Ent- 
deckung durch LICHTENBERG yon verschiedenen Forschern untersucht. Die meisten Ver- 
dffentlichungen, die bisher erschienen, befaBten sich entweder mit der raumlichen Ausbreitung 
der LICHTENBERG-Figuren und ihrer Spannungsabhangigkeit, um die in den elektrischen An- 
lagen auftretenden Uberspannungen auf einfache Weise registrieren und auswerten zu kénnen, 
was zu der Entwicklung des Klydnographen fiihrte — oder aber mit der zeitlichen Entwick- 
lung der LICHTENBERG-Figuren, um Aufschliisse itber den Mechanismus der Gleitentladung zu 
erhalten, was insofern technisch bedeutsam ist, als diese Entladungen langs der Oberflache 
eines Dielektrikums die elektrische Festigkeit der betreffenden Isolieranordnung stark herab- 
setzen. Zusammenfassend berichten dariiber GANGER [3], STRIGEL [4], MEEK und Craccs [5]. 

Der erste Versuch, die Gleitentladungen bzw. die LICHTENBERG-Figuren zum Studium der 
Entladungs- und Jonisierungsvorgange in Gasen zu verwenden, wurde von MERRILL und VON 
HiPpEL unternommen [6]. Die Deutung dieser in verschiedenen Gasen bei unterschiedlichem 
Druck gewonnenen Gleitentladungsfiguren und ihre Interpretation zur Aufklaérung der Ioni- 
slerungsvorgange in der Luft miissen jedoch heute wegen des damals noch unbekannten 
Kanalentladungsaufbaus als unzulanglich angesehen werden. Diese alteren Auffassungen sollen 
deshalb im folgenden an Hand von neuen, systematisch aufgenommenen photographischen 
Aufzeichnungen von Gleitentladungen, die sich von einer positiven Spitze in atmospharischer 
Luft aus entwickeln, gepriift und gegebenenfalls berichtigt oder erweitert werden. 

Die Vorwachsgeschwindigkeit positiver Gleitentladungskanale von einigen Millimetern 
Lange wurde bereits von RoGowski und Mitarbeitern durch Messung der Wirkdauer abge- 
schnittener StoBwellen mit etwa 2-107 cm/s bestimmt [7]. Dabei wurde die ToEPLERsche Gleit- 
anordnung (Spitze senkrecht auf einer auf die Plattenelektrode gelegten Isolierplatte) ver- 
wendet. Dracu ermittelte mit einer ahnlichen Anordnung noch héhere Geschwindigkeiten 
von maximal 1- 10% cm/s [8]. FUNFNER fand mit der von ToEPLER angegebenen Schlieren- 
methode, bei der die in einer Fliissigkeit sich ausbreitende, von der Schaltfunkenstrecke des 
StoBkreises ausgeléste Schallwelle mit dem Licht des, kurzschlieBenden Uberschlagfunkens be- 
strahlt wird, Vorwachsgeschwindigkeiten zwischen 0,5 und 3,0 - 107 cm/s bei 15 bzw. 30 kV [og]. 

PEDERSEN verwendete zwei langliche, auf der Photoschicht einer lichtempfindlichen 
Platte nebeneinander angeordnete Elektroden [10], die iiber unterschiedlich lange Leitungen 
mit dem StoBkreis verbunden waren. Durch Betrachtung der Begrenzungslinie der beiden 
ineinandergreifenden Figuren bei verandertem Leitungsunterschied der Zuleitung schloB 
PEDERSEN auf ein Ausbreitungsgesetz der Figuren von der Form 

fo Ki (dee), 
wobei 7 die jeweilige Lange der vom Pol ausgehenden Strahlen, R die endgiiltige Strahlen- 
lange der fertig ausgebildeten Figur und ¢ die Zeit ist. Fiir « fand er an positiven Figuren bei 
Atmospharendruck 0,26-10° st. Daraus erhalt man fiir die Ausbildung einer annahernd 
vollstandig entwickelten Figur eine Mindestzeit von etwa 1,2 bis 1,5-1077s. SToERK und 
BUNGARDEAN mafen bei Untersuchungen an Klydnographen dagegen Ausbildungszeiten der 
positiven Gleitentladungsfiguren von nur 4-107*s [11]. 


3. Ziel der Arbeit 
Im folgenden wird das Versuchsmaterial, daB in den eingangs genannten Arbeiten iiber 
den Entladungsaufbau in einer Spitze-Platte-Funkenstrecke gewonnen wurde, im Hinblick 
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auf die Entwicklung von Gleitentladungen ausgewertet und mit den erwahnten Messungen 
der alteren Literatur verglichen werden. Eine Erklarung der vorhandenen Abweichungen 
wird dann mit Hilfe der gewonnenen Anschauungen versucht. Die Deutung der Gleitent- 
ladungsbilder und der zeitlichen Abhangigkeit der in ihnen enthaltenen ausgewerteten GréBen 
diirfte zu einer besseren Vorstellung ttber den Aufbaumechanismus der Gleitentladungskanile 
beitragen. 


II. Versuchseinrichtung 

Zum Verstandnis der Entladungsbilder sei die bei den einleitend genannten Arbeiten [1, 2] 
verwendete Apparatur nochmals kurz beschrieben: Sie wurde mit einem auf WINKELNKEMPER 
zuriickgehenden Versuchsaufbau durchgefiihrt. Wie 
in Bild 1 schematisch dargestellt, befindet sich 
dabei ein photographischer Film zwischen den 
Elektroden einer Spitze—Platte-Anordnung, und 
zwar parallel zur Plattenebene, mit der Schichtseite 
der Spitze zugewendet. Die Entfernung s zwischen 
der Filmebene und der Spitze laBt sich durch Ver- 
schieben des Filmes in axialer Richtung dndern. fmm? 
Bei den im folgenden beschriebenen Messungen a 
wurde durchweg eine halbkugelige Spitze von 1 mm 
Durchmesser und eine gut abgerundete Kupferplatte 
von etwa 100 mm Durchmesser verwendet, wobei die 
Schlagweite bei diesen Versuchen 25 mm _ be- 
trug. 

Zur Erzeugung steiler Rechteck-Spannungsst6Be 
verdnderlicher Dauer diente eine symmetrische 
Wanderwellenanlage ; die wirksame Lange der Doppel- Bild 1. Schematische Darstellung der Versuchsanordnung; 
leitung wurde dabei mit Hilfe eines KurzschluBbiigels ¢ Piestoten, ¢ Pies ton tad 
verdndert. Die Stirnzeit der StoBwelle betrug 
infolge Verwendung einer in komprimierten Stickstoff eingebauten Schaltfunkenstrecke etwa 
8. 107s; 


III. Versuchsergebnisse 
1. Auswertung der Entladungsbilder 


Die Entladungsbilder wurden bei verschiedenen Entfernungen s des Filmes von der Spitze 
jeweils durch Veranderung der Leitungslange (d. h. Veranderung der Impulsdauer) bei kon- 
stanter Spannungshohe gewonnen. Zur Erzeugung jedes Entladungsbildes wurde ein Span- 
nungsstoB8 herangezogen. In Bild 3 und 4 sind solche Entla- 
dungsbilder wiedergegeben. Man erkennt, daB die Figuren dann, 
wenn der Film unmittelbar an der Spitze liegt, den mit der 
ToEPLeERschen Gleitanordnung erhaltenen LICHTENBERG-Figuren 
sehr 4hneln, obwohlin der bei den vorliegenden Versuchen ver- 
wendeten Anordnung statt eines festen Dielektrikums (wie bei 
ToEPLER) Luft als Isolationsschicht zwischen Film und Plattevor- 
handen war. Auf einem in einer gewissen Entfernung vor der 
Spitze angeordneten Film entwickeln sich von verschiedenen Auf- 
treffpunkten ausgehend mehrere Gleitfiguren nebeneinander, die Bidimecsen Gecko Decree 
insbesondere bei StoBspannungen langer Wirkdauer sehr unter- RR ICES. 
schiedliche GréBen haben. Bemerkenswert ist, daB diese Figuren sich nicht kreuzen, iiber- 
lappen oder gar beriihren. 

Die Auswertung dieser Aufnahmen wurde auf folgende GréBen beschrankt (s. Bild 2): 

a) die gr6Bte in einem Bild vorhandene Gleitentladungsspur /,, kurz Gleitspurenlange genannt, 

b) den Radius } des Kreises, der alle Auftreffpunkte umfaBt und 

c) die Anzahl der Auftreffpunkte. 
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Die letzten beiden GréBen waren — wie eingangs bereits erwahnt — in friiheren Arbeiten 
in ihrer raumlichen und zeitlichen Abhangigkeit wesentlich, wahrend die erste GroBe — also 
die Gleitspurenlange — fiir die genannten Untersuchungen unwesentlich war. Sie soll nun 
naher ausgewertet und in Abhangigkeit von der Zeit dargestellt werden. 


2. Durchfithrung der Versuche 


Wichtig fiir das Verstandnis der entstandenen Entladungsbilder waren auch noch die 
folgenden Gesichtspunkte bei der Durchfithrung der Versuche. Da die Vorversuche nicht nur 
Bilder groBer Streuung ergaben, sondern manchmal sogar keine Schwarzung auf den Filmen 
hinterlieBen, wurde die filmfreie Versuchsstrecke vor jeder Versuchsreihe durch Anlegen 
hoherer Uberspannungen mehrmals zum Durchschlag gebracht. Auf diese Weise konnten bei 
gleichen Bedingungen Entladungsbilder gewonnen werden, die nur geringfiigige Unterschiede 
aufwiesen. Dies kann darauf zuriickzufiihren sein, daB die Wahrscheinlichkeit des Vorhanden- 
seins eines Elektrons im Gebiet hoher Feldstarke in Spitzennahe wahrend dieses kurzen Zeit- 
raumes sehr gering ist, was auch durch die von STRIGEL [4] gemessenen, relativ langen Streu- 
zeiten von tuber 5-107"s bei einer mit Pappschirm abgedunkelten Anordnung erklart wird. 
Durch die jeder Versuchsreihe vorangehende hoch ionisierte Gasentladung, die auch extrem 
kurzwelliges Ultraviolett ausstrahlt, wurden im Entladungsraum und in seiner naheren Um- 
gebung Ladungstrager in groBer Anzahl gebildet. Dadurch wurde die bis zum Beginn des 
Entladungsaufbaues verstreichende Streuzeit — ebenso wie durch Bestrahlung einer Funken- 
strecke mit Ultraviolett — erheblich verringert. Das Aussetzen der Entladung wurde, solange 
die Zeit zwischen den Uberschlagen und den eigentlichen Versuchen nicht zu lang war, dadurch 
fast véllig vermieden. 

3. MeBergebnisse 

Bei konstanter Elektrodenspannung wurde dann die Gleitspurenlange (Lange des Gleit- 
kanals) 1, aus den bei verschiedener Zeitdauer der Spannungswelle erhaltenen Entladungs- 
bildern gemessen, und zwar bei jeweils kostantem Abstand zwischen Film und Spitze. 


t= 4 10 40-10-8 s 
Bild 3. Gleitentladungsfiguren bei verschiedenen Impulszeiten ¢ (4, 10 und 40+10-*s); Film unmittelbar an der 
Spitzenelektrode (s = 0); StoBspannung 38,5 kV. 


In Bild 3 und 4 sind derartige Aufnahmereihen wiedergegeben. Bild 3 zeigt die Gleitent- 
ladungsfiguren bei einer Spannungsdauer von 4, 10 und 4o- 10~°s und einer Spannungshéhe 
von 38,5 kV, wobei der Film unmittelbar an der Spitze anlag. Man erkennt deutlich, daB die 
Gleitspurenlange mit zunehmender Zeit wachst und daB die Anzahl der vom Gleitpol ausgehen- 
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den Entladungskanale von der Zeit unabhangig ist. Das bedeutet, daB die verschiedenen 
Aste nahezu gleichzeitig entstehen und vorwachsen. In Bild 4 sind die Entladungsaufnahmen 
bei s = 20 mm, einer Sto8spannung von 38,5 kV und einer Wirkzeit von 12,20 und BOwa. oS 
wiedergegeben. Es ist ersichtlich, wie sich die Gleitentladungen von den verschiedenen Auf- 
treffpunkten aus entwickeln. Wenn sie sich einander nahern, stoBen sie sich gegenseitig ab, 
so da’ manche Gleitkanale abgelenkt werden, was besonders in Bild 4 deutlich zu sehen ist. 


pe 20 30-107 8s 
Bild 4. Entladungsfiguren bei verschiedenen Impulszeiten ¢ (12, 20 und 30+ 10~-* s); Filmentfernung von der Spitze 
Ss = 20 mm; StoBspannung 38,5 kV. 


Die Beziehung zwischen J, und der Zeit ¢ (gemessen vom Anlegen der StoBspannung) an 
verschiedenen Stellen des Filmes bei unterschiedlichen Spannungen ist in Bild 5 graphisch 
dargestellt. Aus den Kurven /, = f(t), die alle einen ahnlichen Verlauf haben, kann man er- 
kennen, da der Gleitentladungskanal erst eine gewisse Zeit ¢, nach Anlegen der StoBspannung 
vorzuwachsen beginnt. Diese Zeit ¢, ist naturgemaB von der jeweiligen Entfernung s des Filmes 
von der Spitze abhangig. Wenn der Film unmittelbar an der Spitze anliegt, ist diese Zeit- 
spanne ¢) bei beiden Spannungen nahezu gleich. 

Die Zeit ¢, setzt sich aus den Zeiten ¢) und ¢, zusammen. Diese beiden Zeiten sind physi- 
kalisch folgendermaBen zu verstehen: f, ist die Zeitspanne, die vom Anlegen der StoBspannung 
bis zum Beginn des Vorwachsens einer Kanalentladung verstreicht. ¢)ist in erster Linie durch den 
Anstieg der Spannungswelle bedingt, weniger durch deren maximale Hohe. Wie in einer friiheren 
Arbeit theoretisch und rechnerisch festgestellt wurde [2], betragt die bei der verwendeten Welle 
biszum Erreichen der Einsatzspannung verstreichende Zeit etwa wenige 10-*s. Zu diesem Zeit- 
punkt ist die Feldstarke in der Nahe der Spitze so stark angewachsen, daB sich eine gestartete 
Elektronenlawine bis zur kritischen Vermehrung verstarken kann, wonach dann die Kanalentla- 
dung rasch vorwachst. Die Laufzeit dieser ersten, ,ziindenden“ Elektronenlawine liegt in der Gré- 
Benordnung von einigen 101° s, je nachdem, wo sich das erste Elektron gerade befand. Die 
Zeit t, ist die Laufzeit der Kanalentladung von der Spitze (s = 0) bis zu einer Entfernung s von 
der Spitze. ¢, ist damit die Zeit, die vom Anlegen der StoBspannung bis zu dem Augenblick ver- 
streicht, in dem die Kanalentladung den in einer Entfernung s von der Spitze angeordneten 
Film erreicht hat. Auf Bild 5 sind ¢, und ¢, fiir den Falls = 10 mm mit ¢,, und #,, bezeichnet. 

Aus Bild 5 kann man auch erkennen, da8 alle Kurven fast den gleichen Verlauf haben; sie 
sind im ersten Bereich sehr steil und werden mit zunehmender Zeit stetig flacher, bis sie sich 
dann 4—5- 1078s nach Anlegen der StoBwelle einem Endwert nahern. 

Durch Differenzierung der Lange—Zeit-Kennlinien laBt sich die Vorwachsgeschwindigkeit 
v, der Gleitentladungskandle bei verschiedener Entfernung s des Filmes von der Anode in 
Abhangigkeit von der Zeit f entnehmen. Bild 6 zeigt die Beziehung v, = /(¢) mit s als Para- 
meter bei den Spannungshéhen von 25,0 und 38,5 kV. Aus diesen Kurven erkennt man, 
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wie die Vorwachsgeschwindigkeit des Gleitkanals von 4—5 - 10° cm/s auf etwa 0,1: 108 cm/s 
abnimmt, und zwar zuerst schnell und dann immer langsamer. 

Da aber die Gleitkandle je nach der Entfernung des Filmes von den Elektroden nach ver- 
schiedener Zeit vorzuwachsen beginnen, ist es zweckmaBig, den Verlauf der Vorwachsgeschwin- 
digkeit des Gleitkanals in Abhangigkeit von seiner Lange aufzuzeichnen. Diese Kurven wurden 
aus den obigen Beziehungen ermittelt und sind in Bild 7 dargestellt. Hiernach hat der Kanal, 
wenn der Film unmittelbar an der Spitze liegt, eine Anfangsgeschwindigkeit von etwa 6 - LO” 
cm/s bei 38,5 kV bzw. 5-10%cm/s bei 25,0kV. Die Kurven der Vorwachsgeschwindigkeit 
der Gleitkanale bei Entfernungen s = 10 und 20 mm liegen unterhalb der Kurven bei s = 0. 
Das bedeutet, daB bei gleicher Kanallange die Vorwachsgeschwindigkeit um so niedriger ist, 
je weiter der Film von der Spitzenanode entfernt ist. 
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Bild 5. Gleitkanallinge /, abhangig von der Impulszeit ¢ bei verschiedenen Spannungen 
und Entfernungen s des Filmes von der Spitze. 


Bemerkenswert ist die Feststellung, daB bei der niedrigen Spannung die Gleitkandle nach 
etwa 4+ 10~*s noch ganz langsam vorwachsen (mit einer Geschwindigkeit von etwa 0,1 - 108 
cm/s), bis sie ihre gesamte Ausbreitung erreicht haben. Da die Kandle bei 38,5 kV in ihrem 
Endstadium sehr schnell zum Stillstand kommen, kann darauf zuriickzufiihren sein, daB sie 
dann, wenn sie eine Lange von etwa 45 mm haben, also ein der Plattengr6Be entsprechendes 
Gebiet bedecken, nicht mehr unter den gleichen Bedingungen weiter vorwachsen kénnen, 
da das elektrostatische Feld in diesem auBeren Bereich eine der Richtung der Vorwachs- 
geschwindigkeit entgegengesetzte Komponente besitzt, die die Entladung rasch abbremst. 
Bei Verwendung einer unendlich groBen Plattenelektrode wiirden sich als Beziehung zwischen 
v, und J, bei l, > 35 mm etwa die gestrichelten Kurven in Bild 7 ergeben, was dem Verlauf 
der Kurven bei 25 kV entsprache. 

Aus den Kurven Bild 6 und 7 ist ersichtlich, daB die Vorwachsgeschwindigkeit der Gleit- 
kanale von einigen 108 cm/s im Anfang auf etwa ein Zehntel dieses Wertes im Endstadium 
zurtickgeht. Die von verschiedenen Forschern angegebenen unterschiedlichen MeBwerte 
k6nnen sich daher auf verschiedene Entwicklungsphasen beziehen und brauchen insoweit 
nicht widerspriichlich zu sein. 

Unter der in der Literatur angegebenen Ausbildungszeit einer LICHTENBERG-Figur ist die 
Zeit zu verstehen, die die Figuren bei einer Vollwelle fiir ihre gesamte Entwicklung brauchen. 
Da aber, wie oben erladutert, bei der an sich anzustrebenden Verwendung groBer Platten- 
elektroden die letzten Entwicklungsphasen sehr langsam gegeniiber den ersten vor sich gehen, 
empfiehlt es sich, die Ausbildungszeit einer Gleitentladungsfigur als die Zeit zu definieren, 
die zur Bildung einer Figurenlange von go0% der endgiiltigen, bei langer Beanspruchungs- 
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dauer erhaltenen FigurengréBe erforderlich ist. Demnach betragt die Ausbildungszeit T, 
der Gleitfiguren bei 38,5 kV etwa 30, 35 und 37 - 10~° s beis = 0, 10 und 20 mm, wie aus Bild 5 
hervorgeht. Bei 25 kV ist 7), (d.h. T, bei s = 0) = 3-107 s; bei gréBeren Entfernungen 
des Filmes liegt 7, zwischen 10 und 20 - 19-8 s, weil hier die Gleitentladung nur ganz langsam 
zum Stillstand kommt. 
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Bild 6. Vorwachsgeschwindigkeit vg des Gleitkanals abhangig von der 0 10 20 30 40 mm 60 
Impulszeit ¢ bei verschiedenen Spannungen und Entfernungen des Gleitkanallange ly 


Filmes von der Spitze. 
Bild 7. Vorwachsgeschwindigkeit vg des Gleitkanals abhangig von 
der Gleitkanallange by bei verschiedenen Spannungen und Entfer- 
nungen des Filmes von der Spitze. 


IV. Deutung der MeBergebnisse 

Der Mechanismus der Gleitentladung wurde bereits mehrfach behandelt, wobei die zu- 
grundegelegte Entwicklungsgeschwindigkeit fast um eine Zehnerpotenz niedriger als die in 
diesem Aufsatz mitgeteilte war. Grundsatzlich wurde jetzt der Nachweis erbracht, da8 die 
Vorwachsgeschwindigkeit der positiven Kanale im Feldraum und auf der Oberflaéche eines 
Dielektrikums von derselben GréBenordnung ist. (Bei positiver Spitze wurden im Feldraum 
Vorwachsgeschwindigkeiten von 1—7- 10% cm/s gemessen [1, 2]). Dies legt die Vermutung 
nahe, daB der Vorwachsmechanismus der Entladungskandle im Gasraum und auf der Isolator- 
oberflache gleich ist. Das bedeutet, daB die von den Lawinen abgestrahlten ionisierenden 
Photonen an der Entwicklung der Entladungskandale in beiden Fallen sehr maBgeblich be- 
teiligt sein miissen, da sonst derart hohe Vorwachsgeschwindigkeiten nicht zu erklaren sind. 

Den Mechanismus des auBerst raschen Vorwachsen dieser Entladungskanale kann man 
sich nach unserer heutigen Anschauung wohl etwa folgendermaBen vorstellen: Ein Elektron in 
Spitzennahe vermag unter der Einwirkung des zur Spitze hin rasch ansteigenden Feldes durch 
StoBionisation eine Elektronenlawine zu bilden, die eine starke Anhaufung positiver Ionen 
hinterlaBt, so daB das elektrische Feld an einer vorgeschobenen Stelle erheblich verstarkt 
wird. Die Photoelektronen, die durch die vom Lawinenkopf nach allen Richtungen abgestrahl- 
ten Photonen gebildet werden, werden durch das verstarkte Feld beschleunigt und erzeugen 
bei ihrem Lauf zur Anode Lawinen, die wiederum Anhaufungen positiver Ionen hinterlassen 
und Photonen in allen Richtungen aussenden. Dadurch wachst der Gleitkanal solange weiter 
vor, bis das resultierende Feld nicht mehr zur Ionisation durch ElektronenstoB ausreicht. 

Die Bedingungen fiir das Vorwachsen eines Entladungskanals sind demnach folgende: 
1. Eine ausreichende Erhéhung der Feldstarke durch die zuriickgebliebenen positiven Ladungs- 
traiger und 2. eine gentigende Anzahl ionisierende Photonen mit einer passenden mittleren freien 
Weglange. Diese darf einerseits die Reichweite der erhéhten Feldstarke nicht tiberschreiten, 
weil die gebildeten Photoelektronen dann nicht stoBionisieren kénnen. Wenn andererseits 
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die mittlere freie Weglange der abgestrahlten Photonen zu klein ist, kénnen die Photoelektronen 
zwar durch das verstarkte Feld stoBionisierend Lawinen bilden, die jedoch wegen der zu geringen 
Entfernung zum Kanalkopf die kritische Vermehrung nicht erreichen. Die mittlere freie Weglange 
der Photonen ist lediglich vom Gas und seinem Druck abhangig. Deshalb ist die Entwicklung 
von Entladungskandalen je nach Gasart und Druck sehr verschieden und fithrt zu den unter- 
schiedlichen Gleitentladungsfiguren, wie sie MERRILL und voN HippEr bei verschiedenen 
Gasen und unterschiedlichem Druck erhielten [6]. 

Da diese Vorwachsbedingungen sowohl in der Hauptausbreitungsrichtung des Kanals 
als auch in seitlicher Richtung erfiillt sein knnen, verzweigen sich die Kandle standig wahrend 
ihres Vorwachsens. Dies fiihrt zu der beobachteten Entladungsstruktur und bedeutet, daB 
die zahlreichen Zweige einer positiven Gleitfigur nicht zu verschiedenen Zeiten entstanden 
sind, wie etwa GANGER [3] nach ToEPLER annimmt, sondern sich nahezu gleichzeitig ent- 
wickelt haben. Die Gleitfiguren in Bild 3 mégen diese Feststellung veranschaulichen. Man 
erkennt, daB die Verastelung von Anfang an mit dem Vorwachsen in der Hauptrichtung 
Schritt halt. 

Im Entladungskanal wandern die gebildeten Elektronen rasch zur Anode, und zwar — 
wegen der niedrigen Feldstarke — ohne dabei weitere Ladungstrager zu erzeugen. Die Kanale 
sind also in diesem Stadium keine hoch ionisierten Plasmaschlauche, sondern besitzen eine 
iiberschiissige positive Ladung und stoBen sich wahrscheinlich deshalb gegenseitig ab, wie es 
in Bild 4 deutlich zu erkennen ist. 


V. Zusammenfassung 


An Hand von Versuchsmaterial, das in friiheren Arbeiten zur Untersuchung des raum- 
zeitlichen Entladungsaufbaues diente, wurden die Lange—Zeit-Kennlinien der positiven Gleit- 
entladungskanale aus photographischen Aufzeichnungen bestimmt. Hierbei wurde eine po- 
sitive Spitze/negative Platte-Anordnung verwendet, wobei Photopapier senkrecht zur Haupt- 
feldrichtung an beliebige Stellen des Feldraumes gebracht wurde. Wegen der im letzten Ent- 
wicklungsstadium ganz langsam wachsenden Figur wurde als Ausbildungszeit der Gleitfigur 
die zur Bildung von 90% der vollstandig entwickelten Figur erforderliche Zeit gewahlt. Diese 
Zeit wurde zu 2—20- 108s bestimmt. 

Aus diesen Kennlinien, die bei verschiedener Entfernung des Filmes und unterschiedlicher 
Spannungshohe ermittelt wurden, wurde die Vorwachsgeschwindigkeit der Gleitentladungs- 
kandle in allen Entwicklungsphasen ermittelt. Diese liegt zwischen 1 und 6 - 108 cm/s, d. h. um 
fast eine Zehnerpotenz héher als die bisher gemessene Vorwachsgeschwindigkeit der sich im 
Gasraum entwickelnden positiven Entladungskanale. Deshalb konnte auf einen gleichen 
Aufbaumechanismus geschlossen werden, der dann in Anlehnung an die Vorstellung tiber den 
Kanalaufbau im Gasraum dargelegt wurde. 

Fir die weitgehende Unterstiitzung wahrend der Durchfiihrung der Arbeiten bin ich Herrn 
Prof. R. StRIGEL und Herrn Dr. H. WINKELNKEMPER zu besonderem Dank verpflichtet. 
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ZUR INFORMATION 


Das Buch richtet sich an Studierende der Ingenieurwissenschaften und der Physik und gibt eine Ein- 
fiihrung in die Mathematikin einem Umfang, der etwa den beiden ersten Semestern der mathematischen 
Kursvorlesungen an den Technischen Hochschulen entspricht (Differential- und Integralrechnung samt 
Grundziigen der Vektorrechnung und der analytischen und Differentialgeometrie). Es wird durchwegs 
versucht, die Begriffe und Methoden in einer Weise zu entwickeln, die der aufs Anschauliche gerichteten 
Denkweise des Ingenieurs angemessen ist. Numerische und graphische Methoden nehmen einen verhialt- 
nismaBig breiten Raum ein. Auch eine Einfiithrung in die Nomographie ist eingefiigt sowie ein Exkurs 
in die Getriebelehre. 


Ein spater folgender zweiter Band wird sich mit den Differentialgleichungen und der Funktionentheorie 
befassen. 
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